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Abstract 



Cross sections for suporsymmctric particles production at hadron colliders have been extensively 
studied in the past at leading order and also at next-to-leading order of perturbative QCD. The 
radiative corrections include large logarithms which have to be resummed to all orders in the strong 
coupling constant in order to get reliable perturbative results. In this work, we perform a first and 
extensive study of the resummation effects for supersymmetric particle pair production at hadron 
colliders. We focus on Drell-Yan like slepton-pair and slepton-sneutrino associated production in 
minimal supergravity and gauge-mediated supersymmetry-breaJiing scenarios, and present accurate 
transverse-momentum and invariant-mass distributions, as well as total cross sections. 

In non-minimal supersymmetric models, novel effects of flavour violation may occur. In this case, 
the flavour structure in the squark sector cannot be directly deduced from the trilinear Yukawa 
couplings of the fermion and Higgs supermultiplets. We perform a precise numerical analysis of 
the experimentally allowed parameter space in the case of minimal supergravity scenarios with 
non-minimal flavour violation, looking for regions allowed by low-energy, electrowcak precision, 
and cosniological data. Leading order cross sections for the production of squarks and gauginos at 
hadron colliders are implemented in a flexible computer program, allowing us to study in detail the 
dependence of these cross sections on flavour violation. 



Les sections efficaces de production hadronique de particules supersymctriques out ete largement 
etudiees par le passe, aussi bien a I'ordre dominant qu'a I'ordre sous-dominant en QCD pertur- 
bative. Les corrections radiatives incluent de larges termes logarithmiques qu'il faut resommer a 
tous les ordres afin d'obtenir des predictions consistantes. Dans ce travail, nous effectuons une 
premiere etude dotaillce des effets de resommation pour la production hadronique de particules 
supersymctriques. Nous nous concentrons sur la production de type Drell-Yan de sleptons et sur la 
production associee d'un slepton et d'un sneutrino dans des scenarios de supergravite minimale et 
de brisure de supersymetrie vehiculee par interactions de jauge, et nous presentons des distributions 
d'impulsion transverse et de masse invariante, ainsi que des sections efficaces totales. 

Dans les modeles supersymetriques non minimaux, de nouveaux effets de violation de la saveur 
peuvent avoir lieu. Dans ce cas, la structure de saveur dans le secteur des squarks ne pent pas etre 
deduite directement du couplage trilineaire entre les supermultiplets de Higgs et de fermions. Nous 
effectuons une analyse numerique de I'espace des parametres permis dans le cas de scenarios de 
supergravite minimale avec violation de la saveur non minimale, cherchant les regions permises pax 
les mesures de precision electrofaibles, les observables a basse energie et les donnees cosmologiques. 
La dependance des sections efiicaces a I'ordre dominant pour la production hadronique de squarks 
et de jauginos par rapport a la violation de la saveur non minimale est etudiee en details. 
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Resume 



Le Modele Standard (SM) de la physique des particules [IJ [H [3l H O El E] 
decrit avec succes un grand nombre de donnees experimentales de haute energie. 
Cependant, certaines questions fondamentales restent sans reponse, comme par 
exemple les origines de la brisure de la symetrie electrofaible et des masses des 
particules, la large hierarchie entre I'echelle de Planck et I'echelle electrofaible, 
le mecanisme responsable des oscillations de neutrinos, les origines de la matiere 
sombre et de la constante cosmologique, ou encore le probleme CP lie a I'interaction 
forte. Les tentatives visant a relier differents parametres du SM menent en general 
a des theories plus fondamentales qui resolvent naturellement certains de ces 
problemes ouverts. 

La philosophic generale des Theories de Grande Unification (GUTs) [8l[9l[10] est 
de considerer que les groupes de symetrie du SM emergent de la brisure d'un groupe 
simple de rang plus eleve. A I'echelle GUT, les trois constantes de couplage de jauge 
du SM sont unifiees, et les quarks et les leptons sont decrits par des representations 
communes de ce groupe de jauge plus large. Ces theories predisent en general 
un certain nombre de bosons de jauge additionnels, menant eventuellement a des 
interactions pouvant violer la conservation des nombres baryonique et leptonique. 
II s'agit de I'un des problemes phenomenologiques les plus importants pour les 
GUTs, vu que les baryons sont alors instables, ce qui est contraire aux donnees 
experimentales liees a la non observation de la desintegration du proton. Les theories 
GUTs peuvent expliquer la quantification de la charge electrique et incorporer des 
neutrinos massifs, mais out quelques difficultes pour reproduire la valeur mesuree 
de Tangle de melange electrofaible. De plus, un grand nombre de problemes 
conceptuels deja presents dans le SM demeurent sans reponse. 

Une approche populaire pour resoudre le probleme de hierarchie du SM est 
d'ajouter a I'espace-temps des dimensions supplementaires [HI [12]. Dans ce cadre 
theorique, les interactions de jauge et gravitationnelle sont unies a une echelle 
proche de I'echelle electrofaible, qui est alors la seule echelle fondamentale de la 
theorie, la valeur importante de I'echelle de Planck etant seulement une consequence 
de la presence des nouvelles dimensions. L'espace a quatre dimensions habituel 
est contenu dans une "brane" quadridimensionnelle, elle-meme incluse dans une 
structure plus large contenant N dimensions additionnelles, le "bulk". Dans ces 
theories, chaque champ du SM possede une serie d'excitations de Kaluza-Klein avec 
les memes nombres quantiques, mais une masse differente. Au jour d'aujourd'hui, 
ces excitations n'ont pas encore ete observees, mais I'ordre de grandeur de leur 
masse est le TeV, ce qui les rend tout a fait detectables au futur Grand Collisionneur 
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de Hadrons, le LHC, au CERN. 



Plus recemment, d'autres tentatives pour resoudre ce probleme de hierarchie 
ont ete proposees, comme par exemple les theories "Little-Higgs" ou "Twin-Higgs" , 
qui predisent egalement de nouvelles particules avec des masses de I'ordre du TeV 
[13l[lll[15]. Ces theories incluent des partenaires fermioniques pour les quarks et les 
leptons du SM, et des partenaires bosoniques pour les bosons de jauge. Cela permet 
la stabilisation de la masse du boson de Higgs au-dela de I'ordre dominant grace a 
la realisation d'une symetrie non lineaire reliant les couplages au boson de Higgs 
d'une faqon telle que les divergences venant des corrections quantiques s'annulent. 

Dans cette these, nous nous concentrons sur une autre extension attractive du 
SM, la supersymetrie (SUSY) [El EZl HI [El [20], et plus precisement le Modele 
Standard Supersymetrique Minimal (MSSM) [21], [22j. La supersymetrie a basse 
energie fournit une solution naturelle a plusieurs des problemes conceptuels du 
SM. Reliant les fermions et les bosons, elle permet la stabilisation de la hierarchie 
separant I'echelle de Planck de I'echelle electrofaible \23\ [23] et I'unification des 
couplages de jauge aux hautes energies [2S1 [2B1 [271 [2B]- De plus, la particule SUSY 
la plus legere pent dans certains cas etre vue comme un candidat potentiel pour 
la matiere sombre [Ml [30]. Vu que les partenaires supersymetriques des particules 
du SM n'ont pas encore ete observes jusqu'a present, la supersymetrie doit etre 
brisee a basse energie, mais de fagon douce afin qu'elle reste une solution viable 
pour le probleme de la hierarchie. Les particules SUSY sont done plus massives 
que leurs equivalents du SM, mais leur masse ne devrait pas exceder quelques 
TeV. Une recherche concluante couvrant un large regime de masses allant jusqu'a 
I'echelle du TeV est done I'un des points principaux du programme experimental 
des collisionneurs hadroniques presents et futurs, comme par exemple le Tevatron a 
Fermilab ou le LHC au CERN. 

Les sections efficaces de production des particules SUSY aupres des collisionneurs 
hadroniques ont ete etudiees en detail par le passe, aussi bien a I'ordre dominant 
(LO) [311 [321 |33] qu'a I'ordre sous-dominant (NLO) [Ml [351 [Ml [33 [Ml [Ml [IDl 
en QCD perturbative. II est connu que les corrections NLO QCD [M] et NLO 
SUSY-QCD completes [lO] pour la production d'une paire de sleptons augmentent 
les sections efficaces hadroniques d'environ 35% au Tevatron et 25% au LHC, ce 
qui etend le potentiel de decouverte des sleptons de plusieurs dizaines de GeV. 
Cependant, les corrections SUSY sont bien plus faibles que leurs analogues QCD en 
raison de la presence de squarks et gluinos tres lourds dans les boucles. 

Malgre le succes des premieres collisions proton-proton en mode polarise au 
collisionneur RHIC, les sections efficaces polarisees ont regu bien moins d'attention 
que leurs equivalents non polarises. Les calculs pionniers pour la production de 
squarks et de gluinos non massifs |4H B2] n'ont ete verifies, generalises au cas de 
particules SUSY massives et appliques aux collisionneurs actuels que recemment 
[13]. Concernant la production de sleptons, seulement des calculs negligeant 
les melanges entre les etats propres d'helicite etaient disponibles, et appliques 
uniquement a des experiences d'anciens collisionneurs [44j . 
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En raison de leurs couplages purement electrofaibles, les sleptons sont parmi les 
particules SUSY les plus legeres dans de nombreux scenarios de brisure de super- 
symetrie |45l 146) . Les sleptons et sneutrinos se desintegrent souvent directement 
en la particule SUSY la plus legere (le neutralino le plus leger dans les modeles de 
supergravite minimale (mSUGRA) ou le gravitino pour la brisure de supersymetrie 
vehiculee par interactions de jauge (GMSB)) et le partenaire du SM correspondant 
(un lepton ou un neutrino). Ainsi, un signal relatif a une paire de sleptons produite 
en collisionneur hadronique consistera en une paire de leptons tres energetiques, qui 
sera facilement detectable, et de I'energie manquante associee. 

Dans cette these, nous avons verifie les calculs pionniers pour la produc- 
tion polarisee d'une paire de sleptons |44j . que nous avons ensuite generalises 
afin de prendre en compte le melange des etats propres d'interaction, qui est 
surtout pertinent pour les sleptons de troisieme generation. Nous presentons les 
resultats analytiques pour les courants de sleptons neutres et charges, et predisons 
numeriquement les asymetries simple-spin pour le collisionneur RHIC et pour 
d'eventuelles ameliorations du Tevatron et du LHC ou I'un des faisceaux est 
polarise. Nous avons mis en evidence la sensibilite de I'asymetrie simple-spin a 
Tangle de melange du slepton tau et la possibilite de I'utiliser comme moyen pour 
distinguer le signal SUSY du bruit de fond du SM correspondant a la production 
Drell-Yan d'une paire de leptons |47j . 

Le bruit de fond standard principal lie a la production hadronique d'une paire 
de sleptons vient des desintegrations de paires WW et tt en une paire de leptons et 
de I'energie manquante \48\ I49j. Deux elements cles pour distinguer le signal SUSY 
du bruit de fond standard sont la reconstruction de la masse et la determination 
du spin des particules produites. Pour une paire de sleptons, la masse (s)transverse 
de Cambridge est une observable particulierement utile, puisqu'une connaissance 
precise du spectre en impulsion transverse (qt) suffit alors pour determiner la masse 
[50] et le spin [51] des sleptons. 

Lorsque que Ton etudie la distribution en impulsion transverse d'un systeme 
non colore produit avec une masse invariante M lors d'une collision hadronique, il 
est pertinent de separer les regions cinematiques relatives aux larges et aux faibles 
valeurs de qt- Dans la region des qx importants {qx > M), I'utilisation de la 
theorie perturbative a ordre fixe est parfaitement justifiee, vu que le developpement 
en serie de la distribution en qx est controle par un parametre d'expansion de 
faible valeur, la constante de couplage forte as(M^). Dans la region des petites 
valeurs de qt, les coefficients du developpement perturbatif sont amplifies par 
des termes logarithmiques importants, ln{]VP/q'^), et les resultats bases sur des 
calculs perturbatifs divergent pour qx 0, la convergence de la serie etant alors 
completement detruite. Ces logarithmes proviennent de remission multiple de 
gluons mous par I'etat initial, et doivent etre systematiquement resommes a tous 
les ordres en afin d'obtenir des resultats consistants. La methode pour effectuer 
cette resommation est bien connue [52l[5a[Hl[55llMl[5ZllMl[59l|6ni[6ll[62]. La 
resommation des logarithmes dominants a ete effectuee pour la premiere fois en 
[52] . II a ete montre en [53] que la procedure de resommation est plus naturellement 
effectuee dans I'espace du parametre d'impact b, b etant la variable conjuguee a 
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qt via une transformation de Fourier. En effet, dans ce cas-la, la cinematique de 
remission multiple de gluons factorise completement. Dans les cas particuliers de 
la production Drell-Yan d'une paire de leptons et de la production d'un boson 
electrofaible, la resommation dans I'espace b a ete effectuee au niveau sous-dominant 
(NLL) un formalisme de resommation consistant a n'importe quelle precision 
logarithmique a ete developpe [59j, et les termes d'ordre sous-sous-dominant ont ete 
calcules [BD]. Pour les valeurs de qt intermediaires, le resultat resomme doit etre 
ajuste de fagon consistante avec celui base sur la theorie perturbative, afin d'obtenir 
des predictions d'une precision theorique uniforme sur tout le domaine d'impulsion 
transverse considere. 

Dans ce travail, nous avons implemente le formalisme de resommation en qt 
propose en [611 [62] , et predit le spectre en impulsion transverse pour la production 
d'une paire de sleptons au LHC. Nous avons combine le resultat resomme (valide 
pour les faibles valeurs de qt), calcule au niveau NLL, avec la section efficace 
a ordre fixe (valide pour les larges valeurs de qr), calculee a VO{as) en QCD 
perturbative qui correspond a la production d'une paire de sleptons associee a un 
jet QCD [63]. II s'agit du premier calcul de precision concernant la distribution 
en impulsion transverse pour un processus de production d'une paire de particules 
SUSY aupres d'un collisionneur hadronique. Dans nos resultats numeriques, nous 
avons montre I'importance de la resommation aussi bien pour les faibles valeurs 
de qt que pour les valeurs intermediaires. Par ailleurs, la resommation permet de 
reduire la dependance de la distribution en qt en les echelles non physiques de 
factorisation et de renormalisation. Nous avons egalement etudie I'influence des 
contributions non perturbatives sur le resultat resomme, et observe qu'elle etait 
reduite par rapport a I'effet de la resommation. 

En ce qui concerne les corrections NLO SUSY-QCD, elles ont ete calculees 
uniquement en negligeant le melange des etats propres d'interaction des squarks 
apparaissant dans les boucles [40j . Nous avons generalise ce travail en incluant ce 
melange pertinent pour les squarks de troisieme generation, et avons considere les 
effets de seuil provenant de I'emission de gluons mous par I'etat initial. Lorsque les 
partons initiaux ont tout juste assez d'energie pour produire la paire de sleptons 
dans I'etat final, les corrections virtuelles et I'emission de gluons reels supprimee 
par I'espace de phase menent a I'apparition de termes logarithmiques importants 
a"[ln^'^~^(l — z)/{l — z)]^ a I'ordre n de la theorie perturbative, oii z = M'^/s, 
s est I'energie dans le centre de masse partonique et M la masse invariante de 
la paire de sleptons. Lorsque s est proche de M^, ces logarithmes doivent etre 
resommes a tous les ordres en as- Bien que ces divergences apparaissent de fagon 
explicite dans la section efficace partonique, la section efficace hadronique n'est 
en general pas divergente en raison de la convolution avec les densites de partons 
tres faibles pour les grandes valeurs de la fraction d'impulsion longitudinale du 
proton correspondant aux valeurs de z proches de un. La resommation en seuil est 
done plutot une tentative de quantification de I'effet d'un ensemble de corrections 
bien definies qu'une simple somme de logarithmes d'origine cinematique. Ces effets 
peuvent cependant etre significatifs meme loin du seuil hadronique, et Ton s'attend 
done a des corrections importantes pour la section efficace de production de type 
Drell-Yan d'une paire de sleptons de quelques centaines de GeV au Tevatron et au 
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LHC. 



La resommation en seuil a tous les ordres en equivaut a I'exponentiation des 
radiations de gluons mous, et n'a pas lieu dans I'espace z directement, mais dans 
I'espace A'^ de Mellin, oii iV est la variable conjuguee a z par une transformation 
de Mellin, la region du seuil z ^ \ correspondant a la limite N ^ oo. Ainsi, la 
section efficace resommee dans I'espace z sera obtenue apres une transformation 
inverse finale. La resommation en seuil pour le processus Drell-Yan fut d'abord 
effectuee en |64[ I65j aux niveaux logarithmiques dominant et sous-dominant (NLL), 
correspondant a la resommation des termes de type a" In^" et a" In'^""^ A^. 
L'extension au niveau NNLL (termes de type In^""^ N) a egalement ete effectuee 
a la fois pour le processus Drell-Yan [66] et pour la production d'un boson de 
Higgs [67j. II a ete montre |68l |69] que les contributions dues a remission de 
partons colineaires peuvent egalement etre incluses de fagon consistante dans la 
formule de resommation. Cela correspond au formalisme de resommation "ameliore 
colineairement" , oii des termes contenant un facteur suppressif et une classe 
de contributions universelles independantes de sont egalement resommes. Tres 
recemment, les contributions a I'ordre NNNLL (les termes de type a"ln^"~^ A^) ont 
ete calculees [TDJEIII^. 

Nous presentons ici une etude detaillee des effets de la resommation en seuil 
pour la production de type Drell-Yan d'une paire de sleptons et pour la production 
associee d'un slepton et d'un sneutrino dans le cadre de scenarios mSUGRA et 
GMSB. Nous avons ajuste les resultats resommes a la precision NLL, calcules 
grace au formalisme de resommation ameliore colineairement, avec les resultats 
bases sur la theorie perturbative calcules a la precision NLO. Numeriquement, nous 
avons montre une augmentation non negligeable de la section efficace theorique par 
rapport aux predictions NLO, et une stabilisation de la dependance en les echelles 
non physiques grace a I'apport des termes d'ordres superieurs pris en compte par la 
resommation ^73j. 

L'origine dynamique des contributions logarithmiques intervenant dans les 
formalismes de resommation en impulsion transverse et en seuil est identique, vu 
qu'il s'agit de remission multiple de gluons mous par I'etat initial. Un formal- 
isme de resommation jointe, prenant en compte simultanement les contributions 
des gluons mous dans les deux regions cinematiques concernees {qt <C M et 
~ s) a ete developpe dans la derniere decade [TH [75] • L'exponentiation des 
termes singuliers dans les espaces de Mellin et du parametre d'impact, pour la 
resommation en seuil et en impulsion transverse respectivement, a ete prouvee, 
et une methode consistante pour effectuer les transformations inverses a ete 
introduite afin d'eviter le pole de Landau et les singularites dues aux densites de 
partons. Les applications de ce formalisme a la production hadronique d'un photon 
rapide [76j, d'un boson electrofaible [77], d'un boson de Higgs [78j, et d'une paire de 
quarks lourds [79] montrent les effets de la resommation sur differentes distributions. 

Nous presentons dans ce travail un traitement joint des corrections a faible 
impulsion transverse et des contributions importantes proche du seuil partonique 
pour la production d'une paire de sleptons aupres des collisionneurs hadroniques, ce 
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qui permet une comprehension complete des effets de gluons mous pour le spectre 
en impulsion transverse et pour les distributions en masse invariante. Avec le travail 
sur la resommation en impulsion transverse [63] et la resommation en seuil |73j . 
cette etude [80] complete notre programme ayant pour but de fournir les premiers 
calculs de precision incluant la resommation de gluons mous pour la production de 
sleptons aupres des collisionneurs hadroniques. 

Si les particules SUSY existent, elles doivent aussi apparaitre dans les boucles 
de particules virtuelles et affecter les observables de precision electrofaibles et 
les observables a basse energie. Plus particulierement, les courants neutres a 
changement de saveur qui apparaissent seulement au niveau des boucles dans le 
SM contraignent severement les contributions de nouvelle physique au meme ordre 
perturbatif. Le MSSM se libere de ces contraintes grace aux hypotheses de Violation 
de Saveur Minimale contrainte (cMFV) |8H l82] ou de Violation de Saveur Minimale 
(MFV) [831 mi US]) ou les particules SUSY peuvent intervenir dans les boucles, 
mais les changements de saveur sont soit negliges, soit completement dictes par la 
structure des couplages de Yukawa et par la matrice CKM [861 EZ] ■ 

En SUSY avec MFV, les elements des matrices de masse des squarks violant la 
saveur decoulent des couplages trilineaires de Yukawa entre les supermultiplets de 
Higgs et de fermions et des differentes renormalisations des secteurs des quarks et des 
squarks via les equations du groupe de renormalisation qui induisent des violations 
de saveur supplementaires a I'echelle electrofaible [88l [891 EE IH]. En SUSY avec 
violation de saveur non minimale, des sources de violation de saveur additionnelles 
sont incluses dans les matrices de masse et leurs termes non diagonaux qui ne 
peuvent plus etre simplement deduits a partir de la matrice CKM seule doivent 
etre consideres alors comme des parametres libres. Dans ce travail, nous allons 
considerer le melange des saveurs de squark de deuxieme et troisieme generations, 
car d'une part, les recherches directes de violation de la saveur dependent des 
capacites a determiner la saveur, ce qui n'est experimentalement bien etabli que 
pour les saveurs lourdes, et d'autre part, des contraintes experimentales severes pour 
la premiere generation existent en raison de mesures tres precises des oscillations 

— et des premieres preuves du melange — [921 [931 (Mj • 

Nous avons analyse I'espace des parametres NMFV SUSY, recherchant les 
regions permises par les contraintes venant des mesures de precision electrofaibles, 
des observables a basse energie et des donnees cosmologiques. Nous avons observe 
que le melange des chiralites et des saveurs de deuxieme et troisieme generations est 
fortement contraint, notamment par I'erreur experimentale de plus en plus petite 
sur le rapport d'embranchement 6 — > S7 et la densite relique de matiere sombre. 
Nous avons defini quatre nouveaux points typiques avec leur ligne associee, valides 
a la fois en SUSY avec cMFV, MFV et NMFV, et pour lesquels nous presentons la 
dependance des masses de squarks et de la decomposition des etats physiques de 
squark en la violation de la saveur. 

Considerant la SUSY avec cMFV (le MSSM habituel), les corrections SUSY- 
QCD pour la production de squarks et de gluinos [Mj, de jauginos [IQ], ainsi 
que pour leur production associee [37j ont deja ete calculees. En raison de leur 
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couplage fort, les squarks devraient etre produits abondamment aux collisionneurs 
hadroniques, et I'espace de phase favorise la production des etats propres de masse 
les plus legers. Ainsi, les productions des squarks top [35] et bottom [95] avec 
un grand melange d'helicite ont vequ une attention toute particuliere. Dans cette 
these, nous nous sommes interesses a I'importance des canaux electrofaibles pour la 
production de paires de squarks non diagonales et mixtes de troisieme generation 
aux collisionneurs hadroniques [SO]. NaiVement, Ton s'attend a ce que ces sections 
efficaces, qui sont d'ordre deux en la constante de structure fine, 0{a'^), soient plus 
faibles que celles concernant la production forte d'une paire de squarks diagonale 
d'environ deux ordres de grandeurs. Pour la production non diagonale, I'importance 
des canaux QCD est reduite en raison de la presence de boucles, et celle des canaux 
electrofaibles I'est egalement en raison du couplage faible. L'importance relative de 
ces canaux merite done une etude approfondie. Si I'on considere des squarks bottom 
qui se melangent, leur contribution au niveau des boucles QCD doit egalement etre 
prise en compte. 

Ensuite, pour la premiere fois, nous nous sommes concentres sur les effets 
possibles de la violation de saveur non minimale (NMFV) aux collisionneurs 
hadroniques [97]. A cette fin, nous avons recalcule toutes les amplitudes d'helicite 
pour la production et la desintegration des squarks et des jauginos, en prenant 
en compte les interactions non diagonales des courants charges des jauginos et les 
interactions de Yukawa des Higgsinos, et en generalisant les matrices de melange 
d'helicites bidimensionnelles, souvent supposees reelles, en matrices de melange 
d'helicites et de saveurs, complexes et six-dimensionnelles. Nous avons verifie 
que nos resultats reproduisaient ceux de la litterature existant dans les limites de 
squarks non melanges. 

Dans notre analyse phenomenologique de la production NMFV de squarks et de 
jauginos, nous nous sommes concentres sur le LHC en raison de son energie dans 
le centre de masse elevee et de sa luminosite importante. Nous avons porte une 
attention particuliere a la competition entre les effets lies aux densites de partons 
qui sont domines par les contributions des quarks legers, les contributions fortes du 
gluino qui sont generalement plus importantes que les contributions electrofaibles et 
qui ne doivent pas necessairement etre diagonales en saveur, et la presence de saveurs 
lourdes dans I'etat final, facilement identifiables experiment alement et generalement 
plus legeres que les saveurs de squark de premiere et deuxieme generations. 
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Chapter 1 

Introduction 



The Standard Model (SM) of particle physics [D EJ El H El H [7] provides a 
successful description of all experimental high energy data. However, despite of 
its success many fundamental questions remain unanswered, e.g. the origins of 
electroweak symmetry breaking and particle masses, the large hierarchy between 
the electroweak and the Planck scales, the mechanism leading to neutrino oscilla- 
tions, the origins of dark matter and of the cosmological constant, or the strong 
CP-problem. Attempts to relate different SM parameters lead to more fundamental 
theories, that may at the same time solve some of the open problems of the SM. 

The basic philosophy of Grand Unified Theories (GUTs) [U El |10] is to consider 
the SM symmetry groups as originating from the breaking of a larger simple group. 
At the GUT scale, the three SM gauge coupling constants unify and quarks and 
leptons are embedded in common representations of the unifying gauge group. 
These theories include then number of additional gauge bosons, leading potentially 
to interactions violating the baryon and lepton numbers. This leads to one of 
the major phenomenological problems of GUTs, which predict baryon instability, 
contrary to the experimental non-observation of proton decay. GUT theories can 
explain the quantization of the electric charge and incorporate massive neutrinos, 
but have difficulties in accounting for the measured value of the electroweak mixing 
angle. Besides, many other conceptual SM problems remain unsolved. 

One popular approach to solve the hierarchy problem of the SM is to extend 
space-time to higher dimensions |111 [T2j. In this framework, the gravitational and 
gauge interactions become unified close to the weak scale, which is then the only 
fundamental scale of the theory. The large value of the Planck scale is only a 
consequence of the new dimensions. The usual four-dimensional space is contained 
in a four-dimensional "brane", embedded in a larger structure with N additional 
dimensions, the "bulk". In these theories, each field of the SM possesses a tower 
of Kaluza-Klein excitations with the same quantum numbers, but different mass, 
which have not been observed at the present time, but which should lie in the 
TeV-range. They could then be detected at the future Large Hadron Collider (LHC) 
at CERN. 

Recently, other attempts to solve the hierarchy problem have been proposed, 
e.g. Little-Higgs or Twin-Higgs theories, which predict new particles with masses 
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in the TeV-range as well [131 [T5] . These theories include fermionic partners for 
quarks and leptons and bosonic partners for the SM gauge bosons, which allows for 
stabilization of the Higgs mass beyond tree-level thanks to a non-linearly realized 
symmetry that relates the couplings to the Higgs in such a way that the quantum 
corrections to the Higgs mass cancel. 

In this thesis, we focus on another attractive extension of the SM, supersymme- 
try (SUSY) [161 El UHl [IHl EO], and more precisely the Minimal Supersymmetric 
Standard Model (MSSM) [211 122j . Weak scale supersymmetry provides a natural 
solution for a set of conceptual problems of the SM. Linking fermions with bosons, 
SUSY allows for a stabilization of the gap between the Planck scale and the 
electroweak scale [231 El] and for a consistent unification of SM gauge couplings 
at high energies [251 l26l [271 I28j . In addition, it can include a potential dark 
matter candidate as the stable lightest SUSY particle [291 EOj • Since spin partners 
of the SM particles have not yet been observed and in order to remain a viable 
solution to the hierarchy problem, SUSY must be broken at low energy via soft 
mass terms in the Lagrangian. As a consequence, the SUSY particles must 
be massive in comparison to their SM counterparts, but their mass should not 
exceed a few TeV. A conclusive search covering a wide range of masses up to the 
TeV scale is then one of the main topics in the experimental program at present 
and future hadron colliders, such as the Tevatron at Fermilab and the LHC at CERN. 

Production cross sections for SUSY particles at hadron colliders have been 
extensively studied in the past at leading order (LO) [311 [32l [33] and also at 
next-to-leading order (NLO) of perturbative QCD [SI [SS [Ml [SZl [Ml [Ml [10] • The 
NLO QCD ^39j and full SUSY-QCD ^40j corrections for slepton pair production are 
known to increase the hadronic cross sections by about 35 % at the Tevatron and 
25% at the LHC, extending thus the discovery reaches of sleptons by several tens of 
GeV. However, the presence of massive squarks and gluinos in the loops makes the 
genuine SUSY corrections considerably smaller than the standard QCD ones. 

Despite of the first successful runs of the RHIC collider in the polarized pp mode, 
polarized cross sections have received much less attention. Only the pioneering 
LO calculations for massless squark and gluino production [411 US] have recently 
been confirmed, extended to the massive case, and applied to current hadron 
colliders [l3]. Concerning slepton pair production, only polarized calculations for 
non mixing sleptons were available before, and for old collider experiments only [44j . 

Due to their purely electroweak couplings, sleptons are among the lightest 
SUSY particles in many SUSY-breaking scenarios [451 SS] • Sleptons and sneutrinos 
often decay directly into the stable lightest SUSY particle (the lightest neutralino 
in minimal supergravity (mSUGRA) models or the gravitino in gauge-mediated 
SUSY-breaking models (GMSB)) plus the corresponding SM partner (lepton or 
neutrino). As a result, a slepton signal at hadron colliders will consist in a highly 
energetic lepton pair, which will be easily detectable, and associated missing energy. 

In this thesis, we verify the pioneering polarized calculations for slepton pair 
production [44j and extend them by including the mixing of the left- and right- 
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handed interaction eigenstates relevant for third-generation sleptons. We present 
analytical results for neutral and charged current sleptons and make numerical 
predictions for longitudinal spin asymmetries at RHIC and possible upgrades of 
the Tevatron and the LHC, where one of the beams is considered to be polarized. 
We put particular emphasis on the sensitivity of the asymmetry to the tau slepton 
mixing angle as predicted by various SUSY-breaking mechanisms. Possibilities of 
using asymmetries to discriminate between the SUSY signal and the corresponding 
SM Drell-Yan background are also discussed |47j . 

The main SM background to slepton pair production at hadron colliders is due 
to WW and tt decays to a lepton pair and missing energy [ISl [l9] . Two key features 
distinguishing the SUSY signal from the SM background are the reconstruction of 
the mass and the determination of the spin of the produced particles. For sleptons, 
the Cambridge (s)transverse mass proves to be a particularly useful observable, 
requiring only a precise knowledge of the transverse-momentum (gy) spectrum to 
get their mass [SO] and spin [5T] . 

When studying the qr-distribution of a colourless system produced with an 
invariant-mass M in a hadronic collision, it is appropriate to separate the large-gy 
and small-gr regions. In the large-gy region (qt > M) the use of fixed-order 
perturbation theory is fully justified, since the perturbative series is controlled by a 
small expansion parameter, the strong coupling constant as(M^). In the small-g-r 
region, where the coefficients of the perturbative expansion in asiM"^) are enhanced 
by powers of large logarithmic terms, ln(M^/g|,), results based on fixed-order calcu- 
lations diverge as qx — > 0, and the convergence of the perturbative series is spoiled. 
These logarithms are due to multiple soft-gluon emission from the initial state and 
have to be systematically resummed to all orders in ag in order to obtain reliable 
perturbative predictions. The method to perform all-order soft-gluon resummation 
at small qr is well known [SaESIMlEaESlEZlESlEaEniEIlES]. The resummation 
of leading logarithms was first performed in [52]. It was shown in [53j that the 
resummation procedure is most naturally performed using the impact-parameter [h) 
formalism, where b is the variable conjugate to qx through a Fourier transformation, 
to allow the kinematics of multiple-gluon emission to factorize. In the special case 
of Drell-Yan lepton pair or electroweak boson production, 6-space resummation was 
performed at next-to-leading level in [55] , a resummation formalism consistent at any 
logarithmic accuracy was developed in [59], and the next-to-next-to-leading order 
terms have been calculated in [60]. At intermediate qx the resummed result has to be 
consistently matched with fixed-order perturbation theory in order to obtain predic- 
tions with uniform theoretical accuracy over the entire range of transverse momenta. 

We implement the universal g^-resummation formalism proposed in |6H [62] and 
compute the gT-distribution of a slepton pair produced at the LHC by combining re- 
summation at small qx and the fixed-order cross section at large qx- The resummed 
contribution has been computed at the next-to-leading logarithmic (NLL) accuracy 
and the fixed-order cross section at 0{as) in perturbative QCD, corresponding to 
the production of a slepton pair plus a QCD jet. It is the first precision calculation 
of the gr-spectrum for SUSY particle pair production at hadron colliders. The 
importance of resummed contributions at small and intermediate values of gy, both 
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enhancing the pure fixed-order result and reducing the scale uncertainty, is shown 
in our numerical results i63j . 



Concerning NLO SUSY-QCD corrections [IQ], they have only been computed 
for non-mixing squarks appearing in the loops. We extend this work by including 
the mixing effects relevant for the third generation in the squark sector, and 
we consider the threshold-enhanced contributions of the QCD corrections |39j . 
also due to soft-gluon emission from the initial state. They arise when the 
initial partons have just enough energy to produce the slepton pair in the final 
state. In this case, the mismatch between virtual corrections and phase-space 
suppressed real-gluon emission leads to the appearance of large logarithmic terms 
a"[ln^"~^(l — z)/(l — z)]^ at the n*^ order of perturbation theory, where z = M'^/s, 
s being the partonic centre-of-mass energy and M the slepton pair invariant-mass. 
When s is close to M^, these large logarithms have to be resummed to all orders in 
as- Although they are manifest in the partonic cross section, they do not generally 
result in divergences in the physical cross section since they are smoothed by the 
convolution with the steeply falling parton distributions. Threshold resummation is 
then not really a summation of kinematic logarithms in the physical cross section, 
but rather an attempt to quantify the effect of a well-defined set of corrections 
to all orders, which can be significant even if the hadronic threshold is far from 
being reached. Large corrections are thus expected for the Drell-Yan like produc- 
tion of a slepton pair with invariant-mass of a few 100 GeV at the Tevatron and LHC. 

All-order resummation is achieved through the exponentiation of the soft-gluon 
radiation, which does not take place in z-space directly, but in Mellin A^-space, 
where A^ is the Mellin variable conjugate to z and the threshold region z — > 1 
corresponds to the limit A^ oo. Thus, a final inverse Mellin transform is needed 
in order to obtain a resummed cross section in z-space. Threshold resummation 
for the Drell-Yan process was first performed in |64t [55] at the leading logarithmic 
and next-to-leading logarithmic (NLL) levels, corresponding to terms of the form 
a'^ln^^'N and a^ln^'^-^AT. The extension to the NNLL level (a^ln^'^-^Ar terms) 
has been carried out both for the Drell-Yan process [66] and for Higgs-boson 
production [67J . It was shown in |68[ [69] that contributions due to collinear parton 
emission can be consistently included in the resummation formula, leading to a 
"collinear- improved" resummation formalism where 1/A^-suppressed and a class of 
A^-independent universal contributions are resummed as well. Very recently, even 
the NNNLL contributions (a^ln^""^ A^ terms) became available [70 | [7T ] [72] . 

We present here an extensive study on NLL threshold resummation effects for 
Drell-Yan like slepton pair and slepton-sneutrino associated production in mSUGRA 
and GMSB scenarios, matching the resummed contributions computed within a 
collinear-improved resummation formalism with a fixed-order calculation at NLO 
accuracy. Numerically, we show a non-negligible increase of the theoretical cross 
sections with respect to the NLO prediction and a stabilization of the unphysical 
scale dependences thanks to the higher order terms taken into account in the 
resummed component of the cross section [73]. 

The dynamical origin of the enhanced contributions is the same both in 
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transverse-momentum and threshold resummations, since it comes from the soft- 
gluon emission by the initial state. A joint resummation formalism, embodying 
soft-gluon contributions in both the delicate kinematical regions (gr <C M, ~ s) 
simultaneously, has been developed in the last decade [TUITS]. The exponentiation 
of the singular terms in the Mellin and impact-parameter spaces, for threshold and 
transverse-momentum resummation respectively, has been proven, and a consistent 
method to perform the inverse transforms in order to avoid the Landau pole and the 
singularities of the parton distribution functions has been introduced. Applications 
to prompt-photon [76j, electroweak boson [77j, Higgs boson [75] and heavy-quark 
pair ^79j production at hadron colliders show the substantial effects of the joint 
resummation on the differential cross sections. 

We present a joint treatment of the recoil corrections at small qx and the 
threshold-enhanced contributions near partonic threshold for slepton pair produc- 
tion at hadron colliders, allowing for a complete understanding of the soft-gluon 
effects in differential distributions [SD]. Together with the previous papers on 
transverse-momentum [63] and threshold [73] resummation, this completes our pro- 
gram of providing the first precision calculations including soft-gluon resummation 
for slepton pair production at hadron colliders. 

If SUSY particles exist, they should also appear in virtual particle loops and 
affect low-energy and electroweak precision observables. In particular, flavour- 
changing neutral currents (FCNC), which appear only at the one- loop level even 
in the SM, put severe constraints on new physics contributions appearing at the 
same perturbative order. The MSSM has passed these crucial tests, largely due 
to the assumption of constrained Minimal Flavour Violation (cMFV) }8H [82] 
or Minimal Flavour Violation (MFV) [S31 [Ml EE], where heavy SUSY particles 
may appear in the loops, but flavour changes are either neglected or completely 
dictated by the structure of the Yukawa couplings and thus the CKM-matrix [521 [57] . 

In SUSY with MFV, the flavour violating entries in the squark mass matrices 
stem from the trilinear Yukawa couplings of the fermion and Higgs supermultiplets 
and the resulting different renormalizations of the quark and squark mass matrices, 
which induce additional flavour violation at the weak scale through renormalization 
group running [88[ [89l [90l [9T] . In non- minimal flavour violating SUSY, additional 
sources of flavour violation are included in the mass matrices at the weak scale, 
and their flavour-violating off-diagonal terms cannot be simply deduced from the 
CKM matrix alone, and have to be considered as free parameters. In this work, 
we consider flavour mixings of second- and third-generation squarks, since direct 
searches of flavour violation depend on the possibility of flavour tagging, which is 
established experimentally only for heavy flavours. In addition, stringent experi- 
mental constraints for the first-generation are imposed by precise measurements of 
mixing and first evidence of — mixing |92[ [931 IH] ■ 

Considering SUSY with cMFV (the usual MSSM), NLO SUSY-QCD calculations 
for the production of squarks and gluinos ^34j; gauginos [40j . as well as for their 
associated production [37] are available. Due to their strong coupling, squarks 
should be abundantly produced at hadron colliders. In addition, phase space 
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favours the production of the hghter of the squark mass eigenstates of identical 
flavour. As a consequence, the production of top [35] and bottom [95j squarks with 
large helicity mixing has received particular attention. In this thesis, we investigate 
the importance of electroweak channels for non-diagonal and mixed squark pair 
production at hadron colliders [HS]- Naively, one expects these cross sections, which 
are of 0(0^) in the fine structure constant a, to be smaller than the diagonal strong 
channels by about two orders of magnitude. For non-diagonal squark production, 
the interplay between loop suppression in QCD and coupling suppression in the 
electroweak case merits a detailed investigation, and in the presence of the mixing 
of bottom squarks, their loop contributions must also be taken into account. 

Then, for the first time, we concentrate on the possible effects of non-minimal 
flavour violation (NMFV) at hadron colliders |97j. To this end, we recalculate 
all squark and gaugino production and decay helicity amplitudes, keeping at the 
same time the CKM-matrix and the quark masses to account for non-diagonal 
charged-current gaugino and Higgsino Yukawa interactions, and generalizing the 
two-dimensional helicity mixing matrices, often assumed to be real, to generally 
complex six-dimensional helicity and generational mixing matrices. 

In our phenomenological analysis of NMFV squark and gaugino production, we 
concentrate on the LHC due to its larger centre-of-mass energy and luminosity. 
We pay particular attention to the interesting interplay of parton density functions 
(PDFs), which are dominated by light quarks, strong gluino contributions, which are 
generally larger than electroweak contributions and need not be flavour-diagonal, and 
the appearance of third-generation squarks in the final state, which are easily iden- 
tified experimentally and generally lighter than first- and second-generation squarks. 

This thesis is organized as follows. In the first part of Chapt. [21 we briefly de- 
scribe the MSSM within cMFV, showing several motivating arguments for SUSY and 
defining the model. In the second part of this chapter, we set up the notations that 
we use in the case of NMFV SUSY, and we perform a precise numerical analysis of 
the experimentally allowed NMFV SUSY parameter space with respect to low-energy 
constraints, leading to the definition of four collider-friendly benchmark points for 
which we investigate the corresponding helicity and flavour decomposition of the 
up- and down- type squarks. Finally, we introduce generalized couplings in order to 
compute compact analytical expression for the various cross sections calculated in 
this work. In Chapt. [3l we describe the two gy-resummation formalisms (CSS and 
universal), the threshold-resummation formalism and the joint-resummation formal- 
ism that we have used in the case of slepton pair hadroproduction. Chapts. [4] and 
[S] are devoted to the results, the first one for slepton pair hadroproduction and the 
second one to squark and gaugino production and decays, and we show a large num- 
ber of analytical and numerical results. Our conclusion and outlook are presented 
in Chapt. [H 
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Chapter 2 

The Minimal Supersymmetric 
Standard Model 



2.1 Motivation 



The Standard Model of particle physics is a gauge field theory based on the 
symmetry group SU{3) x SU{2) x U{1), containing the electroweak [U El El Uj and 
the strong O [6l [7| interactions and providing a remarkably accurate description of 
a large class of phenomena. It is well-established by the discovery of all its particle 
content, the Higgs boson excepted, and by precision measurements at colliders. 
However, a new framework will certainly be required, at least at the Planck scale 
Mp = 2.8 X 10^^ GeV, where the quantum gravitational effects become important, 
but more probably at a lower scale, since the absence of new physics between the 
current experimental limit of several hundreds of GeV and Mp is highly improbable. 
Moreover, this large gap leads to what is called the hierarchy problem 



The origin of the masses in the SM is an isodoublet scalar Higgs field 
11001 11011 11021 1103j , yielding electroweak symmetry breaking and one physical Higgs 
boson which couples to each SM fermion of mass nif with a strength driven by the 
Yukawa couplings A/. The quantum corrections Am|^ to the Higgs squared mass 



m^, where 



ml = {mjj)o + Ami, (2-1) 
("T'h)o being a parameter of the fundamental theory, are given by 



2 



Idtt^ 



-2A2 +6mf In — + 



(2.2) 



where A is an ultraviolet cutoff that corresponds to the scale at which new physics 
alters the theory and that regulates fermionic loop-integrals. If A is of order Mp, 
these corrections are some 30 orders of magnitude larger than the expected value of 
about (100 GeV)^, particularly due to the large top Yukawa coupling. Supersym- 
metry [16l [171 UHl IISI l20l EH E2] provides an elegant solution to this problem, since 
a heavy complex scalar particle of mass ms can couple to the Higgs boson with a 
strength A5. The corresponding quantum corrections to ml are 



Ami 



16 7r2 



A^ - 2 ml ln— + 

* ms 



(2.3) 
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Provided that the scalar masses are not too heavy, the systematic cancellation of 
fermionic and bosonic corrections to the Higgs squared mass is then achieved, since 
each fermion of the SM is now accompanied by two scalars, and their couplings to 
the Higgs field are closely related by A5 = |Ajp. The remaining corrections to the 
squared Higgs mass depend only logarithmically on the cutoff A and are thus under 
control. 

Furthermore, a set of conceptual problems of the SM can be solved thanks to 
supersymmetry, such as the unification of the fundamental gauge interactions |25ll26t 
\27\ [28] , since the SUSY particles modify the renormalization-group evolution of the 
gauge couplings with the energy, leading to unification at about 10^^ GeV. SUSY 
also provides a potential cold dark matter candidate, the lightest SUSY particle 
(LSP) |29[ I30j. and can even include gravity, in the framework of local supergravity 
theories [TMfTn5] . 

2.2 Definition of the model 

Most present SUSY models are based on the four-dimensional supersymmetric 
field theory of Wess and Zumino [16], which are free of many of the divergences 
encountered in similar SUSY theories of that time |106i 1107] . The simplest model is 
the straightforward supersymmetrization of the SM with the same gauge interactions, 
called the Minimal Supersymmetric Standard Model |2H 122] . As for any SUSY 
model, it postulates a symmetry between fermionic and bosonic degrees of freedom 
in nature, predicting thus the existence of a fermionic (bosonic) SUSY partner for 
each bosonic (fermionic) SM particle. A complete introduction to supersymmetric 
field theories and the MSSM can be found in Refs. pIM fTIM fTTn] . 

2.2.1 Field content 

Due to the various conserved quantum numbers of the known bosons and 
fermions, a minimal supersymmetric model cannot be built up with the SM 
particles alone, and new particles have to be postulated. Quarks and leptons 
get scalar partners called squarks and sleptons, while electroweak bosons and 
gluons get fermionic partners referred to as gauginos and gluinos. Since in the 
SM, the left- and right-handed parts of the fermionic fields transform differently 
under the gauge group, it is required that fermions get two superpartners, named 
left- and right-handed sfermions. Finally, to preserve the electroweak symmetry 
from gauge anomaly and to give masses to both up- and down-type fermions, 
the MSSM requires two Higgs doublets and their fermionic superpartners, the 
Higgsinos. The field content of the MSSM is shown in Tab. 12.11 Let us note that the 
local version of supersymmetry includes a spin-2 state and a spin-3/2 state which 
can be interpreted as the spin-2 graviton and its spin-3/2 superpartner, the gravitino. 

All of these particles are organised in chiral and gauge supermultiplets, containing 
an equal number of fermionic and bosonic degrees of freedom. Chiral supermultiplets 
contain one on-shell Weyl fermion ip (i.e. the left- or the right-handed part of a 
fermionic field) and its associated scalar complex field (p, corresponding to a total 
of two fermionic and two bosonic real degrees of freedom. Off-shell Weyl fermions 
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Table 2.1: Field content of the MSSM. 



Names 


particle 


spin 


superpartner 


spin 


(s) quarks 


Q 


{uLdi) 


1/2 







(x3 families) 


u 


a. 


1/2 


~ * 







d 




1/2 


d* 





(s)leptons 


L 




1/2 


(z^cl) 





(x3 families) 


e 




1/2 


e*R 





Higgs(inos) 


Hu 







{Hi Hi) 


1/2 




Hd 









1/2 


gluon/gluino 


9 


1 


9 


1/2 


W bosons/winos 




1 




1/2 


B boson / bino 


B 


1 


B 


1/2 



having two additional real fcrmionic degrees of freedom, an auxiliary complex scalar 
field F is introduced, preserving supersymmetry off-shell, but being eliminated when 
one goes on-shell by imposing its equations of motion. Gauge supermultiplets contain 
a massless gauge boson and its on-shell associated fermionic partner A", which 
corresponds as well to two fermionic and two bosonic real degrees of freedom. If one 
goes off-shell, the fermionic field gets two additional real degrees of freedom, while 
the vector boson only gets one. As for chiral supermultiplets, an auxiliary field D 
with one real bosonic degree of freedom is introduced, preserving SUSY off-shell and 
being eliminated on-shell through its equations of motion. 

2.2.2 Lagrangian density 

The MSSM gauge interactions are the same as those of the SM and are deter- 
mined by the gauge group SU{3) x SU{2) x U{1). The full Lagrangian density for 
a renormalizable SUSY theory is then given by 

C = >Cchirai + >Cgauge - V25(0*r»A« - V2gX^%i,^T''cl>) + g {ct>* T'^ cj>) D'' . (2.4) 

-Cchirai ^.nd >Cgauge contain the kinetic terms and the gauge interactions for the chiral 
and gauge supermultiplets, 

>Cgauge = -^F;,F'^'^"-iAtV'D^A» + ^L>»D«, (2.6) 
where a are the Pauli matrices 
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The last terms of Eq. (|2.4p are interactions whose strength is given by the usual 
gauge couplings, but which are not gauge interactions from the point of view of an 
ordinary gauge theory. The covariant derivatives 



make the Lagrangian gauge-invariant, being hermitian matrices corresponding 
to the representation in which the chiral supermultiplets transform under the gauge 
group and satisfying \T°',T^] = if°'^'^T'^, where /"^^^ are the totally antisymmetric 
structure constants defining the group. Finally, in Eq. ()2.6p . F^^^ is the usual Yang- 
Mills field strength 

= d^Al - d,Al + gr'^A^Al. (2.9) 

The interactions between the chiral and the gauge supermultiplets are embodied in 
derivatives of the superpotential W 

W^ = ^ and W'^ = (2.10) 

0(j)i d(t)id(t)j 

W being analytic in the complex fields (f). Let us note that the auxiliary fields Fi do 
not explicitly appear in the Lagrangian since they are expressed in terms of VFj, using 
the equations of motion Fi = —W^ and F^* = —W^. This Lagrangian is obviously 
invariant under a global infinitesimal SUSY transformation e, the transformation 
rules being 

5(j)i = eipi, Sipia = i{a^e^)a^^i4>i + ^aFi, 
= (e%X- + At^e) , '^A^ = ^(a^a-6)„F;, + ^e„Z)^ 

5Fi = ieta'^D^Vi + V2g{T''<P)ie^ X^^ , SD'' = (eV'D^A'^ - D^At^^e) . (2.11) 

The MSSM is specified by the choice of its superpotential 

W = -ey^LHd-uy^QHu-dyAQHd + iiHuHd, (2.12) 

where e, u, d, L, Q, and Hd are the chiral and Higgs superfields described in 
previous subsection. The 3x3 Yukawa matrices y give rise to the masses of the 
quarks and leptons when the Higgs fields acquire their vacuum expectation values 
(vevs). As in the SM, different rotations are needed to diagonalize both the up- and 
down- type quark Yukawa matrices yu and yd, leading to the usual flavour mixing 
driven by the CKM matrix [861 I87j . Finally, the ^u-term provides the Higgs and 
Higgsino squared-mass terms. This superpotential is minimal since it is sufficient to 
produce a phenomenologically viable model. However, other gauge-invariant terms 
could be included in the superpotential, 

Wfi = ^XLLe + X'LQd+ ^X"udd + n' L Hu, (2.13) 

violating either the total lepton number L, or the baryon number B. In principle, we 
could just postulate B and L conservation, forbidding then the terms of Eq. ()2.13p . 
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as there are no possible renormahzable terms in the SM Lagrangian violating B or L. 
But neither B nor L are fundamental symmetries of nature since they are violated 
by non-perturbative electroweak effects Therefore, an alternative symmetry is 

rather imposed, forbidding the B- and L- violating terms of Eq. (|2.13p . the i?-parity 
[TO] , It is defined by 

R=(^_lfB+L+2S^ (2.14) 

S being the spin of the particle. The SM particles then have a positive i?-parity, 
while their SUSY counterparts have a negative one. Due to i?-parity conservation, 
interaction vertices have to contain an even number of SUSY particles, leading to 
SUSY particle production by pairs at colliders and to decays into states containing 
an odd number of stable LSPs, which can only interact via annihilation vertices. 
Furthermore, if we assume the LSP to be electrically and colour neutral, it can 
even be a potential dark matter candidate [30]. In this work, we assume -R-parity 
conservation, but an introduction to SUSY models with i?-parity violation can be 
found in Refs. [mifTT^ . 

The scalar potential is already included in the Lagrangian, via the F- and D- 
terms, which can be expressed as a function of the scalar fields, 

r) = F'* F, + ^D-D- = w^W: + ^Yl ali^T'^'Pf, (2.15) 

a 

where we sum over the different gauge groups. We should note that V is entirely 
determined by the other interactions of the theory, contrary to the SM potential 
containing free parameters. 



2.2.3 Soft SUSY-breaking 

SUSY particles still remain to be discovered, and their masses must therefore 
be considerably larger than those of the corresponding SM particles, so that super- 
symmetry must be broken. In order to remain a viable solution to the hierarchy 
problem, SUSY can, however, only be broken via soft terms in the Lagrangian, i.e. 
with positive mass dimension, which prevents us from introducing new quadratic 
divergences in the quantum corrections to the Higgs squared mass. Since we do not 
know the SUSY-breaking mechanism and at which scale it occurs, we usually modify 
the Lagrangian at low energies by adding all possible terms breaking explicitly the 
SUSY. In the MSSM, we get [IH] 

Csoft = -^(^M3gg + M2WW + MiBB + c.c}j 

- iuqQ + mf^ L + uin^u^ — dm^d -|-em|e^^ 

- {ml,^ HI Hu + Hi H^ + {h H^ H^ + c.c.)) 

- {ust.^iQHu-'daiAQHd-eaieLHd + c.c}j. (2-16) 

The first line contains the gluino, wino and bino mass terms and the second line 
the squark and slepton mass terms, iHq, m^, m^, m^, m§ being 3x3 hermitian 
matrices in family space. In the third line of £softj we have the mass terms for the 
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Higgs fields contributing to the scalar potential and in the fourth line the trilinear 
scalar interactions, a^, ad, and ag, which are also 3x3 matrices in generation space. 
Contrary to the supersymmetric part of the MSSM Lagrangian, which has only one 
new free parameter ^, the SUSY-breaking Lagrangian contain 105 masses, phases 
and mixing angles, which cannot be rotated away by redefining the field basis and 
which have no counterpart in the SM [115j. Most of them are strongly constrained, 
since they introduce new sources of flavour mixing and CP violation, which could 
enhance processes severely restricted by experiment, such as — , — and 
mixing, flavour-changing neutral-current (FCNC) B, fi or t decays, and so 
forth pTB] . 

Present models assume that SUSY is broken in a hidden sector, containing 
particles that have no or small coupling to the visible sector, and SUSY-breaking is 
mediated to the visible sector via an interaction shared by the two sectors. Let us 
note that if the mediating interaction is flavour-blind and the related parameters 
are real, we get automatically conditions which evade the flavour and CP violating 
terms in £soft at low energies. In local SUSY theories, the Goldstone fermion related 
to SUSY breaking, the goldstino, is absorbed by the gravitino which acquires a 
mass through the super-Higgs mechanism jll7l 1118] . analogously to the usual Higgs 
mechanism where the electroweak gauge bosons acquire a mass by absorbing the 
Goldstone bosons associated to electroweak symmetry breaking. We have studied 
two SUSY-breaking scenarios, supergravity [1191 11201 fT2l J and gauge-mediated 
SUSY-breaking [Ml |123l UMl 1125] . 

In the framework of supergravity, SUSY-breaking is mediated to the 
MSSM through gravitational interactions, appearing in the Lagrangian via non- 
renormalizable terms suppressed by powers of the Planck mass. Assuming minimal 
supergravity, the soft terms in £soft are completely determined by five parameters, 
the universal scalar and gaugino masses mo and the universal trilinear coupling 

the ratio of the vevs of the two neutral Higgs fields, tan f3, and the sign of /x. At 
the SUSY-breaking scale, we get the relations 

Mg = M2 = Mi = mi/2, (2.17) 
ml^=ml = m^ = ml = ii4 = mll; m]j^ = m]j^= ml, (2.18) 
au = ^oYu, ad = ^oYd, ae = vloye- (2.19) 

Assuming that all of these parameters are real, the problematic flavour and CP 
violation terms of the Lagrangian are automatically suppressed. Low-energy 
parameters are deduced from renormalization-group evolution of the high scale 
parameters down to the electroweak scale. 

SUSY-breaking can also be mediated by usual gauge interactions. To this aim, 
new chiral supermultiplets are introduced in the Lagrangian, the messenger fields, 
which carry SU{2>) x SU{2) x U{1) quantum numbers and couple to the particles of 
both the visible and hidden sectors. Virtual loops of messengers generate gaugino 
and sfermion masses, i.e. the mass terms in £soft! in a completely renormalizable 
framework. At the messenger scale Mmesj the trilinear couplings are generated via 
two-loop diagrams, and we can then neglect them with respect to the masses gener- 
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ated by one-loop diagrams, 



ad 



0. 



(2.20) 



Gauge interactions being flavour-blind and assuming real SUSY-breaking param- 
eters, undesired large FCNC and CP violation effects are again avoided. SUSY- 
breaking is parameterized by one gauge-singlet chiral superfield, whose scalar and 
auxiliary components acquire the two vevs {S) and {Fs), respectively, defining the 
scale A = {Fs)/{S). At the scale Mmes, the soft SUSY-breaking mass parameters 
are given by 



4tt 



M„ 



2(iV5 + 3iVio)A2/ 



A 



Mr, 



ai(M„ 



1 2 



(2.21) 



q, (2.22) 



where the coefficients Cj are quadratic Casimir invariants and Ni is the multiplicity 
of the messengers in the 5-1-5 and 10 -|- 10 vector-like supermultiplets, assuming 
that messengers come in complete multiplets of SU{5) to preserve gauge coupling 
unification. The threshold functions are given by \126\ I127j 



9[x) 



^ [(1 + x) log(l + x) + {l-x) log(l - x)], 



X 

1 + x 

x"^ 



log(l + x) - 2Li2 



X \ 1^ / 2x 
l + x) ^ 2 ^HlTx 



(2.23) 

+ {x ^ -x). (2.24) 



The different terms of >Csoft are then determined via the evolution of six parameters, 
the numbers of (s)quark and (s)lepton messenger fields rig and rij- needed for the 
calculation of the multiplicities, the messenger scale Mmes) tan/3 and the sign of 
/i. 



2.2.4 Electroweak symmetry breaking and particle mixing 

The scalar potential of the MSSM is given by 



V 



+ mljH^l^ + (|^|2 + mjjJ\Ha\^ + b {H^H^ + c.c.) 



1 



+ -(g^ + g'^){\H^\ 



\H,\y + \g^\H:HT + HlH, 



(2.25) 



The terms proportional to ^, g and g' come from the F- and D-terms of Eq. ()2.15p . 
while the terms proportional to mn^, rnH^ and h come from the SUSY-breaking 
Lagrangian of Eq. ()2.16p . SU{2) gauge transformations allow to rotate away any 
possible vev of one of the charged Higgs fields, and we simply take {H^) = 0. 
Since the minimum of the potential satisfies dV/dH^ = 0, we automatically get 
{H'^) = 0, which leads to electric charge conservation in the Higgs sector. As for 
U{1) gauge transformations, they allow to redefine the phases of and Hd, so 
that all complex phases of the Lagrangian can be absorbed, making the vevs real 
and positive. CP is thus not spontaneously broken by the Higgs scalar potential, 
and the Higgs physical states are also CP eigenstates. 
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To get electroweak symmetry breaking, the potential has to have a minimum, 
which is the case if 



M^+mj,J + M^ + mj,J>2b and + mlj{\fi\' + mj,J < b\ (2.26) 

These two conditions can be satisfied if ni'jj^ 7^ ^Hd' i^plyi'^S that electroweak sym- 
metry breaking is not possible without SUSY-breaking, since in unbroken SUSY, the 
two Higgs mass terms do not exist and the masses are then equal to zero. Imposing 
the stationary conditions dV/dH^ = dV/dH^ = leads to the relations 

2 2 
+ m]j^ -bcotp- cos(2/3) = |/up + m|^^ -6tan/3 + ^ cos(2/3) = 0. (2.27) 

The particle masses are calculated by expanding the potential about the min- 
imum, once the Higgs fields get their vevs, Vu and v^, leading to various bilinear 
terms with fields with the same quantum numbers, contributing to the off-diagonal 
terms in the different mass matrices. Let us first note that the electroweak gauge 
bosons mix as in the SM. Among the eight degrees of freedom of the two Higgs dou- 
blets, three of them are the Nambu-Goldstone bosons G^*^'^^, remaining massless 
after electroweak symmetry breaking and being absorbed by the electroweak gauge 
bosons which become massive. The five remaining ones represent two CP even (the 
light ho and the heavy Hq), one CP odd (Aq) and two charged Higgs bosons (H^) 

\cosP sin (3 /\lm(//0) 



sin P — cos P\ I 
cos P sin P 




(2.29) 



r- I COS a sino \ / Ref-ff!?) - uu 1 

V2 I 1 ' (2-30) 

sin a COS a/ \Ke{H^)- 

where the mixing angle a and the tree-level masses are given by 

tan(2a) = tan(2/3) {m\ + m|) / {m\ - m|) , (2.31) 
m\ = 2 6/sin(2/3), (2.32) 
mjj± = m\ + m^, (2.33) 



1/2 + m| =F \J {m\ + m|)2 - cos(2/3)^ , (2.34) 



mz and mi^/ being the Z and W boson masses. The negatively charged Goldstone 
and Higgs bosons are defined by G~ = {G~^)* and H~ = {H~^)*. Let us note that 
the upper limit on the lightest neutral Higgs mass 

rriho <mz\cos{2p)\ (2.35) 

has already been exceeded by the current experimental lower bound of 114.4 
GeV [129j . However, the tree-level Higgs mass formulas above receive significant 
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one-loop corrections, the upper limit for m^jO being then shifted to about 140 GeV 
[TMimilT3m[T33l[TMIT35] . 



In the most general case, sfermion mass eigenstates are obtained by diagonalizing 
6x6 mass matrices, except for sneutrinos with their 3x3 mass matrix, since there 
exist only left-handed sneutrinos. In cMFV SUSY, the mixing between different 
generations, which could enhance strongly constrained FCNC processes [Ml I136t 
I137j . is neglected, and the 6x6 matrices are decomposed into several 2x2 matrices 
for squarks and charged sleptons, describing the mixing of scalars of a specific flavour 
[22lfT38] 

= [ "^^^ "^^"^^^ 1, (2.36) 
m/ mRL m%ji 



with 



tuIl = m| (r| - 6/ sin^ Ow) m% cos 2/3 m^, (2.37) 

o 

''^RR — mpi + efsin 9w itt-z ^^^'^P ~^ "i^fy (2.38) 

cot P for up — type sfermions. 
tan P for down — type sfermions. 



rriLR = rn*jiL = A} - fi < 



(2.39) 



The soft SUSY-breaking mass terms for left- and right-handed sfermions are mp 
and mp, respectively, and ^/ is the trilinear Higgs-sfermion-sfermion interaction, 
i.e. the entries of the diagonal matrices af of Eq. ()2.16p . The weak isospin quantum 
numbers are Tj = ±1/2 for left-handed and Tj = for right-handed sfermions, their 
fractional electromagnetic charges are denoted by e/, and 6w is the weak mixing 
angle. Ai"^ is diagonalized by a unitary matrix S^, Ai"^ S^^ = diag (mfjmg) and 
has the squared mass eigenvalues 



m?2 



I (mii + m%ji T ^ {ml^ - m\^Y + Amj Im^Rp) , (2.40) 



where by convention mi < m2. For real values of mm, the sfermion mixing angle 
9 J, 0<ej< -ir/2, in 



sf=\^ ^' / I with r; 1 =5/| 1 (2.41) 

can be obtained from 



smt/j cosf7j 





tan 20. = — (2.42) 

^ mi^ - mfj^ 

Let us note that for purely left-handed sneutrino eigenstates a diagonalizing matrix 
is not needed. 

The neutral Higgsinos and gauginos {B^ , W^, H^, H^) also mix when SU{2) x 
U{1) is spontaneously broken. The diagonalization of the generally complex mass 
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Y 



matrix [22l fT28] 

^ Ml — ?Ti^ sin cos P mz sin sin /? ^ 

M2 mz cos 9w cos P —mz cos 6w sin/3 

mz sin cos f3 mz cos ^vv cos P — 

y sin 6]y sin /? — mz cos 9w sin /? — /i J 

(2.43) 

by a unitary matrix leads to four neutral mass eigenstates, the neutralinos Xi 
(i = 1,2,3,4), Xi being the lightest one. The resulting diagonal matrix Njj = 
N*YN~^ = diag{m^o,m^o,m^o,m^o), has four real entries, and the neutralino 
fields are given by 

X° = %^°, i,j = l,...,4, with il^'^ = {-iB^,-iW^H^a,Hy. (2.44) 

Analytical expressions for the diagonalizing matrix and the mass eigenvalues can 
be found in Refs. |139t I140j . In general, all the possible complex phases can be 
absorbed by a redefinition of the fields, making the masses real and non-negative. 

The charged analogues of the neutralinos are the two charginos, Xi^ ^^id X2, 
resulting from the mixing of the charged wino and Higgsino fields {W^ , H^^). Two 
unitary matrices U and V are needed to diagonalize the generally complex mass 
matrix [22l ri28j 

/ M2 mw V^sinp \ 

X = \ , (2.45) 

\ mw \/2 cos P ji j 

since X 7^ X'^ . The eigenvalues of the diagonal matrix Mc = U*XV^^ = 
diag (m^ , m^ ) can be chosen to be real and non-negative, absorbing all complex 
phases by a suitable redefinition of the fields 

1 

'Xt2 2 



ml± = -\ M|+/x2+2m^^ J(M|+/x2 + 2m2^,)2-4(/iM2-m2^, sin(2/?))2 L(2.46) 



where by convention m-± < m-±, while the U and V matrices 

,0+ ifdetX>0 / cos6'± sin6'±\ 

U = O- and V = { ~ , with 0±=\ (2.47) 

asO^ if det X < ^ - sin e± cos e±J 

are determined by the mixing angles 9± with < 9± < it/2, 

2y/2mw{M2 sin P+fi cos P) 2y/2mw cos P+fi sin P) 

tan(2t^+) = ^ 7^ ——and tan(2t^_) = ^ ^ , . 

^ +^ M| -^2 + 2m2^, cos (2/3) ^ ' - /i2 _ 2m2^ cos(2/3) 

(2.48) 

The chargino mass eigenstates are given by 

U I ] and f «n = W I . 
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2.3 The MSSM with non-minimal flavour violation 



2.3.1 The model 

When SUSY is embedded in larger structures such as grand unified theories, new 
sources of flavour violation can appear [141J. In addition, SUSY parameter space is 
today becoming more and more constrained by electroweak precision measurements, 
direct searches for Higgs and SUSY particles at colliders, and precise data on the 
cosmological relic density of (possibly) SUSY dark matter. This encourages the 
investigation of non-minimal models, e.g. with additional sources of flavour violation. 
Non-minimal flavour violation in SUSY is well parameterized in the super-CKM basis 
[142j . where the up- and down- type squark mass matrices are 



Q 



^LL,l 
a21 

a31 



a21 

a31 
^LR 



Al2 


Al3 


miUlLRA 




a23 
^LL 




a32 






Al2 

^RL 


Al3 

^RL 


2 

^RR,1 


"1-2 m-i?L,2 


A 23 
^RL 


a21 
^RR 


a32 
^LR 


m3 mRL,3 


a31 
^RR 



Al2 

^LR 

•m2 mLR,2 

A 32 
^RL 



Al2 

^RR 
2 

'^RR,2 

a32 
^RR 



Al3 

^LR 

A 23 
^LR 

ms mLR,3 



Al3 

^RR 

A 23 
^RR 



m 



RR,3 



J 

(2.50) 



the indices (1,2,3) being the flavour indices {u,c,t) for the up-type and {d,s,b) for 
the down-type mass matrix, and niLL, mRR, tulr and niRL being defined by Eqs. 
(|2.37p . (|2.38p and (|2.39p . The flavour-changing elements of the mass matrices are 
usually normalized with respect to their diagonal entries 



(2.51) 



and obey to the relations ^^I/^rr = ^llrr- 
not the mass eigenstates, all the charged current and interactions couple 
with a strength given by the CKM matrix [86^ I87j . as their SUSY counterparts 
do. M? and M? are then diagonalized via two additional 6x6 matrices and 



hfi* 



In this basis, even if the fields are 



i?" M? and diag( 



2 

m- 



R^, diag(m|^ , . . . , 
by convention, the masses are ordered as rriq^ < 
eigenstates are given by 



R'^M'^^R'^\ where 
< niq^. The physical mass 



(ni,n2,tt3,n4,tt5,U6)'^ 
{di,d2,d3,di,d5,dQ)'^ 



-R"(ttL, CL, tL, UR, CR, Ir)'^, 

R'^idi, SL, II, dR, SR, IrY . 



(2.52) 



In the limit of vanishing off-diagonal parameters, the matrices R'^ become flavour- 
diagonal, leaving only the well-known helicity-mixing already present in cMFV. 



2.3.2 Experimental constraints on flavour violating SUSY 

The scaling of the flavour violating entries Ajj with the SUSY-breaking scale 
MsusY implies a hierarchy A^l > Alr,rl » ^rr |141lll43j . Note also that SU{2) 
gauge invariance relates the A|^ of up- and down-type quarks through the CKM- 
matrix, implying that a large difference between them is not allowed. Experimental 
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Table 2.2: The 95% probability bounds on \X'f'-'^'\ obtained in Ref. |144j . 





LL 


LR 


RL 


RR 


kl 










12 


1.4x10^2 


9.0x10"^ 


9.0x10-5 


9.0x10-3 


13 


9.0x10^2 


1.7x10^2 


1.7x10-2 


7.0x10-2 


23 


1.6x10-1 


4.5x10-3 


6.0x10-3 


2.2x10-1 



bounds coming from the neutral kaon sector (on Amx, e, ^V^)) ^' {^"^b) and D- 
meson oscillations (Am^)), various rare decays (BR(6 57), BR(/i — > 67), BR(t — > 
67), and BR(r — > //7)), and electric dipole moments (d„ and de) can be used to 
set constraints on non-minimal flavour mixing in the squark and slepton sectors 
[1441 1145^ 11461 1147j . As example, we show the 95% probability bounds on |A^j''^'| 
in Tab. 12.21 [144j . In our own analysis, we take implicitly into account all of the 
previously mentioned constraints by restricting ourselves to the case of only one real 
NMFV parameter, 

A ^ XtL = y^tL- (2.53) 

Let us note that in addition, direct searches of flavour violation depend on the 
possibility of flavour tagging, established experimentally only for heavy flavours, 
comforting us in our restriction to consider only mixing between the second and the 
third generations in our analysis. 

Allowed regions for this parameter are then obtained by imposing several low- 
energy electroweak precision and cosmological constraints. We start by imposing the 
branching ratio 

BR(6 ^ 57) = (3.55 ± 0.26) x 10-^, (2.54) 

obtained from the combined measurements of BaBar, Belle, and CLEO [148] , which 
affects directly the allowed squark mixing between the second and third generation. 
A second important consequence of NMFV in the MSSM is the generation of large 
splittings between squark-mass eigenvalues. The splitting within isospin doublets 
influences the Z- and M^-boson self-energies at zero-momentum S2,vk(0) in the elec- 
troweak p-parameter Ap = S^(0)/M| — Svy(0)/M^ and consequently the Vl^-boson 
mass My/ and the squared sine of the weak mixing angle sin^ Qy^. The latest com- 
bined fits of the Z-boson mass, width, pole asymmetry, W^-boson and top-quark 
mass constrain new physics contributions to T = —0.13 it 0.11 [149j or 

Ap = -aT = 0.00102 ± 0.00086, (2.55) 

where we have used a{Mz) = 1/127.918. A third observable sensitive to SUSY 
loop-contributions is the anomalous magnetic moment = (g'^ — 2)/2 of the muon, 
for which recent BNL data and the SM prediction disagree by |149j 

Aa^ = (22 ± 10) X 10-1°. (2.56) 
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For cosmological reasons, i.e. in order to have a suitable candidate for non-baryonic 
cold dark matter [30], we require the LSP to be stable, electrically neutral, and a 
colour singlet. The related dark matter relic density is constrained to the region 

0.094 < ncDAih^ < 0.136 (2.57) 

at 95% {2a) confidence level, which has been obtained from the three-year data of 
the WMAP satellite, combined with the SDSS and SNLS survey and Baryon Acous- 
tic Oscillation data and interpreted within an eleven-parameter inflationary model 
|150| . which is more general than the usual six-parameter "vanilla" concordance 
model of cosmology. 



We impose both the BR(5 — > 57) bounds on the two- loop QCD/one-loop SUSY 
calculation \li7\ I151j and the Ap bounds on the one-loop NMFV and two-loop 
cMFV SUSY calculation |152j at the 2a-level. We take into account the SM and 
SUSY contributions to up to two loops [ 1531 1154j and require them to agree with 
the experimental region within two standard deviations. i^cDM is calculated using 
a modified version of DarkSUSY 4.1 [155] . that manages NMFV. 



2.3.3 Scan of the mSUGRA parameter space 

The above experimental limits are now imposed on the minimal supergravity 
models with their five free parameters (mo, ^1/21 ^O) tan/3, sign(|u)) at the grand 
unification scale. Since our scans [97] in the (mo, mi/2) plane depend very little 
on the trilinear coupling Aq, we set it to zero. Furthermore, we fix a small (10), 
intermediate (30), and large (50) value for tan/3, and investigate the impact of 
the sign of n for tan P = 10 only, before we set it to > for the two other 
cases (see below). We solve the renormalization group equations numerically to 
two-loop order using the computer program SPheno 2.2.3 [156] . evolving the five 
parameters at grand unification scale in order to compute the soft SUSY-breaking 
parameters at the electroweak scale (see Eq. (|2.16p ) with the complete one- loop 
formulas, supplemented by two-loop contributions in the case of the neutral Higgs 
bosons and the /^-parameter. At this point we generalize the squark mass matrices, 
introducing the A parameter, and compute the low-energy, electroweak precision, 
and cosmological observables with the computer programs FeynHiggs 2.5.1 [157] 
and DarkSUSY 4.1 [T55] . 

For the masses and widths of the electroweak gauge bosons and the mass of 
the top quark, we use the current values of mz = 91.1876 GeV, mw = 80.403 
GeV, mt = 174.2 GeV, Fz = 2.4952 GeV, and Tw = 2.141 GeV. The CKM-matrix 
elements are computed using the parameterization 



^ C12C13 S12C13 si3e~'^ \ 

V = -S12C23 - Ci2S23Si3e*'' C12C23 - Sl2S23Sl3e*'' S23C13 

\ S12S23 - Cl2C23Sl3e*'^ -C12S23 - Sl2C23Sl3e*'^ C23C13 / 



(2.58) 



where Sij = sin 6ij and Cij = cos 9ij relate to the mixing of two specific generations i 
and j and 6 is the SM CP violating complex phase. The numerical values are given 
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by 

si2 = 0.2243, S23 = 0.0413, S13 = 0.0037, and <5 = 1.05. (2.59) 

The squared sine of the electroweak mixing angle sin^ 9]y = 1 — / and 
the electromagnetic fine structure constant a = \f2G pm^ 'svc? Q-s^ j tx are cal- 
culated in the improved Born approximation using the world average value of 
Gf = 1.16637 • 10"^ GeV-2 for Fermi's coupling constant [149J. 

Typical scans of the mSUGRA parameter space in mo and ^^^2 with a relatively 
small value of tan /J = 10 and tIq = are shown in Figs. 12.11 and 12.21 for /x < and 
yu > 0, respectively. All experimental limits described in the previous section are 
imposed at the 2(T-level. The b ^ sj excluded region depends strongly on flavour 
mixing, while the regions favoured by — 2 and the dark matter relic density are 
quite insensitive to variations of the A-parameter. Ap constrains the parameter 
space only for heavy universal scalar masses niQ > 2000 GeV and heavy universal 
gaugino masses mi/2 > 1500 GeV, so that the corresponding excluded regions are 
not shown here. The dominant SUSY effects in the calculation of the anomalous 
magnetic moment of the muon come from induced quantum loops of a neutralino 
or a chargino and a slepton, while squarks contribute only at the two-loop level. 
This reduces the dependence on flavour violation in the squark sector considerably. 
Furthermore, the region /i < is disfavoured in all SUSY models, since the one- loop 
SUSY contributions are approximatively given by [158| 

^SUSY, 1-ioop ^ 13 ^ 10-10 ( HL^y tan/3 sgn(/.), (2.60) 

V ^^SUSY / 

where MgusY is a typical SUSY mass scale. Negative values of p would then 
increase, not decrease, the disagreement between the experimental measurements 
and the theoretical value of a^. Furthermore, the measured 6 — > 57 branching ratio 
excludes virtually all of the region favoured by the dark matter relic density, except 
for very high scalar SUSY masses. We therefore do not consider negative values of 
fj, in the rest of this work. 

In Figs. 12.31 and 12.41 we show the {mQ,mi/2)-p^8knes for larger tan/3, namely 
tan /3 = 30 and tan /3 = 50, and for p > 0. The regions which are favoured both 
by the anomalous magnetic moment of the muon and by the cold dark matter relic 
density, and which are not excluded by the 6 — > 57 measurements, are stringently 
constrained and do not allow for large flavour violation. 



2.3.4 NMFV benchmark points and slopes 

Restricting ourselves to non-negative values of p, we now inspect the {mo, 1711/2)' 
planes in Figs. 12. 2112. 41 for scenarios that are allowed or favoured by low-energy, 
electroweak precision, and cosmological constraints, that permit non-minimal 
flavour violation among left-chiral squarks of the second and third generation up to 
A < 0.1, and that are at the same time collider- friendly, with relatively low values 
of m-o and mi/2- We propose [97j the four benchmark points given in Tab. 12.31 We 
also attach a model line (slope) to each point in Tab. 12.41 These slopes trace the 
allowed/favoured regions from lower to higher masses and can, of course, also be 
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A^=0. kw^=IO. \1<0, X=0 A, 1=0, iMi^=10, \1<0. X=0.03 A^,=0. urn ^=10. n<0,X=a05 A„=0. hiii^=JO. \1<0, 'A=O.I 




nij^ [Gevf mj,j [GeV]"*^ '"^"^ "™ nij^ [Gevf" nij^j [Gevf 



Figure 2.1: The (mo, mi/2)-pla'nes for tan/3 = 10, = GeV, ^ < 0, and A = 0, 
0.03, 0.05 and 0.1. We show WMAP (black) favoured as well as 6 ^ 57 (blue) 
and charged LSP (beige) excluded regions of mSUGRA parameter space in minimal 
(A = 0) and non-minimal (A > 0) flavour violation. 



A„=0. kw^=IO. \1>0. X=0 A, 1=0. kiiiP=IO. \1>0. X=Om A„=0. ?(»iP=/0, \i.>Q.\=0.05 A,^=0, iai!^=10. \i>0. 'k=O.I 




m,^ [GeV] m,„ [GeV] m,„ [GeV] m,^ [GeV] 



Figure 2.2: Same as Fig. EH for tan/? = 10, Aq = GeV, /i > 0, and A = 0, 0.03, 
0.05 and 0.1. 

A^=0. kiji^=JO. \1>0 ^.=0 A^=0. tan ^=30, \\>0. X=0.03 Ag=0. t<.in^=30. ji>0.X=0.05 A^=a taTi^=30, \i>0. X=0.1 




i 



Figure 2.3: Same as Fig. ED for tan /? = 30, = GeV, ^ > 0, and A = 0, 0.03, 
0.05 and 0.1. 

A „=0, P=50, M->0. A.=0 A „=0. km ^=50. \1>0. 'k=0.03 










l\ 












Figure 2.4: Same as Fig. EH for tan /? = 50, tIq = GeV, /i > 0, and A = 0, 0.03, 
0.05 and 0.1. 
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Table 2.3: Benchmark points allowing for flavour violation among the second and 
third generations for Aq = 0, fj, > 0, and three different values of tan p. 





mo [GeV] 


mi/2 [GeV] 


Ao [GeV] 


tan/? 


sgn(/i) 


A 


700 


200 





10 


1 


B 


100 


400 





10 


1 


C 


230 


590 





30 


1 


D 


600 


700 





50 


1 



Table 2.4: Model lines allowing for flavour violation among the second and third 
generations for Aq = 0, > 0, and three different values of tan/3. 



A 


180 GeV 


< 


^1/2 


< 


250 GeV , 


mo = 


-1936 GeV+ 12.9 m^/a, 


B 


400 GeV 


< 


^1/2 


< 


900 GeV , 


mo = 


4.93 GeV + 0.229 mi/2 


C 


500 GeV 


< 


^1/2 


< 


700 GeV , 


mo = 


54 GeV + 0.297 mi/2 


D 


575 GeV 


< 


mi/2 


< 


725 GeV , 


mo = 


600 GeV. 



used in MFV scenarios where the off-diagonal terms are expressed in function of the 
CKM matrix and the Yukawa couplings, and in cMFV with A = 0. 

Starting with Fig. 12.21 and tan /? = 10, the bulk region of equally low scalar 
and fermion masses is all but excluded by the 6 — > 57 branching ratio. This 
leaves us with two favoured regions. Our benchmark point A lies in the so-called 
focus point region of low fermion masses mi/2, where the lightest neutralinos are 
relatively heavy, have a significant Higgsino component, and annihilate dominantly 
into pairs of electroweak gauge bosons. Our values for the universal masses are 
smaller than those of the pre-WMAP point SPS 2 051 iB] (m-o = 1450 GeV, 
mi/2 = 300 GeV) and post-WMAP point BDEGOP E' {mo = 1530 GeV, 

mi/2 = 300 GeV), which lie outside the region favoured by (grey-shaded) and 
lead to collider-unfriendly heavy squark and gaugino masses. Our benchmark point 
B lies in the co-annihilation branch of low scalar masses mg, where the lighter fi 
mass eigenstate is not much heavier than the lightest neutralino, the two having a 
considerable co-annihilation cross section. This point differs from the points SPS 
3 (mo = 90 GeV, mi/2 = 400 GeV) and BDEGOP C' (mo = 85 GeV, mi/2 = 400 
GeV) only very little in the scalar mass. 

At the larger value of tan /? = 30 in Fig. 12.31 only the co-annihilation region 
survives the constraints coming from 6 — > 57 decays, where we choose our point C, 
which has slightly higher masses than both SPS lb (mo = 200 GeV, mi/2 — 400 
GeV) and BDEGOP V {mo = 175 GeV, mi/2 = 350 GeV). Finally for the very large 
value of tan /3 = 50 in Fig. 12.41 the bulk region reappears at relatively heavy scalar 
and fermion masses. Here, the couplings of the heavier scalar and pseudo-scalar 
Higgses and A^ to bottom quarks and tau-leptons and the charged-Higgs 
coupling to top-bottom pairs are significantly enhanced, resulting e.g. in increased 
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dark matter annihilation cross sections through s-channel Higgs-exchange into 
bottom-quark final states. So as tan/3 increases further, the so-called Higgs-funnel 
region eventually makes its appearance on the diagonal of large scalar and fermion 
masses. We choose our point D in the concentrated (bulky) region favoured 
by cosmology and at masses, that are slightly higher than those of SPS 4 
(mo = 400 GeV, mi/2 = 300 GeV) and BDEGOP L' (mo = 300 GeV, mi/2 = 450 
GeV). In this scenario, squarks and gluinos are very heavy with masses above 1 TeV. 

Let us now turn to the dependence of the precision variables discussed in Sec. 
12.3.21 on the fiavour violating parameter A in our four benchmark scenarios. As 
already mentioned, we expect the leptonic observable to depend weakly on 
the squark sector, which is confirmed by our numerical analysis of the SUSY 
contribution a^^^'^^\ We find constant values of 6, 14, 16, and 13x10"^'' for the 
benchmarks A, B, C, and D, all of which lie well within 2a (the latter three even 
within la) of the experimentally favoured range (22 it 10) x lO"^'^. 

The electroweak precision observable Ap is shown first in Figs. I2.5ti2.8l for the 
four benchmark scenarios, where only the experimental upper bound of the 2(T-range 
is visible as a dashed line. While the self-energy diagrams of the electroweak gauge 
bosons depend obviously strongly on the helicities, fiavours, and mass eigenvalues 
of the squarks in the loop, the SUSY masses in our scenarios are sufficiently small 
and the experimental error is still sufficiently large to allow for relatively large 
values of A < 0.57, 0.52, 0.38, and 0.32 for the benchmark points A, B, C, and D, 
respectively. As mentioned above, Ap conversely constrains SUSY models in cMFV 
(A = 0) only for masses above 2000 GeV for mo and 1500 GeV for mi/2- 

The next diagram in Figs. 12.5112.81 shows the dependence of the most stringent 
low-energy constraint, coming from the good agreement between the measured 
6 ^ S7 branching ratio and the two-loop SM prediction, on the NMFV parameter 
A. The dashed lines of the 2(T-bands exhibit two allowed regions, one close to 
A = (vertical green line) and a second one around A ~ 0.57, 0.75, 0.62, and 0.57, 
respectively. As is well-known, the latter are, however, disfavoured by 5 ^ s/i"^/i~ 
data constraining the sign of the 6 — > 57 amplitude to be the same as in the SM 
|160j . We will therefore limit ourselves later to the regions A < 0.05 (points A, C, 
and D) and A < 0.1 (point B) in the vicinity of cMFV. 

The 95% confidence-level (or 2a) region for the cold dark matter density is 
shown as a dashed band in the upper right part of Figs. 12.5112.81 However, only 
the lower bound (0.094) is of relevance, as the relic density falls with increasing 
A. This is not so pronounced in our model B as in our model A, where squark 
masses are light and the lightest neutralino has a sizable Higgsino-component, so 
that squark exchanges contribute significantly to the annihilation cross sections. 
For models C and D there is little sensitivity of ^cDMh'^ (except at very large A), 
as the squark masses are generally larger. The rapid fall-off of the relic density 
for very large A < 1 can be understood by looking at the resulting lightest up- 
and down-type squark mass eigenvalues in the lower left part of Figs. 12. 5112. 81 
For maximal fiavour violation, the off-diagonal squark mass matrix elements are 
of similar size as the diagonal ones, leading to one squark mass eigenvalue that 
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Up-Type squark masses 



Down-type squark masses 




Figure 2.5: Dependence of the precision variables Ap, BR(6 — > 57), and the cold 
dark matter relic density QcDM^'^ (top) as well as of the lightest SUSY particle, up- 
and down-type squark masses (bottom) on the NMFV parameter A in our bench- 
mark scenario A. The experimentally allowed ranges (within 2a) are indicated by 
horizontal dashed lines. 
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Figure 2.6: Same as Fig. 12.51 for our benchmark scenario B. 
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Figure 2.7: Same as Fig. 12.71 for our benchmark scenario C. 
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Figure 2.8: Same as Fig. 12.51 for our benchmark scenario D. 
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approaches and finally falls below the lightest neutralino (dark matter) mass. 
Light squark propagators and co-annihilation processes thus lead to a rapidly 
falling dark matter relic density and finally to cosmologically excluded NMFV 
SUSY models, since the LSP is assumed to be electrically neutral and a colour singlet. 

An interesting phenomenon of level reordering between neighbouring states can 
be observed in the lower central diagrams of Figs. I2.5ti2.8l for the two lowest mass 
eigenvalues of up-type squarks. As A and the off-diagonal entries in the mass matrix 
increase, the splitting between the lightest and highest mass eigenvalues naturally 
increases. It is first the second-lowest mass that decreases up to intermediate values 
of A = 0.2. ..0.5, whereas the lowest mass is constant, and only at this point the 
second-lowest mass becomes constant and takes approximately the value of the 
until here lowest squark mass, whereas the lowest squark mass starts to decrease 
further with A. These "avoided crossings" are a common phenomenon for Hermitian 
matrices and reminiscent of meta-stable systems in quantum mechanics. At the 
point where the two levels should cross, the corresponding squark eigenstates mix 
and change character. This level-reordering phenomenon occurs also for other 
scenarios, as well as for down-type squarks, even if it is not so pronounced for the 
latter, as shown in the lower right diagrams of Figs. 12. 5112. 81 

As shown in Eq. ()2.52p in NMFV, squarks exhibit mixing of left- and right- 
handed helicities and the three generations. For our benchmark scenario A, the 
helicity and flavour decomposition of the six up-type (left) and down-type (right) 
squark mass eigenstates is shown in Fig. 12.91 for the full range of the parameter 
A G [0; 1] and in Fig. 12.101 for the experimentally favoured range in the vicinity 
of MFV, A G [0;0.1]. Left- and right-handed flrst-generation and right-handed 
second-generation squarks remain, of course, helicity- and flavour-diagonal, with 
the left-handed and flrst-generation squarks being slightly heavier due their weak 
isospin coupling (see Eqs. (j2.37p . (j2.38p and (|2.39p ) and different renormalization- 
group running effects. The lightest up-type squark ui remains the traditional 
mixture of left- and right-handed stops over a large region of A < 0.4, but it 
shows at this point the flavour transition expected from the level reordering 
phenomenon (central plot of Fig. 12. 5p . The transition happens, however, above 
the experimental limit of A < 0.1. Below this limit, it is the states U2, uq, di, 
and in particular and dg that show, in addition to helicity mixing, the most 
interesting and smooth variation of second- and third-generation flavour content (see 
Fig. I2.10"|) . Note that at very low A ~ 0.002 the states d^ and sl rapidly switch levels. 

For the benchmark point B, whose helicity and flavour decomposition is shown 
in Fig. 12. m level reordering occurs at A ~ 0.4 for the intermediate-mass up-type 
squarks. Close inspection of Fig. 12.121 shows, however, that also dji and sji switch 
levels at low values of A ~ 0.02. At A ~ 0.01, in addition spt and switch levels, 
and at A ~ 0.002 it is the states ul and cl- The lightest up-type squark is again 
nothing but a mix of left- and right-handed stops up to A < 0.4. Phenomenologically 
smooth transitions below A < 0.1 involving taggable third-generation squarks are 
observed for u^, uq, di, and dg- 

For our scenario C, shown in Fig. 12. 13^ just below A = 0.1, ur and cr as well as 
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Figure 2.9: Dependence of the chirality (L, R) and flavour {u, c, t; d, s, and b) content 
of up- (uj) and down- type (di) squark mass eigenstates on the NMFV parameter 
A G [0; 1] for benchmark point A. 




Figure 2.10: Same as Fig.[22]for A G [0;0.1]. 
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Figure 2.11: Same as Fig. 12.91 for benchmark point B. 




Figure 2.12: Same as Fig. EH] for A G [0;0.1]. 
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Figure 2.13: Same as Fig. 12.91 for benchmark point C. 




Figure 2.14: Same as Fig. EH for A G [0;0.1]. 
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dij and sn rapidly switch levels, and ul and cl switch levels at very low A ~ 0.002. 
These changes were already visible upon close inspection of the lower central and 
right plots in Fig. 12.71 On the other hand, the lightest squarks ui and di only 
acquire significant flavour admixtures at relatively large A ~ 0.2... 0.4, whereas 
they are mostly superpositions of left- and right-handed stops and sbottoms in the 
experimentally favourable range of A < 0.1 shown in Fig. 12.141 Here, the heaviest 
uq and show already smooth admixtures of third-generation squarks as it was 
the case for the scenarios A and B discussed above. The most interesting states 
are U2, u^, d2, and d^, respectively, since they represent combinations of up to four 
different helicity and flavour states and have a significant, taggable third- generation 
flavour content. 



The helicity and flavour decomposition for our scenario D, shown in Fig. 12.151 
is rather similar to the one in scenario A, and only the mixed down-type state ^4 is 
now lighter and becomes d2- The lightest up- type squark ui is again mostly a mix 
of left- and right-handed top squarks up to A ~ 0.4, where the level reordering and 
generation mixing occurs (see lower central part of Fig. 12. 8p . At the experimentally 
favoured lower value of A < 0.1, the states U2, di, d2, and exhibit some smooth 
variations, shown in detail in Fig. 12.161 albeit to a lesser extent than in scenario A. 
At very low A ~ 0.004, it is now the up-type squarks ul and cl that rapidly switch 
levels. 



2.4 Generalized couplings 

Considering the strong interaction first, it is well known that the interaction of 
quarks, squarks, and gluinos can in general lead to flavour violation in the left- and 
right-handed sectors through non-diagonal entries in the matrices R'', 

{^qjqk'g^^qjqka] = {-^jfc' ~-^i(fc+3) } ' (^•^^) 

Of course, the involved quark and squark both have to be up- or down-type, since 
the gluino is electrically neutral. 



For the electroweak interaction, we define the square of the weak coupling con- 
stant = e^/sin^6'vy in terms of the electromagnetic fine structure constant 
a = / (47r) and the squared sine of the electroweak mixing angle xw = sin^ 6y/ = 
sl, = l- cos2 ew = l- cl,. The_Ty± - - xf , Z - xf - Xj , and Z - x° - X° 
interactions are proportional to 



_ 1 (xr T/* I 7\r _TA* „„j /^R ^ 



O^. = --=iV,4y;2 + iV.2V,l and Ol"^ = ^N*^Uj2 + N*^U,i, (2.62) 



V2 ' V2 

O'^ = -V-iV}*i - lVi2V;2 + Sijxw and O'^ = -U*^Uji - \u*2Uj2 + <5.jXh/,(2.63) 

0'l^ = -\NizN*^ + \N,^N*, and 0,'f = ^A^i^iV.s - ^iV^iV,4 • (2.64) 
In NMFV, the coupling strengths of left- and right-handed (s)quarks to the elec- 
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troweak gauge bosons are given by 

{Lqq'Z, Rqq'z} = (2 — 2 6^ Xw) X Sqq> , 



Jqq' 

3 



(2.65) 



\J-'qiqiZ-,RqiqjZ) 
{Lqq'W , Rqq'w} 



(2r| -2eqxw) X Y,{I^kRrk^Rr(3+k)Rli3+k)h (2.66) 



k=l 



{V2CwVgg:,0}, 



Y,{V2cwVu,d,I^kRjt,0}, 



(2.67) 
(2.68) 



k,l=l 



where Vki denotes the elements of the CKM-matrix. To simphfy the notation, we 
have introduced flavour indices in the latter, di = d, d2 = s, d^ = b, ui = u, 
U2 = c, and ^3 = t. The SUSY counterparts of these vertices correspond to the 
quark-squark-gaugino couplings. 



{eq-T^) swNa+T^ cos ewNi2 



2 mw cos (3 



-R 



7?* 

UjUkXi 



djUiXi 



^Ujdixf 



J?* 
UjdiXi 



j(K^o) 2mw cosp 



{eq-T^)swNa+Tj^cwN^2 



TJU* I 
^ik + 



2 mw sin /? 



eg sw Nil -RJ(fc+3) 
3 



m^, cw Rl, 



■jk 



= E 



2m^ sin/3 

'i(fc+3) 



3 

E 



mw cos (5 J 

mw sin/3 



3 r m V* R'^ 

y "^"fe ''i2 -"i(fc+3) 

^ L V^mw' sin/3 

6^ V2 cos/3 



(2.69) 
(2.70) 
(2.71) 
(2.72) 

(2.73) 
(2.74) 
(2.75) 
(2.76) 



These general expressions can be simplified by neglecting the Yukawa couplings ex- 
cept for the one of the top quark, whose mass is not small compared to mw- The 
usual MSSM couplings, and the couplings involving sleptons and sneutrinos are eas- 
ily obtained by replacing the mixing matrices, the masses, and the CKM matrix in 
a proper way from the above equations. 
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Chapter 3 



Resummation formalisms 



3.1 General points 

3.1.1 Main features of the resummation procedure 

The finite energy or angular resolution of any particle detector implies that 
physical cross sections are always inclusive over arbitrarily soft produced particles. 
For example, a single quark jet cannot be distinguished from a quark jet accom- 
panied by some collinear gluons or by partons with vanishing momentum. This 
implies that after the renormalization procedure, the remaining divergences in any 
physical observable are soft (due to parton radiation with small four-momentum) 
and /or collinear (due to emission of partons moving in parallel to the emitting one) . 

In fixed-order perturbation theory, infrared divergences of virtual gluons are 
exactly cancelled by the emission of real undetected ones, but this cancellation 
can leave large finite terms. In specific kinematics configurations, real and virtual 
contributions can be highly unbalanced, spoiling the cancellation mechanism, and 
only resummation to all order restores the balance. At the exclusive boundary 
of the phase space, when for example the tagged final state carries almost all the 
available energy or when its transverse momentum tends to zero, the real emission 
is strongly suppressed, producing a loss of balance with the virtual contribution, 
and spreads out up to the phase space kinematical limits through finite logarithmic 
terms L which become large. 

Let us take a given physical observable i?, whose perturbative expansion in the 
strong coupling constant is 



The logarithmic corrections L can be large even if Us is small, enhancing thus 
the coefficients of the perturbative expansion. The ratio of two successive terms 
in Eq. ()3.ip is O(asL^), which means that any higher order contribution is 0{1) 
with respect to the previous terms, and the reliability of perturbative predictions is 
thus spoiled. However, all-order resummation gives a satisfactory solution to this 
problem in the context of perturbation theory. 



R 




(3.1) 
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The key features of the resummation procedure are the dynamical and kine- 
matical factorizations. Dynamical factorization comes from gauge invariance 
and unitarity and is completely general. The multi-gluon radiation amplitude 
factorizes, in the soft-limit, into the product of single-gluon emission probabilities. 
It is similar to multiple soft-photon emission in QED, apart from the fact that 
radiations of photons are uncorrelated, since they are electrically neutral. In QED, 
soft divergences related to degenerate final states cancel out after summation, as 
stated by the Bloch-Nordsieck theorem [TBlj. In general, this theorem is broken in 
QCD, since the collinear degeneracy leads to a divergence that does not necessary 
cancel out in the transition rates due to colour correlation in gluon interactions. 
The Kinoshita-Lee-Nauenberg theorem \162\ I163j asserts that all the infrared 
divergences, soft and collinear, exactly cancel if the summation over initial and final 
degenerate states is carried out, and the resulting observables are infrared safe. 
Although it is not possible to sum over degenerate initial states, the corresponding 
uncancelled soft divergences are under control if the energy scale is high enough, 
because they are power-suppressed, while the QCD factorization theorem guarantees 
that the collinear divergences can be absorbed in the scale dependence of the parton 
densities, and they thus factorize from the hard-scattering process. This eventually 
leads to an exponentiation due to coherent rather than independent emission as in 



The second factorization concerns kinematics and occurs if phase space can be 
factorized in terms of single-parton phase spaces. It does not generally occur in the 
space where the shape variable is defined, because it depends in a complicated way 
on the multi-parton configuration, but rather in a conjugate space introduced via 
Mellin or Fourier transformation, in which momentum conservation is more easily 
implemented. 

Both factorizations lead to exponentiation, and the observable R defined in Eq. 
(|3.ip can be rewritten as 



where each gi function corresponds to the resummation of a specific class of loga- 
rithms and contains all QCD effects due to gluon colour charge and running coupling. 
The ratio of two successive terms is now 0{1/L), and the perturbative expansion is 
converging in the L ^> 1 region, provided that a^L <C 1. All the Os-L™ terms with 
n + 1 < m < 2n appearing in Eq. (j3.1|) are taken into account by the exponentiation 
of the lowest order terms and reappear in the expansion of the exponent at a given 
power of as- 

3.1.2 Soft-photon resummation in QED 

Let us first describe multiple soft-photon emission in QED and its exponentiation 
[55 1ll64j with a generic process fi{pi) f2{P2) F{M)+X, describing the production 



QED. 



oo n+1 




m 



n=l m=l 




(3.2) 
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q q q 
Figure 3.1: Emission of n soft photons of momenta qi. 



of a final state F of mass M, plus some unobserved states X relative to the number of 
emitted photons. The one-photon emission partonic cross section da^^^M; q), where 
q is the photon momentum, is given by the product of Born cross section do"^'^^ (M) 
and the single photon emission probability dw^^^ (q) 

dcjW(M;g) = dcT(°)(M)du;W(g) =dcj(°)(M)-^ (3.3) 

with restriction to fast fermions (with negligible mass with respect to their mo- 
mentum) and to small photon energy uj compared to the scale M. j^{q) = 
(P2/-P2 ■ Q. — Pi/pi ■ q) is the conserved emission current of a soft-photon with mo- 
mentum g by a fermion of momentum pi and antifermion of momentum p2- The 
amplitude for n soft-photon emissions of momenta qi (see Fig. 13. ip is given, in the 
soft-limit, by 

n " 

M(") = ^TjPllIl, (3.4) 
n! ^J^ qi ■ pi 

where Si is the polarization vector of the i"^ photon. By summing over the different 
polarization states, we see that the squared matrix element factorizes 

1 " 

dcT(")(M;gi,...,gJ = daW(M)-TTdt^;«(g,). (3.5) 

n! 

1=1 

This means that soft-photon emission is uncorrelated and that infrared singular- 
ities exponentiate. Indeed, kinematical factorization occurs naturally in impact- 
parameter b space, h being conjugate to the shape variable, where the (5-function 
implementing four-momentum conservation exponentiates. Summing Eq. (|3.5p to 
all orders, we get 

da(M) = da(°)(M)^y"ldu;(^)(gi)...d«;W(g„)5^^9,, + P-pi-p2^ 

= d^7W(M) I ^e'HP~Pi-P2)y\L [ dw^i)(q)e^i-> 
J (27r)4 ^ [n! J 



n 

dc7(°)(M) / _0_e*(-P-Pi-P2) exp / dw^^\q) e'"'^ 



(3.6) 



where dw is the Fourier transform of the single-photon emission probability dw and 
P the four-momentum of the final state F. By using the usual expression of dw and 
performing the integral, we can deduce an exponentiated form factor, containing all 
the logarithms. 
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Figure 3.2: Emission of two soft-gluons. 



3.1.3 Soft-gluon resummation in QCD 

Contrary to photons which are electrically neutral, gluons carry colour charge. 
However, the infrared divergences related to gluon correlations cancel among them- 
selves and an exponentiation is still possible [1651 116b'] . Considering the same generic 
process as in the previous section, the cross section for single-gluon emission can be 
written as in Eq. (|3.3p . provided that we introduce colour operators, 

dcT(i)(M;g) =dcJ(°)(M)d^l;«(g) =dcj(°)(^)7^ VJ^^Q)? , (3-7) 

where u (q) denotes the gluon energy (momentum) and the current J^{q) = Tp^jf^{q) 
includes now the colour charge operator T^^- 

For the emission of n soft-gluons of momentum qi and energy LOi , diagrams where 
gluons are emitted by harder ones have to be taken into account, as in Fig. 13.21 for a 
two-gluon emission. The difference with QED comes from the second diagram, not 
present in the case of photons. If we assume strong ordering in gluon energies uJi <C 
UJ2 <^ ■ ■ ■ <^ uJn ^ M, the use of the soft limit is justified, and we can recursively 
define the amplitude M^"') , since the emission of the softest gluon factorizes 

(ai,...,a„ I m(")) = (a2,...,a„ | J^Hq^) \ 6i, . . . , (6i, . . . , | M^""!)), (3.8) 

where we sum over repeated indices and where (^Uj) are the colour (Lorentz) 
indices of the emitted gluons. M^"^ can then be computed through an iterative 
insertion scheme, and performing the summation over n leads to the exponentiation 
of the gluonic radiations. As an example, we will show the exponentiation of the 
gluon coherent state operator, whose application on the vacuum leads to the usual 
resummation formulas, given e.g. in |167j for transverse-momentum resummation, 
in [Blj for threshold resummation, and in [75| for joint resummation in the case of 
Drell-Yan lepton pair production. 



Let us consider the emission of two gluons. Using the usual definition of the 
current J^{qi) = Z]"=2'^i^^ + ToJ/i(g) where the index labels the hard /1/2 
contribution, we get 

(ai, a2 I M(2)) = (a2 | T^^ e, • j(9i) I bi){h \ M^'^) + {a^ \ | 61) (61 1 M^)) 

= Cllfll^+Cl^;;\ (3.9) 

after the introduction of the polarization vectors in Eq. (|3.8|) . The correlation oper- 
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ators C are recursively defined by 

Cll^ = {ai\M^'^) = {b\TS'e,-j{q)\b) 

n 

Cl^-- = (a2...a„|5]r,-^|6i...6„_i)C^(t-?r. (3-10) 

i=2 

The first term of Eq. (j3.9p corresponds to the emission of the softest gluon by one 
of the external legs (left diagram of Fig. 13. 2p , and the second to the emission by the 
harder gluon (right diagram of Fig. 13. 2p . Generalization to any number of gluons is 
straightforward , 

n 

(ai,...,a„|M("))=J^ q„,)e„,...q„,)e„,e(l,...,A;)+perm., (3.11) 

fc=l {ni,...,nfe} 

where the dependence on the colour indices is not explicitly shown in order to 
simplify the expression. For each value of k, the indices {oi, . . . , an} are partitioned 
into k classes of rii indices, corresponding to k chains of rij gluons. Inside a given 
class, the gluons are ordered in energy through the 0„- functions so that a specific 
gluon can only be emitted by a harder one. Eventually, the 0(1,... ,k) function 
orders the energies of the hardest gluon of each class. 

Let us now define the coherent state operator for real emission, containing the 
correlation structure between the emitted gluons 

oo „ 

U'\A+) = Y.Ur{au...,an\M^-^)At:-{q,)...Aty{qn), (3.12) 

where = d[gi] . . .d[g„,] is the relativistic n-gluon phase space and 74^"'(g) the 

creation operator of a gluon with colour (polarization) index {oii). Eqs. (|3.10p 
and (IXTTI) lead to 

oo 71 „ 

71=1 k=i {m} 

X WC(^n,-i)Qn, [ei ■ Jmi)] . . . G(l, ...,k), (3.13) 

where is the product of rij gluon creation operators and J^V. = yZl-o T,-^— . 

After summing over {n,}, each square bracket factor reproduces the gluon coherent 
state operator U'^^{A~^) where Qi is the momentum of the fastest gluon in each class. 
Eventually, we get exponentiation 

00 „ 

Um{A+) = l + J^^s' / d[gi]...dfe] 9(1,..., fc) (£"'=• J^y...(e"i.4}i) 

k=i 

= exp 1 d[q] Q{M - w) A+^iq) [W{A+)] (e° • 

(3.14) 
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where Pw orders the colour matrices appearing in Jmq in such a way that the harder 
operator acts first. 

Eq. (j3.14p is not unitary, since we have only creation operators. The unitary 
is restored through the introduction of destruction operators, replacing A'^{q) by 
T^aiq) = (9) — Aa{q). It corresponds to the inclusion of the virtual corrections 
that we have neglected in the first place. 

3.2 Transverse-momentum resummation 

3.2.1 The Collins-Soper-Sterman (CSS) formalism 

Let us consider the process 

hi{pi)h2{p2)^F{M^,q^) + X, (3.15) 

where two hadrons hi and /12 with momenta pi and p2 collide with a centre-of-mass 
energy to produce an observed colourless final state system F with a mass 
M and a transverse momentum qt, while X is unobserved. As F does not carry 
any colour charge, the lowest order partonic mechanism is initiated either by qq 
annihilation or by gluon fusion. Isolating the divergences from the fixed-order 
production cross section, we can derive the resummed cross section to get a finite 
expression for the (/T-spectrum of the system F. 

On the basis of the QCD factorization theorem, the most general form of the 
differential cross section. 



dV 



dMMq^ 



Y] [ dXa f dXb fa/hiiXa, fJ-p) fb/hii^b; f^p) d-ab{M,qT,S; flR, flp), (3.16) 
„ . JO Jo 



where s = XaXb Sh is the partonic centre-of-mass energy, can be computed by con- 
voluting the hard scattering function with universal parton densities fa/hi ^-nd 
fb/h2 of partons a,b in the hadrons /ii,/i2. The PDFs depend on the longitudinal 
momentum fractions of the two partons Xa^b and on the unphysical factorization 
scale np. The hard part dab has a perturbative expansion in asi^R)-, fJ-R being the 
unphysical renormalization scale, 

aab{M,qT,s]HR,iiF) = '^{— — al^\M,qT,s;^iF), (3.17) 

n=0 V ^ / 

where p is the minimum power of ag needed for the process to occur at leading 

(n) 

order. Isolating the divergent terms as the transverse momentum qx —>■ 0, a^^, can 
be written as [52] 

a^\M,qT,s;f,F) = R^^^ {M, qr, s; ^xp) + ^'^^''^M, s; fip) 6{qT) 

+ Z^<6 iM,s;fiF) 2 ' (3-18) 

m=0 

(n) 

where R^f^ is the regular part of the hard function, containing all the terms less 
singular than S{q^) and as qT — > 0. All the logarithmically enhanced and the 
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(5(g|n) contributions are resummed to all orders in a^, and the remaining terms 
correspond to the usual perturbative series after the removal of the logarithmic and 
^{1t) contributions. 

We define then the finite component by 



d^a 



fin 



V / AXa I dXbfa/hi{Xa;f^F) fb/h2i^b■,^J'F) 
a,b 

xf^f^^Y^' Rab{M,qT,S■,^.F), 
n=l V / 



(3.19) 



which is exactly the difference between the full perturbative result and its asymp- 
tote at low qt, containing the terms which are at least as singular as or 
proportional to S{q^). At low qx, we are subtracting two terms dominated by their 
singularities, and their difference is thus not significant. In the large-gy region, 
the asymptote does not play a significant role, and the fixed-order theory is recovered. 



For the resummed component, we rather work with its Fourier transform W with 
respect to the impact-parameter b, the variable conjugate to qx, because as stated 
in the previous section, the kinematic factorization of the multiple gluon emission 
is more naturally performed in impact-parameter space |53i . Using the singular 
part of the partonic cross section aab, i-e. the two last terms of Eq. ()3.18p . and Eq. 
(|3.16p . we get 



oo 



a,b n=0 







W''ib;M,s;^R,^XF) 

/o 

-M/nj \xf ^ , ,ln"(MVg^ 

m=0 



2 \ -\ 



(3.20) 



The dependence in the two partons a and b factorizes, and W obeys an evolution 
equation where b and M dependences are separated, so that Eq. (|3.20|) can be sim- 
plified to the usual CSS resummation formula [56| Wl\ ^59], taking into account the 
coherent gluon emission described in section r3.1.1l [167J. Omitting the dependences 
on the unphysical scales, we get 



dV 



a, 6 



W[,ib;M,s) 



xW^,{b;M,s), (3.21) 
E /' — /' — C^aMCs/b), z,) Cg{as{C,/b), Z2) S^{M, b) 

~^J0 ^1 Jo ^2 



xa^^°)^(M)5(M2 



Zl Z2 S) 



(3.22) 



where Jq{x) is the zeroth order Bessel function of the first kind. By convention, the 
constants Ci, C2, and C3 are chosen to be Ci = C3 = 60 = 2e~^^ and C2 = 1, where 
= 0.5772... is the Euler-Mascheroni constant. is computed perturbatively 

and is the process-dependent partonic cross section. It embodies the all-order re- 
summation of the large logarithms ln(M^ b"^) (conjugate to InM"^ /q^), and Ucg*"^^^ is 
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the total cross section for the LO partonic subprocess cc F. The resummation of 
the logarithms is completely achieved by the exponentiation in the quark (c = q) or 
gluon (c = g) Sudakov form factor , which can be expressed as [551 EH EZl I168j . 



5f (M, b) = exp 




r,2 



,(3.23) 



depending only on the two functions Ac, relative to soft-radiation, and relative 
to flavour conserving collinear radiation, and which can be perturbatively computed 



AJa. 



E 

n=l 



and 



S(^^i)"B<»). (3.24) 



n=l 



The coefficients of these series can be obtained by comparing the expansion of Eq. 
()3.2ip at a given order in ag with the expression of the fixed-order cross section in 
the small-g^^ limit. The lowest order coefficients, needed for resummation at NLL 
accuracy, have been computed both for the quark [551 1168j and the gluon |169j form 
factors, 



^(1) 
g 



Af)- 



Cf 
-CfK 



-Cf 



and 


^(1) 
g 


and 


Af 


and 





Ca 
1 



CaK 



(llC7A-2n/), 



(3.25) 
(3.26) 

(3.27) 



where K 



~ IT ) C'a 



, and where Cp = 4/3 and Ca = 3 are the usual 

(1\ F 

QCD colour factors. The process-dependent coefficient Be has been computed 
for different processes, as for example Drell-Yan pair [60] or Higgs production |170j . 

Finally, the coefficient functions of Eq. ()3.22p can as well be computed per- 
turbatively 



Cab{as,z) = Jab 5(1 -2) + ^ 



n=l 



TT 



(3.28) 



where z = M'^/s. They contain a collinear contribution which has its origin in the 
particularities of the MS scheme, where the full splitting functions are not factorized 
into the parton densities, and a hard process-dependent contribution coming from 
the finite part of virtual-loop corrections. The first coefficient C^^-* is known for a 
large number of processes [5S1 [Ml lIZOl EUl JTSl JTSl IMl 175] . 



3.2.2 Non-perturbative effects in the CSS formalism 

The gT^-distribution is affected by non-perturbative effects associated with the 
large-5 region, b > 1/Aqcd, where the strong coupling and the PDFs at the 
scale bo/b enter the non-perturbative regime. To take these effects into account, 
the function is evaluated at a new variable b^^ = b/ y/l + (6/^max)^, far from 
the region where non-perturbative effects are relevant, the latter being included in a 
function F^^ . Practically, the function of Eq. (|3.21|) is replaced by 

wS,''^{b;M,xuX2) = W^,{b*;M,s)FS^{b;M,xuX2). (3.29) 
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The value of 6max and the shape of F^^ have to be chosen so that 
W^jj^{b; M, xi,X2) — Wab{b; M, s) when b < 6max- In the original formahsm [561 EH] 

Pab'{b,M,xi,X2) = exp [-ln{M^bl,^Jgi{b) - ga/hA^u^) - 56/^2 (^2, &)] , (3.30) 

where the three g-functions are assumed to vanish as 6 ^ 0. The predictive power of 
the CSS formalism relies on the universality and the scale dependence of these func- 
tions. As the PDFs, they can be evolved from a given energy until the required scale. 

The main source of non-perturbative effects comes from partons with a non-zero 
intrinsic transverse-momentum already inside the hadron and from unresolved gluons 
with qT < l/&max- Global fits of experimental Drell-Yan data allow for different forms 
for this non-perturbative function [6OI I176|. I177|. I178|. I179j 



FT^''^''\b,M,x,,x,)=e^V 



fJ('^)(6,M,xi,X2) =exp 
fJ('''"^)(6,M,xi,X2) =exp 
fJ('=")(6,M,xi,X2) =exp 



b 51 + 52 In ■ 



-b'^{ 91+92 In ^"""^^ ]_5^^g,3 ln(100xi2;2; 



-b^[ 91 + 92 In ^"""^^ + 91 93 ln(100 X1X2] 



6 I ai + 02 In ■ 



2 
M 



+ 03 ln(100 xi X2] 



3.2GeV 

The most recent values of these parameters can be found in [1781 1179] . 



(3.31) 
,(3.32) 
,(3.33) 
.(3.34) 



3.2.3 Disadvantages of the CSS formalism 

Although the Collins-Soper-Sterman resummation formalism has been used in 
several processes at various levels of perturbative accuracy, it presents some disad- 
vantages, in particular regarding the process dependence of the various coefficients 
of Eqs. ([3211 [3221 [323]). For instance, the first coefficient of the C-functions, C^J^^ 
can be written in terms of the one-loop matrix element of the considered process 
|173j . Even the Sudakov form factor, supposed to be universal and to depend only 
on the quark or gluon nature of the emitting particles, is process-dependent through 
the coefficients involving also loop diagrams |170j . Moreover, the dependence 
on renormalization and factorization scales, parameterized through the arbitrary 
coefficients Ci, C2 and C3 does not correspond to the usual procedure followed in 
fixed-order perturbative calculations, and the parton distribution functions called at 
the scale bo/b rather than fxp involve extrapolation in the non-perturbative region. 
Eventually, expanding the resummed expression at a given power in as in the 
large-gr region leads to unwanted factorially growing coefficients, with oscillating 
signs [TgO] . 

Some of the difficulties can be overcome by performing resummation in g^-space 
instead of 6-space |18H I182| . but at the cost of non fulfilment of the transverse- 
momentum conservation [53]. However, all these problems can be avoided by using 
the recently proposed universal resummation formalism in 6-space |6H [52] 1183] . 
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3.2.4 Universal resummation formalism 



General expressions of the first C-function coefficient and of the second B- 
function coefficient show exphcitly the process-dependence of these terms |170j 



Ci^)^(z) = -P^^{z) + 6ab6{l-z){Ca^ + ^Aa{ct>)], (3.35) 



1 

6^2' 

= -2 5P(2)+/3o0C7„vr2 + X(0)), (3.36) 
where -P^fc(-z) is the C'(e)-term in the one- loop Altarelli-Parisi splitting kernel 

(2) 

(e = (4 — D)/2), 6Paa the coefficient of the 6{1 — z) term in the two-loop splitting 
functions, and Aaifj)) is the finite part of the one-loop virtual contributions. Contrary 
to the universal coefficients A^'^ and B^^^ determined only by the Altarelli-Parisi 
splitting functions, Ba and C^j, contain both collinear and hard contributions, 
for which should be evaluated at two different scales, the same scale as the 
parton densities bo/b for the process- independent collinear contribution and the 
hard scale Af for the process-dependent hard contribution. 

Making the replacement 

)^(M) ^ cT^f )^(M, a,(M)) = 4^°^^(M) i/f («.(M)), (3.37) 

we include now all the hard process-dependent contributions in cx^g^^^ through the 
function , which has a perturbative expansion 

00 

n=l 

Eqs. (j3.22p and (j3.23p can then be rewritten in their universal form, 
W[,{b;M,s) = V [' — Cca{as{bo/b),zi)C,b{as{bo/b),Z2)S,{M,b) 

^ Jo ^1 Jo ^2 

xa^g°)-^(M) 6{M^ - zi Z2 s), (3.39) 



Sc{M,b) = exp J - / 



2 



M2 

A^{a,{q)) In ^ + B^ia.^q)) 



(3.40) 



where the subscript F denotes the process-dependent factors. The universal Sudakov 
form factor Sc contains only real and virtual contributions due to soft and flavour 
conserving collinear radiation at scales M > qt ^ 1/b, the Cab coefficients include 
real and virtual contributions due to collinear radiations at very low qt ^ 1/6, the 
hard contributions produced by virtual corrections at scales qx ^ M are embodied 
in , and eventually, soft radiations at very low qt < 1/6 are not considered, since 
real and virtual soft contributions cancel in this kinematic range due to infrared 
safety. All of these coefficients depend now only on the flavour and the colour 
charges of the radiating partons. 
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The two versions of the resummation formula can be related by the use of the 
renormalization-group identity 



H^{as{M)) = exp 



H^{asibo/b)), (3.41) 



h^iasiM)) = /?(«.)^^^^^, (3.42) 

where P{as) is the QCD /3-function. The relations between the process-dependent 
and process-independent coefficients are then given by 



C^bias,z) = ^HF{as)Cab{as,z) (3.43) 

i?f(a.) = i?,(a,)-/3(a,)iM^. (3.44) 

dmos 

This identity can be use to show that the universal resummation formula is invariant 
under the transformation 

H^{as{M)) ^ H^{as{M))[g{as{M))]-\ (3.45) 

where g is an arbitrary perturbative function. The universal resummation coeffi- 
cients are then not unambiguously determined, which is a consequence of the fact 
that the transverse-momentum distribution is not a collinear-safe observable. This 
ambiguity is similar to the one encountered in the parton distribution functions, 
renormalized by fixing an arbitrary factorization scheme. The resummation coeffi- 
cients then have to be defined after choosing a specific resummation scheme. In this 
work, we set {as) = H^^ {as) = 1, corresponding to the Drell-Yan resummation 
scheme. 

In Eq. (j3.2ip . the PDFs are called at a scale embodying a 6-dependence, which 
eventually leads to an extrapolation into the non-perturbative regime. Besides, the 
dependence on the factorization scale cannot be directly computed as for the 
fixed-order theory. These problems can be solved |184[ I185t I186t I187j by using the 
scale dependence relation for parton densities 

fa/h{xM/b) = ^ \ —Uad{z]bo/b,iJ,F)fd/h{x/z,iJ,F), (3.46) 

d ''^ ^ 

where Uab is the evolution operator matrix obtained by solving the DGLAP evolution 
equations to the required perturbative accuracy. To avoid dealing with convolution 
integrals, we rather express the resummation formula in Mellin A/^-space, taking 
moments with respect to r = M"^ / Sh- Eq. (|3.2ip can then be rewritten, restoring 
the dependence on the factorization scale, as 



dVjv 
dM2 d4 



res ab 



f-OO I 

fa/h, {N+l; fip) fb/h, {N+l; /iF) / db - Mb qT)W[,{N, 6; M, ^^F). 

ab 

(3.47) 
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The moments of the PF-function are given by [61j 
W[,iN,b;M,t,F) = ^a^J^''^^iM)H^{asiM))S,iM,b)Y,(Ccd{N;aMb)) 



c de 

X Uda{NM/h,lJ^F)C-,e{N-as{bo/b))Ueb{N-M/b,fiF)). (3.48) 

where we consider only a single parton species to simplify. The generalization to the 
multiflavour case is treated in App. A of [02]. This expression can be written in a 
resummed form, where constant and logarithmic terms factorize, after restoring all 
scale dependences, 



tbi^^b;M, iiF,I^R) = 'Hab(^^ OisinR)] — ,— ,— )exp|a(iV,L,as(/xij); — ,77)), 

(3.49) 

with L = InQ^^^/^Q. The resummation scale Q is introduced due to the degree of 
arbitrariness involved by the factorization [188J, since the argument of the large 
logarithms can be rescaled as ln(M^6^) = ln((5^6^) + ln(M^/(5^), provided that Q is 
independent of b and ln(M^/Q^) = 0{1) when 6M S> 1. Similarly to the case of fiR 
and ^F, one should set Q = M and estimate uncalculated subleading logarithmic 
corrections by varying Q around the central value M. 

The function does not depend on the impact-parameter b and therefore 
contains all the perturbative terms that behave as constants in the limit 6 — > 00. Its 
evaluation can be done perturbatively 



Ff^r M M M ^ (LO)F.„. 



^\ IT J fiF Q 

(3.50) 



Let us remark that the notation ab —> cc is a compact way of writing the full 
subprocess ab — > cc+X F+X. The first coefficients for Drell-Yan pair production, 
i.e. those corresponding to the emission of one gluon, Tif^^gq, and of one quark, 
n'„^ , are given by [62J 



rtgg^gg l^V , , , ^ J Hg ^lig + ^ Jig m ^2 y Q2 



+2C«(7V) + 2 7S)(iV)ln^, (3.51) 

I^F 

= CgHN) + i^iN)ln9l, (3.52) 

flR fXF Q fJ-F 



with 



^2 



The exponent Qj^ includes all the logarithmically divergent terms when 6 — > cxd and 
can be systematically expanded as 
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g{N, L; as, — = L (i L) + g^^^ {N, - /3o a,(/i|) L; — + . . . , 

(3.54) 

where Po = {UCa — 2nj)/12 is the first coefficient of the QCD (3 function, Uf 
being the number of flavours. The term Lg^^'^ cohects the leading logarithmic (LL) 
contributions, the function g^'^^ the NLL ones, and so forth. The explicit expressions 
for the gi functions needed to perform NLL resummation are given by 



.(2) 



5«(A) 
iV,A;— ,— 

m Q J 



^(1) A + ln(l-A) 
To A ■ 



(3.55) 
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+ 



ln(l - A) 



Pi 
/3o 



-ln2(l 
2 ^ 
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Pi 
ln(l 



1 - 
A) 



1 - A 



A 



+ 



+ ln(l - A) 



A 



1 - A 



1 - A 



ln(l- A) ) In^, 



(3.56) 



where /3i = {11 C\ — bCAnj — 3Cjnj)/24 is the second coefficient of the QCD (5 
function. In the small-6 (large-gr) region, resummation should not be applied since 
the perturbation theory is reliable. A slight modification of the expansion parameter 
L is then introduced [1831 1185j , 



L ^ L = In 



2 u2 



+ 1 



(3.57) 



so that the logarithmic terms are now suppressed for small 6-values, reducing the 
impact of the unjustified resummed logarithms in the small-6 region, but leading to 
an equivalent behaviour in the large-6 region. 



3.2.5 Inverse transforms and matching procedure 

Once resummation has been achieved in N- and 6-space, inverse transforms have 
to be performed in order to get back to the physical x- and gr-space. Special 
attention has to be paid to the singularities in the resummed exponent, related to 
the divergent behaviour near A = and A = 1. For A^-space, an inverse Mellin 
transform is performed following a contour inspired by the minimal prescription 
[189| and the principal value resummation |190j . 



dV 



dM2 dql 



dN 

2TTi 



-N 



dM2 dq^ 



(3.58) 



where the contour Cat is chosen in such a way that all the singularities related to the 
A^-moments of the PDFs are to the left of the integration contour in the complex 
A-plane 



AT = C + ze^''f', 



(3.59) 



with < z < oo, vr > 
integral. 



> tt/2 and C > 0. This leads to an exponentially convergent 
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For the inverse transform of 6-space, we should first note that the functions g)^ 
are singular when = {b'^/Q'^ exp{7r/(/3oas)}, which is related to the divergent 
behaviour of the running perturbative coupling near the Landau pole, corresponding 
to very large values of b. In Sec. 13.2.21 these singularities were regularized through 
the variable 6*, preventing the variable b to become too large. Alternatively, we can 
deform the integration contour of the b integral in Eq. (j3.47p . continuing it in the 
complex plane [76 ^ [77] . We define then two integration branches 



b = (cos ip ±1 sin ip)t 



(3.60) 



with < t < oo, and the Bessel function is replaced by the auxiliary functions 



hi{z,v) 



1 

vr 
1 

vr 



-7v~\-ivn 



ivn 

— iVTT 



(3.61) 



Their sum is always hi{z,v) + h2{z,v) = 2Jq{z), but they distinguish positive and 
negative phases of the 5-contour, being then associated with only one of the branches. 

The parameter C is chosen in such a way that all the singularities related to the 
A^- moments of the parton densities are to the left of the integration contour. It has 
to lie within the range < C < exp[7r/(26oas) — 7e] in order to have convergent 
inverse transform integrals for any choice of (f) and 



The resummed component of the cross section dominates in the smsdl-qT region, 
and the finite component defined in Eq. (|3.19|) dominates at large values of qt- In the 
intermediate-gT region, both components have to be consistently matched in order 
to obtain uniformly accurate theoretical predictions. To this aim, we express the 
finite component as the difference between the usual fixed-order result at a specific 
order f.o. in as and the expansion of the resummed component at the same order 



r dV 1 




r dV 1 




r dV 1 




dM2d4 


fin 


dM2d4 


f.o. 


dM2d4 


res 



(3.62) 



f.o. 



Besides we impose the condition 



r dV 1 






" d^a " 








reSjl.a. 


f.o. 


dM2dg2 


res 


f.o. 



(3.63) 



meaning that the expansion at a given order f.o. of the resummed component evalu- 
ated at a specific logarithmic accuracy La. equals the expansion of the full resummed 
component at the same order f.o. in Og- The complete cross section at a given loga- 
rithmic accuracy is then given by 



dV 



dM2dg2 



d2a 



dM2dg2 



+ 



resj.a. 



d'a 



dM2dg2 



(3.64) 



fin 



Our matching procedure guarantees that we retain the full information of 
the perturbative calculation up to the specified order, plus the resummation of 
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logarithmically enhanced contributions from higher orders, without double-counting 
any term. Moreover, after integration over qx, it allows us to exactly reproduce 
the fixed-order calculation of the total cross section, which is not the case for the 
CSS-formalism due to non-vanishing contributions of the resummed component in 
the high-QT region. 

The fixed-order truncation of the resummed component [d^c]^^^ |f o. is obtained 
by perturbatively expanding the VF-coefRcient in Eq. (|3.39p and performing the b- 
integral of Eq. ([337]), 



dVgfe 

dM2dg2 



exp c L n=l 

~(n) ( ~ M M M 



TT 



(3.65) 



where ab = qq, qg. Considering Drell-Yan lepton pair production, we get for the NLL 
resummation matched with the O(a^) fixed-order calculations 162J 



~ M 



Q 



with 



F(l;2) 
qq^qq 



(N) 



E 

■t=i 



M 



and (iV) = 0, 



(3.66) 



(3.67) 



/ N- — 

m^m \ Q 



'qq 



and (iV) 



-7£\(3.68) 



In{QT/Q)=Q' db-UhqT)W^ 
Jo ^ 



1 



(3.69) 



3.3 Threshold resummation 

3.3.1 Formalism 

Integrating Eq. (|3.16p over qT, the differential cross section can be written as 

dcj 



/ dXa dXb fa/h:^iXa;iJ'F) fb/hiiXblfJ'F) ^abiz,M;asiHR),llR,I^F), (3.70) 



dM2 



where the partonic cross section aab is expanded in powers of a 

M M 



At the n*^ order, the mismatch between virtual corrections and phase-space 
suppressed real-gluon emission leads to the appearance of large logarithmic terms 
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ag[ln^'^~^(l — z)/{l — 2;)]+, which have to be resummed when s is close to M^. 
Although these large logarithms are manifest in the partonic cross section, they 
do not generally result in divergences in the physical cross section, contrary to the 
QT-spectrum, because they are smoothed by the convolution with the steeply falling 
parton distributions in Eq. (j3.7U|) . Threshold resummation is then not a summation 
of kinematic logarithms in the physical cross section, but rather an attempt to 
quantify the effect of a well-defined set of corrections to all orders, which can be 
significant even if the hadronic threshold is far from being reached. 

The hadronic cross section of Eq. ()3.70p is more conveniently written in Mellin 



A^-space [191] . 

^-^ MM 

a{N, M) = V f,/h. {N + 1, ^1F) f,/h,{N + 1, fip) ^ab{N, M; as, — , — ). (3.72) 

After performing the resummation of the radiative corrections, the moments of the 
partonic cross section are given by 

air)(A^, as) = C,fe(a,) exp [s{N, as)] , (3.73) 

where the scale dependences are suppressed for brevity and where a^^^^ represents 
the Born cross section. In Mellin space, the Caf,-functions, collecting the hard con- 
tributions, i.e. the A^-independent terms in Mellin space or the terms proportional 
to 5(1 — z) in the physical space, can be written as a perturbative series in the strong 
coupling, 



00 

n 



n=l 



Cabias) = Sab + Y.(^) ^i?- (3-74) 



For Drell-Yan pair production, the first coefficients are given by 
cS^ = ^^^(^-4+^ln^) and 0^=0. 



(3.75) 



The universal Sudakov form factor S is given by integrals over functions of the 
running coupling. 



^2 



S{N,as) = 2 dz— / J^A{as{q^)). (3.76) 

^0 J- - ^ Jul Q 

This form is valid for processes where the final state particles do not carry any colour 
charge. The function A embodies the contributions related to the collinear emission 
of soft gluons by initial-state partons. It is a series expansion in the strong coupling 
constant, 

00 

^(^s) = Y.{^y An, (3.77) 

n=l 

whose coefficients are perturbatively computable through a fixed-order calculation. 
In particular, it has been proven [192j that in the MS factorization scheme the 
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coefficients A„ are exactly equal to the large- coefficients of the diagonal splitting 
function 



/' 

Jo 



-A{as) lniV + 0(l) 



(3.78) 



where N = N exp[7£;]. Performing the integration in Eq. ()3.76p and using Eq. ()3.77p . 
we obtain the form factor up to NLL accuracy, 



S{N,as)=gi{X) \nN + g2{X). 



(3.79) 



The functions gi and g2 resum the LL (a^ ln"+i N) and NLL (a^ In" N) contribu- 
tions, respectively, and are given by [661 ED] 



51 (A) 

52 (A) 
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[2A-F(l-2A)ln(l-2A)], 

2\ + ln(l - 2 A) ^ ln2(l - 2 A) 



(3.80) 



^(2) 

— 2-[2A + ln(l-2A)] 



, , , ^2 2A« A M2 

+ -^[2A + ln(l-2A)]ln^ In^, 

Po 



/3o 



(3.81) 



where A = [/3o c^s InA^J/vr. Thus, the knowledge of the first two coefficients of the 
function A{as) [H5lfT69]. 



^(1) = Cf and ^(2) = ^ '^^ 



Ca 



67_ ^ 
18 6 



9 



■N 



(3.82) 



together with the first coefficients of the C-functions allows us to perform resumma- 
tion up to NLL. 



3.3.2 Improvements of the resummation formalism 

In the limit of large N ^ the cross section is clearly dominated by terms of 
0{\.v? N)^ O(lnA^) and 0{1). It seems thus reasonable to neglect terms suppressed 
by powers of 1/A^ in the resummation formalism. Actually these last terms are mul- 
tiplied by powers of \n.N and could as well provide a non-negligible effect in the 
threshold limit. In [681 [69] it has been shown that these contributions are due to 
collinear parton emission and can be consistently included in the resummation for- 
mula, leading to a "collinear-improved" resummation formalism. The modification 
simply amounts to the introduction of an A^-dependent term in the and C^g^ 
coefficients. For the Drell-Yan case, Eq. ()3.75p is replaced by 
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■In 
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7^ 
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(5(1) = (-(1) _ 
qg aa ^ 



IniV- iln^ 



qg 
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(3.83) 



(3.84) 



Furthermore, for Drell-Yan processes and deep inelastic scattering, the exponen- 
tiation of the contributions embodied in the C-function has been proven in |193j . 
leading to the following modification in Eq. (|3.73p : 



'^ir^(A^,«s) = ^^^°^ exp 



exp 



S{N, as 



(3.85) 
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As the authors of Ref. [193j recognize, this exponentiation of the A^-independent 
terms is not comparable to the standard threshold resummation in terms of predictive 
power. While in the latter case a low-order calculation can be used to predict 
the behaviour of full towers of logarithms, in the former case it is not possible to 
directly get information on the behaviour of constant terms at, say, n loops, but a 
complete calculation at the n^^ perturbative order is still necessary. Nonetheless, the 
comparison of the numerical results obtained with and without the exponentiation 
of the constant terms can at least provide an estimate of the errors due to missing 
higher-order corrections. 



3.3.3 Inverse Mellin transform and matching procedure 

As for transverse-momentum resummation, once resummation has been achieved 
in Mellin space, an inverse transform has to be performed in order to get back to 
the physical x-space. The customary way to perform this inversion, avoiding the 
singularities of the A-moments, is the "Minimal Prescription" of |189j . 



1 



2TTi 



dNr-^ a{N,M). 



(3. 



The constant Cmp has to be chosen so that all the poles in the integrand are to the 
left of the integration contour in the complex A-plane except for the Landau pole 
at A = exp[7r/(2 Pq which should lie far to the right on the real axis. 

Finally, a matching procedure of the NLL resummed cross section to the NLO 
result has to be performed in order to keep the full information contained in the fixed- 
order calculation and to avoid possible double-counting of the logarithmic enhanced 
contributions. A correct matching is achieved through 
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(F.O.: 



+ 
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2TTi 



dAr 



-N 



a 



('^'=^)(A,M) -a(^^P)(A,M) 



(3.87) 



where 

^(F.O.) is 

the fixed-order perturbative result, the resummed cross 

section, and the truncation of the resummed cross section to the same 

perturbative order as o"^^'^'^. In the Drell-Yan case, and taking into account the 
improvement of Eqs. (|3.83p and ()3.84p . the expansion of the resummed partonic 
cross section up to order reads 



a(-P)(A,M) 
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1 + ^ [Cp[2hx^N 



a, 



a.. 



L vr 



2 InA^ 



M_2_\ Ml) 
2 j + ^gg 



, (3.88) 
(3.89) 



Let us note that in Mellin space, the fixed-order NLO cross sections for Drell-Yan 
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(3.91) 



with S'i(iV) = Y!j=i l/i- In the large-TV hmit, Si{N) ~ IniV + 1/(2 A^), and we get 



^i^-°-^(iV,M) 



a 
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1 + ^ (cf (2 ln2iV-2 lniVln^)+C'(;-) 
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,(3.92) 
(3.93) 



Comparing Eqs. (|3.88p . (|3.89p . (|3.92p . and (|3.93p . we see that the expansion of the 
resummed cross section at order as correctly reproduces the fixed-order result in the 
large- limit, including even terms that are suppressed by 1/A^. 



3.4 Joint resummation 
3.4.1 Formalism 



Let us go back to the QCD factorization theorem of Eq. ()3.16p . writing the 
unpolarized hadronic cross section as 
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dM2 dql 



^ / AXa j dXbfa/hi{Xa,HF) fb/h2iXb,fJ-F) 



AM^dqi, ^if fJ-R 



(3.94) 



as the convolution of the partonic cross section aab with the universal distribution 
functions fa,b/hab' have explicitly shown the lower integration limits depending 
on the quantity r = M'^/s^, which approaches the limiting value r = 1 when the 
process is close to the hadronic threshold. In Mellin A^-space, the hadronic cross 
section naturally factorizes 

,,1^, 2 = / ^T~^fa/hii]^ + 1, fJ'F)fb/h2i^ + 1> fJ-F)^abiN; OsiHR), —, — ), 

dM^dgf, ^ Jc2ni piF f^R 

(3.95) 



where the contour C in the complex A^— space will be defined later and the A^- 
moments of the various quantities are defined according to the Mellin transform 



F{N) = I dxx^~^F{x), 
Jo 



(3.96) 
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for X = Xa,b,z,T and F = fa/ha,b/hi^^^<^- The jointly resummed cross section in 
A^-space is usually written, at the NLL accuracy, as [75l [771 178] 



(iV) 



d^b 



47r 



X exp 



■•^^ C,/;,,(M,6,Ar;a,(/iij),/ii.) 
Cg/;,,(M,6,A^;a,(^/j),;UF), (3.97) 



where the parton densities and the sums over the different possible initial partonic 
states are included in the definition of the C-coefficients. The exponent is given by 



Er\NMo.,{^R),—) 
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Ac{as{n)) In + Bc{as{f^)) 
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(3.98) 



The function x organizes the logarithms of N and b in joint resummation and can 
be defined as 



x{N,b) = b + 



N 



l + rjb/N' 



(3.99) 



where b = b M e^^ /2, N = Ne^^ and r/ = 1/4. The only requirement to be satisfied 
when introducing a particular form of x is that the leading- and next-to-leading 
logarithms in TV and b are correctly reproduced in the limits TV — > oo or 5 — > oo, 
respectively. The choice of Eq. ()3.99p avoids the introduction of sizeable subleading 
terms into perturbative expansions of the resummed cross section at a given order 
in a<j which are not present in fixed-order calculations [77]. As for QT-resummation, 
the A-term resums the soft radiation, while the i?-term accounts for the difference 
between the eikonal approximation and the full partonic cross section in the threshold 
region, i.e. the flavour-conserving collinear contributions. In the large- limit, these 
coefficients are directly connected to the leading terms in the one-loop diagonal 
anomalous dimension, calculated in the MS factorization scheme (see Eq. I3.78P [192] 



7cc(iV,a,) = -A,{as) IniV - + 0(l/iV), 
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Performing the integration in Eq. (j3.98p . we obtain the form factor up to NLL, 
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with A = Po/TTasifJ-R) Inx, and the (^-functions are given by 
a'-c^ 2 a + In (1 - 2 a) 

2A + ln(l-2A) 
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(3.104) 



The C-coefhcients are chosen to correspond to transverse-momentum resumma- 
tion for b ^ oo, N being fixed, and are defined by 

^-^ MM 

Cd/H{M, b, N; asM,f^F) = Cd/a, {N; as{ — ))Ua,a{N; —, flF)fa/H{N+l, lip). 

a,ai ^ ^ 

(3.105) 

The C-functions contain contributions to the singular behaviour at vanishing qx 
(but not at threshold) of the fixed order cross section, while Ujk represents the 
evolution of the parton densities from scale hf to scale M/x and the //-function 
of Eq. (j3.97p absorbs hard virtual contributions. The C- and //-functions can be 
expressed perturbatively in powers of Og, 
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(3.106) 
(3.107) 



3.4.2 Reorganization of the resummed cross section 

In order to explicitly factorize the dependence on the parameter X: it is possible 
to organize the resummation of the logarithms in analogy to the case of transverse- 
momentum resummation [61^ I62j . The hadronic resummed cross section can be 
written as 



dM2 dq^ 
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y^T~(-ab 



. . M M \ ( ^ ,^ , , M , 1 , , 

cc\N;as{fiR), — , — exp <^ Gdlnx; as{fiR), — )>. (3.108) 



The function Tiab^cc does not depend on the parameter x, and contains all the terms 
that behave as constant in the limits 6 ^ oo or ^ oo, 
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Using Eqs. ()3.97p . ()3.105p and the QCD evolution operator Uab{N; fi, fio), fulfilling 
the evolution equation 

"^^"d lI'/T^' = ^ 7ac(a.(/.)) U^N- /X, ^o), (3.110) 



which has the solution 



Uab{N;ii,no) = exp 



(3.111) 



we get the first order coefficient 



(3.112) 

The x-dependence appearing in the C-coefficients and in the evolution operator of 
Eq. (|3.1U5|) is factorized and included in the exponent Qc, which is equal to the 
exponent Ec defined in Eq. ()3.2p . provided we make the replacement 

B,{as) ^ B,{as) = B{as) + 2/?(a,) +2j{as), (3.113) 

din as 

At NLL accuracy, Eq. (j3.103p remains almost unchanged, since only the coefficient 
gc of Eq. (|3.104p has to be slightly modified, through the replacement 

S«-^«(iV) = i3«+27«. (3.114) 

Although the first-order coefficients C^!^ and Hc^^ are in principle resummation- 



scheme dependent [61] . this dependence cancels in the perturbative expression of 
TCab-^cc [62]. In the numerical code we developed for slepton pair production, we 
chose to implement the Drell-Yan resummation scheme and take Hq{as) = 1, the 
C-coefhcients being the same as for gy-resummation (see Eq. p.53p ) 

3.4.3 Inverse transform and matching 

Once resummation has been achieved in N- and 6-space, inverse transforms 
have to be performed in order to get back to the physical spaces. Again, special 
attention has to be paid to the singularities in the resummed exponent, related 
to the divergent behaviour near x = exp[7r/(2/3oas)], i.e. the Landau pole of the 
running strong coupling, and near h = —2N and b = —4N, where X = and 
infinity respectively. The integration contours of the inverse transforms in Mellin 
and impact-parameter space must therefore avoid all of these poles. 

As for g^-resummation, the 6— integration is performed by deforming the integra- 
tion contour with a diversion into the complex 6-space [76| , defining two integration 
branches 

6 = (cos (/9 lb z sin (/9)t, with < t < oo, (3.116) 
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under the condition that the integrand decreases sufficiently rapidly for large 
values. The Bessel function Jq is replaced by two auxiliary functions hi^2{z-, v) related 
to the Hankel functions (see Eq. ()3.6ip ) 



hi(z,v) = / d^e 

vr 



iz sinO 



vr 



—ivn 

— iVTT 



h2{z,v) = -i / dee-'''''"'. (3.117) 



Their sum is always hi{z,v) + h2{z,v) = 2Jq{z), but they distinguish positive and 
negative phases of the 6-contour, being then associated with only one of the two 
branches. 

The inverse Mellin transform is performed following a contour inspired by the 
Minimal Prescription [189j and the Principal Value Resummation |190j . where one 
again defines two branches 

N = C + ze^'^ with < z < oo,tt xj) > ^. (3.118) 

The parameter C is chosen in such a way that all the singularities related to the 
A^- moments of the parton densities are to the left of the integration contour. It has 
to lie within the range < C < exp[7r/(2/3oas) — 7^;] in order to have convergent 
inverse transform integrals for any choice of (j) and (p. 

Finally, a matching procedure of the NLL resummed cross section to the NLO 
result has to be performed in order to keep the full information contained in the 
fixed-order calculation and to avoid possible double-counting of the logarithmically 
enhanced contributions. A correct matching is achieved through the formula 



dV . . d2cT(f'-0-) 



(i2^(rcs) |-j2^(cxp) 



, , / diV „^ /• 6d6 ^ , ^, 

(3.119) 



dM^d^^'^'^'^^^-dM^^^'^'"^) 

where d2a(F-0-) is the fixed-order perturbative result at a given order in a^, 
d^a^^"'^ is the resummed cross section at a given logarithmic accuracy, 
the truncation of the resummed cross section to the same perturbative order as 
^^2^(F.o.)^ and where we have removed the scale dependences for brevity. 

The expansion of the resummed results reads 

dfj^^'^'P) MM y-^ 

{N,b;as{fiR),—,—) = 2^fa/hAN + l,iJF)fb/h,{^ + ^,I^F) 



dM2 dql ' ' ' ' I^F —r 

M M 



X(Tir^)(iV,6;a,(/x«),— ,— ), (3.120) 
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where cr^i^^^ is obtained by perturbatively expanding the resummed component 



M M. 



aifP) (N, b; as M ,—,—) = ^ (M) { 6,J,b + ^ 



n=l 



vr 



''ab—*cc 



M M\ 



.. (3.121) 



The perturbative coefficient S^") is a polynomial of degree 2n in Inx and H*^") 
embodies the constant part of the resummed cross section in the limit of 6 ^ oo or 
— > oo. The first order coefficient fP-'^ is given by 



with 



^{i;2) 



ab^cc 

The integrals 



KLA^Mx) = KC2c-c^^'x + t':^2Anx, (3.122) 
2 A^p6ca5-cb and siScc = "2 (^^'^-^ca^cfe + 5ca7a + 55fe7£0(-3.123) 



/n(gT,A^) 



hdh 



.h{bqT)\n^[b + 



N 



l + 7?6/iv) 



(3.124) 



appearing when performing the inverse 6— transform, must be computed numerically. 



3.5 Comparison between g^, threshold and joint resum- 
mations 

In this section, we summarize the main differences between the three resum- 
mation formalisms described in the previous sections. For all resummations, the 
Sudakov form factor at the NLL accuracy can be written as 

g(iV,L)=5(^)(A)L + 5(')(A), (3.125) 

where we have removed all scale dependences for brevity. The logarithm L, A and 
the two ^-functions are given in Tab. 13.11 



Let us first note that M corresponds to the final state invariant-mass in the 
joint resummation formalism, but to the resummation scale in the gy-resummation 
formalism (see Eq. ()3.49p ). Such a scale is not needed in joint resummation, since 
the logarithm cannot be rescaled with a factor that is A^- and 6-independent. 
Finally, 6o is defined as for gy-resummation, 6o = 2e~'^^ . 

For all three resummations, the exponent is taken to be universal and process- 
independent, since the enhanced logarithmic contributions have the same dynamical 
origin, i.e. the soft-gluon emission from the initial state, which cannot depend on 
the hard-scattering process. All the process-dependence is factorized outside of the 
exponent, i.e. in the "H-coefficient for joint and transverse-momentum resummation 
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Joint 


Threshold 


L = ln(...) 




M b , iV 


iV 


A 








5«(A) 


^(1) A+ln(l-A) 
/3o A 


2^(1) A+ln (i-a) 
/3o A 


1^ [2A+(l-2A)ln(l-2A)] 




+(C'°S-^)(A+'°('-*)) 


^^[2A+ln(l-2A) + iln2(l-2A)] 
_^[2A+ln(l-2A)]-^lnf 
+ ^[2A + ln(l-2A)]ln^ 



Table 3.1: Comparison of the exponent in the (universal) transverse-momentum, 
joint, and threshold resummation formalisms. 



and in the C-coefHcient for the threshold resummation. In the conjugate Mellin- and 
6-spaces, the resummed partonic cross section can then be written as 

aifj)(iV,L) = HS,(iV)cxp{g(Ar,L) + Fr}, (3.126) 
aif(Ar,L) = a^^O)c^^^^^)^^^^g^N,L)]. (3.127) 

-F^^ is the non-perturbative form factor parameterizing the non-perturbative effects 
relevant for soft-gluon emission with very small qr- The expressions for all the 
coefficients can be found in the previous sections. 

As a side remark, let us note that the '+l'-term introduced in the logarithm 
of the transverse-momentum resummation in order to minimize the impact of the 
resummation in the small-6 region (where the perturbative theory is reliable) is not 
necessary for joint resummation, since for small b the corresponding logarithm is 
not tending towards zero. 
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Chapter 4 



Slepton-pair production at 
hadron colliders 

The LO cross section for the production of non-mixing slepton pairs (see Fig. 
I4.ip has been calculated in [311 [32l [33l HI]. We generalize these results, including 
the mixing effects relevant for third generation sleptons, and we calculate single- 
and double-spin asymmetries, for both neutral (7, Z^) and charged (VF^) currents 
[47] , making numerical predictions for the RHIC collider and for possible upgrades 
of the Tevatron [195j and the LHC [196J as well. 

Sleptons being among the lightest super symmetric particles in many SUSY- 
breaking scenarios [l6], they often decay into the corresponding SM partner and 
the lightest stable SUSY particle, a possible signal for slepton pair production at 
hadron colliders consisting thus in a highly energetic lepton pair and associated 
missing energy. The corresponding SM background is mainly due to WW and ti 
production [48\ I49|. An accurate calculation of the transverse- momentum spectrum 
allows us to use the Cambridge (s)transverse mass to measure the slepton mass 
|50j and spin [51] and to distinguish then the signal from the background. When 
studying the transverse-momentum distribution of a slepton pair produced with an 
invariant-mass M in a hadronic collision, it is appropriate to separate the large-gy 
and small-gr regions. In the large-f^T region (qx > M), the use of fixed-order pertur- 
bation theory is fully justified, since the perturbative series is controlled by a small 
expansion parameter, as(M^), but in the small-gr region, where the coefficients of 
the perturbative expansion in as(M^) are enhanced by powers of large logarithmic 
terms, ln(M^/g|,), results based on fixed-order calculations are completely spoiled. 
However, a precise gr-spectrum can be obtained by systematically resumming these 
logarithms to all orders in a^, using the universal formalism described in Sec. I3.2.4[ 

The NLO QCD corrections have been calculated in [39], and the full SUSY-QCD 
corrections have been added in [40J, considering however only massive non-mixing 
squarks and gluinos in the loops. We extend this last work by including the mixing 
effects relative to the squarks appearing in the loops [73]. We also consider the 
threshold-enhanced contributions, due to soft-gluon emission from the initial state, 
which arise when the initial partons have just enough energy to produce the slepton 
pair in the final state. In this case, the mismatch between virtual corrections 
and phase-space suppressed real-gluon emission leads to the appearance of large 
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Figure 4.1: Feynman diagram for slepton pair {V = 7, Z") and slepton-sneutrino 
associated {V = W^) production at leading order in perturbative QCD. 

logarithmic terms a"[ln^"~^(l — z)/{l — z)]+ at the n**^ order of perturbation 
theory. When s is close to M^, the large logarithms have to be resummed, which 
is achieved through the exponentiation of the soft-gluon radiation, within the 
collinear-improved formalism described in Sec. 13.31 

Since a complete understanding of the soft-gluon effects in differential distribu- 
tions requires a study of the relation between the recoil corrections at small qx and 
the threshold-enhanced contributions, we finally present a joint treatment of these, 
within the method described in Sec. 13.41 



4.1 LO unpolarized cross section and spin asymmetries 

4.1.1 Analytical results 

In the following, we define the square of the weak coupling constant g'^ = 
e^/sin^Ow in terms of the electromagnetic fine structure constant a = e^/(47r) 
and the squared sine of the electroweak mixing angle xw = sin^ dw ■ The coupling 
strengths of left- and right-handed (s)leptons to the neutral electroweak current given 
in Sec. 12.41 are 

{Lff,z,Rff'z} = {2Tf -2efXw)>^Sff, (4.1) 
i^kirz^^ijrz} = {Lii'z SjiSli,Rii'z Sj2Sl2*}, (4.2) 

{Lqq'W,Rqq'w} = {V^ Cw Vgg' , 0} , (4.3) 

{hnw^^ku,w} = {V2cwSll,0}. (4.4) 

We express (un)polarized cross sections for the electroweak 2 — > 2 scattering pro- 
cesses 

in terms of the conventional Mandelstam variables, 

s = {Pa+Pb)^ , t = {pa-pif , and u={pa-p2f (4.6) 
and the masses of the neutral and charged electroweak gauge bosons mz and mw- 
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The neutral current differential cross section for the production of mixing slepton 
pairs in collisions of quarks with definite helicities ha^b is given by |47j 



27rQ^ 



dt 



ut — nij rrij 



3s2 
eq ei 6ij Re 



+ 



+ 



Eg ef5ij{l - hahb) 

l-ha)il + hb) Lggz + {l + ha)il-hb) RggZ 



4 ( 1 — X VK ) ( 1 — "T-l /■§) 
{l-ha){l + hb)LlgZ + {l + ha){l-hb)RlgZ 



32x2,(l-xiy)2(l-m|/s)2 



(4.7) 



the three terms representing the squared photon-contribution, the photon-Z'' inter- 
ference and the squared Z'^-contribution, respectively, while the purely left-handed, 
charged current cross section is 



dal 



dt 



ut — m? m? 



[l-ha){l + hb) 



Lqq'W Lj.f^fW 



32x2 (l-m2 



■ (4.8) 



Averaging over initial helicities, 
daW : 



dafi + dati + d<TW 1 + d^W 



(4.9) 



we obtain the unpolarized partonic cross sections 



da 
dt 



27ra^ 
~3^ 



ut — rrir m? 



4 ef6ij + - 



Cq ei 5ij Re 



LqqZ + RqqZ 



4 (1 — a;H/) (1 — w-l/s) 



+ 



da^ 
dt 



32x2^(l-XH/)2(l-m|/s)2 



3s2 



ut — m? m? , 



Lqq'W Ll-i:,^W 



(4.10) 
(4.11) 



which agrees for mass-degenerate non-mixing sleptons with the neutral current 
result of Ref. [31j and with the charged current result of Ref. [33]. Note that for 
invariant final state masses close to the Z^- or H^-pole, the Z^- and VT-propagators 
must be modified to include the decay width of the corresponding electroweak boson. 

Prom Eqs. ()4.7p and ()4.8p . one can easily calculate the polarized cross sections 

(4.12) 
(4.13) 



dAaW 


d^a 


-d^i:u 


-da% + da%_. 








dAaW 


_ d^a 


+ dat. 


-da%-da%_. 
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which then read [47j 



W 

LL 



dA(TL 



27ra2 
L 



ut — rrij rriy 



eq ei 5ij Re 



(4.14) 

LqqZ — RqqZ 



^qqZ ^qqZ 



^2xl^{l-xw?{l-ml/s) 



dAai 



-a'. 



(4.15) 
(4.16) 



Let us note that the single-polarized cross section corresponds to a polarized initial 
quark, the corresponding quantity for a polarized initial antiquark being obtained 
from the following definition, 



dAa 



(') 



da 



1,1 



da 



(') 



da 



(') 

-1,-1 



(4.17) 



These expressions show that it will be interesting to study the dependence of the 
neutral current single-spin asymmetry 



Aj 



dAoL 
da 



(4.18) 



on the SUSY-breaking parameters, since it is the single quantity remaining sensitive 
to these. Furthermore the squared photon contribution, insensitive to them, is elim- 
inated. Finally, this scenario may also be easier to implement experimentally, e.g. at 
the Tevatron, since protons are much more easily polarized than antiprotons |195j . 
Let us note that our polarized results agree with those of Ref. [44j for non-mixing 
sleptons after integration over t, provided that we put parentheses around the 
terms of Eqs. (5) and (7) of Ref. [Hj, and if we replace the index g by e in the first 
occurrence of (a^ + ebq) in Eq. (7) of Ref. 



4.1.2 Unpolarized cross section 

For the masses and widths of the electroweak gauge bosons, we use here 
older values of mz = 91.1876 GeV, mw = 80.425 GeV, Tz = 2.4952 GeV, and 
Tw = 2.124 GeV |197j . which are slightly different from the current ones (see Sec. 
I2.3.3p . Since the sfermion mixing is proportional to the mass of the corresponding 
SM partner (see Eq. (j2.42p ). it is numerically only important for third generation 
sleptons. Consequently, the lightest slepton is the lighter stau mass eigenstate fi 
and we focus our numerical studies on its production. 

The f mass limits imposed by LEP depend strongly on the assumed SUSY- 
breaking mechanism, the mass difference between the f and the LSP, and the 
mixing angle 9r. The weakest limit of 52 GeV is found for GMSB models and stau 
decays to gravitinos, if no constraints on their mass difference are imposed [198j . 
This is the scenario that we will study for the RHIC collider, which has the most 
restricted hadronic centre-of-mass energy, equal to 500 GeV. For the Tevatron and 
the LHC, with their considerably larger centre-of-mass energies of 1.96 and 14 TeV, 
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Figure 4.2: Unpolarized hadronic cross sections for non-mixing f pair production of 
at the RHIC, Tevatron, and LHC colliders as a function of . The vertical lines 
indicate the two different f mass limits of 52 |198j and 81.9 GeV [199j. 

respectively, we will, however, impose the stricter mass limit of 81.9 GeV |199j . which 
is valid for staus decaying into neutralinos, with a mass difference of at least 15 GeV. 

Thanks to the QCD factorization theorem, the unpolarized hadronic cross section 

Cr = X]/ / (^Xbfa/haiXa,I^F)fb/htiXb,I^F)^ab (4.19) 

a,b -^^Z^- 

can be written as the convolution of the relevant partonic cross section dab (in- 
tegrated here over t) with the universal distribution functions fa^b/hat °^ partons 
a,b inside the hadrons /ia,f), which depend on the longitudinal momentum fractions 
of the two partons Xa^b and on the unphysical factorization scale fip- In order to 
employ a consistent set of unpolarized and polarized parton densities (see next 
subsection), we choose the LO set of GRV98 |200j for our unpolarized predictions 
at the factorization scale = ■ 

In Fig. 14. 2[ we show the unpolarized hadronic cross sections for non-mixing f 
pair production at the RHIC, Tevatron, and LHC colliders as a function of the f 
physical mass. The observation of tau sleptons will be difficult at RHIC, which is the 
only existing polarized hadron collider, but they will be detectable at the Tevatron, 
extending considerably the discovery reach beyond the current limits. In contrast, at 
the LHC, f pair production will be visible up to masses of 1 TeV. Before application 
of any experimental cuts, the SUSY signal cross section is at least three orders of 
magnitude smaller than the corresponding SM background coming from tau lepton 
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Figure 4.3: Longitudinal single-spin asymmetry ^4^ as a function of the cosine of the 
stau mixing angle for fi pair production at RHIC. 



pair production. Evaluated using the physical tau mass of rrir = 1.77699 GeV as 
factorization scale for the GRV98 LO parton densities, this cross section is equal 
to 1.7, 3.4, and 8.3 nb for the RHIC, Tevatron, and LHC colliders, respectively. 
Imposing an invariant-mass cut on the observed lepton pair and a minimal missing 
transverse energy will, however, greatly improve the signal-to-background ratio. In 
addition, as we will see in the next section, asymmetries may provide an important 
tool to further distinguish the SUSY signal from the corresponding SM process. 

4.1.3 Single-spin asymmetries 

Using again the QCD factorization theorem, we calculate the hadronic cross 
section for longitudinally polarized hadrons ha with unpolarized hadrons hb 

^^L = ^ [ dXa f dXbAfa/h,{Xa,fJ^F)fb/ht{^b,f^F)dAa^ab (4.20) 

through a convolution of polarized {A fa/A) and unpolarized (fb/s) parton densities 
with the singly polarized partonic cross section A&l (integrated over t). As 
mentioned above, we employ a consistent set of unpolarized |200j and polarized 
[20 Ij LO parton densities. The theoretical uncertainty due to the less well known 
polarized parton densities is estimated by showing our numerical predictions for 
both the GRSV2000 LO standard (STD) and valence (VAL) parameterizations. 

The physical masses of the SUSY particles and the mixing angles are computed 
with the program SuSpect 2.3 |202j . including a consistent calculation of the Higgs 
mass, with all one-loop and the dominant two-loop radiative corrections in the 
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Figure 4.4: Dependence of the longitudinal single-spin asymmetry Al on the cosine 
of the stau mixing angle for fi pair production at the Tevatron (top) and at the LHC 
(bottom). 
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renormalization group equations that link the restricted set of SUSY-breaking 
parameters at the gauge couphng unification scale to the complete set of observable 
SUSY masses and mixing angles at the electroweak scale. We choose three SPS 
benchmark points [IHlllB], the GMSB point SPS 7 with a light tau slepton decaying 
to a gravitino for RHIC and its very limited mass range, the typical mSUGRA point 
SPS la' with an intermediate value of tan (3 = 10 and a slightly reduced common 
scalar mass of rriQ = 70 GeV for the Tevatron, and the mSUGRA point SPS 4 with 
a large scalar mass of mo = 400 GeV and large tan /? = 50, which enhances mixing, 
for the LHC with its larger mass range. Since O-f depends directly on the universal 
soft SUSY-breaking mass scale A in the GMSB model and on the trilinear coupling 
Aq in the mSUGRA models, we test the sensitivity of the single-spin asymmetry on 
these parameters (see Figs. I4.3ll4.4p . 

In Fig. 14.31 we show the single-spin asymmetry for fi pair production at RHIC 
as a function of the cosine of the stau mixing angle. The asymmetry is quite large 
and depends strongly on the stau mixing angle. However, very large values of 
cos^f, corresponding to f masses below 52 GeV may already be excluded by LEP 
|198| . while small values of cos^f may be unaccessible to RHIC limited luminosity, 
which is not expected to exceed 1 fb~^. Polarization of the proton beam will also 
not be perfect, and the calculated asymmetries should be multiplied by the degree 
of beam polarization Pl ~ 0.7. The uncertainty introduced by the polarized parton 
densities increases considerably to the left of the plot, where the stau mass 41 GeV 
< rrir < 156 GeV and the associated values of the parton momentum fractions 
Xa^b — 2mf /y/sh become large. While the SM background cross section is then still 
two orders of magnitude larger than the SUSY signal after imposing an invariant- 
mass cut on the observed tau lepton pair, the SM asymmetry of -0.04 (-0.10) for STD 
(VAL) polarized parton densities can clearly be distinguished due to its different sign. 

An asymmetry measurement at an upgraded Tevatron would be extremely 
valuable, since the predicted asymmetry is very sizeable in the entire viable SUSY 
parameter range, and depends strongly on the parameter Aq, and then the stau 
mixing angle (top figure of Fig. 14. 4p . Unfortunately, the parton density uncertainties 
are still large, but will be reduced considerably in the future through more precise 
measurements at the COMPASS, HERMES, PHENIX, and STAR experiments. As 
a recent experimental study demonstrates, events with tau leptons with associated 
missing energy larger than 20 GeV can be identified with the CDF-II detector, 
considering hadronic tau decays |203j . Again, the negative SM asymmetry of -0.09 
(for both polarized parton densities) would be clearly distinguishable due to its 
opposite sign. 

For the LHC, where studies of tau slepton identification with the ATLAS detector 
|204] and tau tagging with the CMS detector j205j have recently been performed, 
SUSY masses should be observable up to the TeV-range. We show the predicted 
asymmetry for a possible polarization upgrade of this collider in the bottom panel of 
Fig. 14.41 It is still comfortably large and has again the opposite sign with respect to 
the SM asymmetry of -0.02 (for both polarized parton densities). The dependence 
of the asymmetry on the stau mixing angle is, however, also reduced, while the 
uncertainties from the polarized parton densities, which are not yet well known 
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at the small x values relevant for the large LHC centre-of-mass energy, are quite 
enhanced. 



4.2 Transverse-momentum spectrum 

4.2.1 Analytical results 

As said in the introduction, the main SM background to slepton pair produc- 
tion comes from WW and tt production. However, an accurate calculation of the 
transverse-momentum spectrum would allow us to distinguish the signal from the 
background. A transverse-momentum distribution for slepton pair production at 
leading order is induced by a final state containing one slepton pair and one QCD 
jet. The relative partonic cross section can be divided into a resummed and a finite 
part, as in Eq. ()3.64p . 

dMHq^ dMMql dMMq^' ^' ' 

where a, b label the partons which take part in the hard process. In Mellin A^-space, 
the (partonic) resummed component is deduced from Eqs. ()3.47p and ()3.49p . 

-^j0^{qT, M, s;asinR), fiR, fip) = db - Jo(bqT)W[t,{N,b; M, fip, fJ,R), 

(4.22) 

W[^{N,b■,M,^lF,^^R)=n^b(N, asM;—,—,^)exp\g{N,L,asM;—,^)}, 

(4.23) 

where at the NLL accuracy, the functions Ti and Q of Eqs. (|3.50p and (j3.54p read 



p M M M ^ (LO)F.„. 



TT / ""^'"'^ [IR HF Q 



(4.24) 



GiN, L; as, — ,^) = L g('^ {- Po a^im) L) + g^'^ (N, - (3o c^sM U ^,^),(4.25) 

flR Q TT TT ^IR Q 

with 

L = ln(^ + l), (4.26) 

this definition of L allowing us to avoid the unwanted resummed contributions in 
the small-6 region. The explicit expression of g^^^ and g^"^^ in terms of the univer- 
sal perturbative coefficients A^q \ A''q \ and sj^^ (see Eqs. ([325]), (l3:26H and[32Zl)) 
are given by Eqs. (|3.55p and (j3.56p . Slepton pair production being similar to SM 
lepton pair production, except for the final state, we can use the Drell-Yan coeffi- 
cients from Eqs. (|33n) and (I332D . but the Born cross section a m 
Eq. (|4.24p corresponds obviously to the slepton pair production one, which can be 
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obtained for neutral and charged currents by integrating Eqs. ()4.10p and ()4.1ip over 
t, respectively, 



(LO) 



(M) 



9M2 



4 ef d^j + 



eq ei Sij{LqqZ + Rqqz) Re 



+ 



qqZ) 



4xw {1 - xw) {I - ml/M^) 

2 



32x2^(l-XH/)2(l-m|/M2)2 



9M2 



(4.27) 
(4.28) 



_32x2^(l-x^)2(l-m2,/M2)2_ 
where the slepton-mass dependence is factorized in the velocity 

P = yjl + mf/M^ + m^/M^ - 2{mf/M'^ + m2/M2 + m2/M4), (4.29) 
mj and mj being the two masses of the final state particles. 

The second term {da'^^'^ /dM'^dq'j^) in Eq. ()4.2ip is free of divergent contributions 
and can be computed by fixed-order truncation of the perturbative series, as in 
Eq. (j3.62p . In order to be consistently matched with the resummed contribution at 
intermediate qt {qt — M), this term is evaluated starting from the usual perturbative 
calculation of the partonic cross section from which we subtract the expansion of the 
resummed component of the cross section at the same perturbative order. In order 
to consistently perform a NLL+LO matching, we need the fixed-order cross section 
relative to the production of a slepton pair with non- vanishing transverse-momentum 
at order as, i.e. one slepton pair plus one parton \39\ I40j . The different channels 
read 



daqg 




dM2d4 


2tts 


dagg 




dM^dq^ 


2tts 






dM2dg2 


2tts 



Agg{s,t,u;M')a'^i:,>{M), 



(4.30) 
(4.31) 
(4.32) 



where s, t, and u denote the partonic Mandelstam variables relative to the 2 — > 2 
subprocess ab ^ ^, Z^, + one parton. 



S — XaXhSh, 

t = M2-^S;,(M2+4)xfee^ 

u = -Jsh{M^ + q^)xae-y, 



(4.33) 
(4.34) 

(4.35) 



y being the rapidity of the slepton pair. We have defined the ^at-functions as [206], 

Agg{s,t,u;M^) = -I - + - + —— \ (4.36) 
Agg{s,t,u;M^) = -Agg{u,t,s;M'^). (4.37) 
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It is known [57] that the transverse-momentum distribution is affected by 
non-perturbative effects which become important in the iarge-6 region. In the 
case of electroweak boson production, these contributions are usually parameter- 
ized by multiplying the function in Eq. ()4.23p by a NP form factor F'^^{b) 
[n2\ 11761 \m\ \nS\ I179j . whose coefficients are obtained through global fits to 
DY data. We include in our analysis three different parameterizations of non- 
perturbative effects corresponding to three different choices of the form factor: the 
Ladinsky-Yuan (LY-G) [176J, Brock-Landry-Nadolsky-Yuan (BLNY) [178j, and the 
recent Konychev-Nadolsky (KN) jl79| form factor. The explicit expressions can be 
found in Eqs. ([SSSD, <^^^ and dOi]) . 

4.2.2 Numerical results 

We present quantitative results for the ^T-spectrum of slepton pair (slepton- 
sneutrino associated) production at NLL-I-LO accuracy at the LHC collider. For the 
masses and widths of the electroweak gauge bosons, the electroweak mixing angle 
and the electromagnetic fine structure constant, we use the same values as those 
given in Sec. I4.1.2^ while the CKM-matrix elements are computed as in Sec. 12.3.31 
As in the previous section, we focus our study on the lightest slepton mass eigenstate 
fi and thus we consider the processes 



We use the MRST (2004) NLO set of parton distribution functions [207] and 
Us is evaluated at two-loop accuracy. We fix the resummation scale Q (see Eq. 
(|3.49p ) equal to the invariant-mass M of the slepton (slepton-sneutrino) pair and 
we allow ^ = = to vary between M/2 and 2M to estimate the perturbative 
uncertainty. We also integrate Eq. ()4.2ip with respect to M^, taking as lower limit 
the energy threshold for fif* {tiu-t) production and as upper limit the hadronic 
energy = 14 TeV at the LHC). 

In the case of fif* production (neutral current process, see the top panel of 
Fig. 14. Sp . we choose the SPS 7 GMSB benchmark point [l5] which gives, after the 
renormalization group evolution of the SUSY-breaking parameters performed by 
the SuSpect computer program [202j . a light fi of mass nif^ = 114 GeV. In the case 
of fiu* + flvr production (charged current process, see the bottom panel of Fig. 
14. 5p . we use instead the SPS la mSUGRA benchmark point which gives a light fi 
of mass nif^ = 136 GeV as well as a light Uj- of mass m/^^ = 196 GeV. 

In both cases we plot the LO result (dashed line) and the total NLL-I-LO 
matched result (solid line) with their uncertainty band relative to scale variation 
(yellow band for LO result, and green band for NLL-I-LO result). The expansion 
of the resummation formula at LO (dotted line) is also shown. We can see that 
the LO result diverges, as expected from the generic form of the fixed-order cross 
section in Eq. (|3.18p . for both processes as ^ 0. The asymptotic expansion of 
the resummation formula at LO is in very good agreement with LO both at small 



qq' 



qq -^-i,Z^ ^ fif*. 



(4.38) 
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pp -^y, Z^tjXj* at the LHC 




Figure 4.5: Differential cross section for the processes pp — > fifj^ (top) and pp 
fiv* + flvr (bottom). NLL+LO matched result, LO result, asymptotic expansion 
of the resummation formula and A-parameter are shown. 
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and intermediate values of qt-, from which we can conclude that the cross section 
is clearly dominated by the logarithms that we are resumming in this kinematical 
region. The effect of resummation is clearly visible at small and intermediate values 
of qT, the resummation-improved result being nearly 39% (36%) higher at = 50 
GeV than the pure fixed-order result in the neutral (charged) current case. When 
integrated over qt, the former leads to a total cross section of 66.8 fb (12.9 fb) in 
good agreement (within 3.5%) with the QCD-corrected total cross section at 0{as) 

m- 

The scale dependence is clearly improved in both cases with respect to the pure 
fixed-order calculations. In the small and intermediate gT-region (up to 100 GeV) 
the effect of scale variation is 10% for the LO result, while it is always less than 5% 
for the NLL-I-LO curve. 



Finally, we study the dependence of the total NLL-I-LO matched result on non- 
perturbative effects. We show the quantity A 

A = ^ ^ — (4.39) 

as a function of the transverse momentum of the slepton pair. The parameter A gives 
thus an estimate of the contributions from the different NP parameterizations (LY- 
G, BLNY, KN) that we included in the resummed formula, which are under good 
control, since they are always less than 5% for qx > 5 GeV and thus considerably 
smaller than the resummation effects. 



4.3 Invariant-mass distributions 

4.3.1 Next-to-leading order calculations 

Restoring the various scale dependences in the QCD factorization theorem of Eq. 
(|4.19p . the invariant-mass distribution of a given process can be written as 



0- = X]/ "^^"^ / (^Xbfa/haiXa,fJ-F)fb/htixb,fJ'F)^abiz,M;asifJ.R), 
a,b •^^Z^" ^ 



M M\ 



(4.40) 



where M is the invariant-mass. The partonic cross section is usually expanded in 
powers of ag 

Concerning slepton pair hadroproduction, we have computed the LO (n = 0) coef- 
ficients in the previous sections, which read, showing explicitly the dependence on 
the scaling variable z = M'^ /s, 

-Sl)^^;-'-) = '^'^^^{l-z), (4.42) 
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ctq and (Tq being defined in Eqs. ()4.27p and ()4.28p for the neutral and charged currents, 
respectively, 



MM) 



+ 



9M2 



el ef 6ij + 



Cq ei 6ij{LqqZ + Rqqz) Re 



Axw {I - xw) {I - ml/M^) 

2 



?,2xl,{l-xw?{l-ml/M^) 
9M2 



32x2^(l-XH/)2(l-myM2)2 



(4.43) 



(4.44) 



The NLO QCD and SUSY-QCD corrections, i.e. the n = 1-coefficients in the 
perturbative expansion of Eq. (j4.41|) . have been studied for non- mixing sleptons in 
[39\ I40j . At NLO, the quark-antiquark annihilation process receives contributions 
from virtual gluon exchange (see upper part of Fig. 14. 6|) . real gluon emission (see 
Fig. 14. 7p diagrams, and we also have to take into account the quark-gluon initiated 
subprocesses (see Fig. 14. 8p . 



The infrared and collinear singularities of the three-parton cross sections are 
extracted using the dipole subtraction formalism |208j . and the virtual corrections 
have been evaluated in the MS renormalization scheme. For the SM QCD diagrams 
one has the well-known results |39l I209j 
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(4.45) 

3z 



(4.46) 



which expose the LO cross-sections (Tg^(M) in factorized form. Cf = 4/3 and Tr = 
1/2 are the usual QCD colour factors, and Pqq^qg are the Altarelli-Parisi splitting 
functions [21Uj 



Cf 
2 
Tr 



■6il-z) + 



(l-z) 



(1 + ^) 



and 



[z' + (1 - zf] . 



(4.47) 
(4.48) 



We remind the reader that our normalization corresponds to a perturbative 
expansion in powers of Os/vr. 
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Figure 4.6: Contributions of virtual diagrams for slepton pair {V = 7, Z^) and 
slepton-sneutrino associated {V = W^) production at next-to-leading order in per- 
turbative QCD. The first and second lines show the QCD and SUSY-QCD cor- 
rections, respectively. In the SUSY-QCD case, one has to sum over squark mass- 
eigenstates k,l = 1,2. 




Figure 4.7: Contributions from real gluon emission diagrams for slepton pair (V = 
7, Z^) and slepton-sneutrino associated {V = W^) production at next-to-leading 
order in perturbative QCD. 




Figure 4.8: Contributions from qg diagrams for slepton pair {V = 7, Z^) and slepton- 
sneutrino associated {V = W^) production at next-to- leading order in perturbative 
QCD. 



73 



The three lower diagrams of Fig. 14.61 contain SUSY-QCD corrections. We gener- 
ahze the results from [40j to the case of mixed squark mass eigenstates k,l = 1,2 in 
the virtual loop diagrams, obtaining for neutral and charged currents [73j 
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(4.49) 
.(4.50) 



where now only the diagonal squared photon contribution to the Born cross section 
factorizes. The virtual loop coefficients f^z-, fz and fw are 



i=l,2 



/7 = 2+5: 

+ 

+ 2 



2m?-2m?+M2 



M2 



m| + (M^ — 2 m? ) m? + 
M2 



Co/ (0,M2,0,m|,m?,m|) 



(4.51) 



/7Z = 2 (L^gg + RZ. 
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(4.54) 



The functions -Bo/(p )"^i)"^2)' Bqj:{p^ ,1711,7712) and CQf{pf,{pi + 
P2)'^ ,P2,rn'l, 1712,771^) are the finite parts of the scalar two- and three-point 
functions 
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Our results agree with those of Ref. [21 Ij in the case of mass-degenerate non-mixing 
squarks. Note that the quark mass, which appears in the off-diagonal mass matrix 
elements of the squarks running in the loops (see Eq. ()2.36p ). corresponds to a linear 
Yukawa coupling in the superpotential and can not be neglected, even if it is much 
smaller than the total centre-of-mass energy of the colliding partons allowing for a 
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massless factorization inside the outer hadrons. The full NLO contributions to the 
cross section are then given by 



(1) 

■(') 



z,M;—,— )-- 



(1;QCD) 



(') 



M M 



(7, 



M M\ 



(1;SUSY) 



■(') 



z,M: 



M M\ 



4.3.2 Threshold-enhanced contributions 



) ,(4.58) 
(4.59) 



Eqs. ()4.45p and ()4.46p explicitly show logarithmic terms of the form as[ln(l — 
z)/(l — z)]+. When the initial partons have just enough energy to produce the slepton 
pair in the final state, i.e. z is close to one, these terms can become large and have 
to be resummed to all order in in order to correctly quantify the effect of this 
set of corrections, corresponding to soft-collinear gluon emission. Collinear parton- 
emission can also be taken into account, using the collinear-improved resummation 
formalism which is described in Sec. 13.31 for Drell-Yan pair production, which can 
directly be applied to slepton pair production by replacing the Born cross sections 
by those of Eqs. (|4.43|) and (j4.44|) . In Mellin space, we have then 
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(4.62) 



the Sudakov form factor S being defined by Eqs. (|3.79p . (|3.80p and (j3.8ip . The errors 
due to missing high-order corrections can be estimated by comparing our results with 
those obtained by taking the resummed contribution given by Eq. (j3.85p . where the 
hard contributions exponentiate as well. 
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The matching with the fixed-order results is achieved through Eq. (|3.87p 
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The truncation of the resummed cross section to the same perturbative order as 
^(F.o.)^ reads for a NLL-I-NLO matching (i.e. at order as) as in Eqs. (j3.88p and 
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Let us note that in Mellin-space, the large- hmit of the fixed order NLO cross 
sections of Eqs. (|4.45p and (j4.46p is correctly reproduced by the expansion of the 
resummed cross section at order as, including even terms that are suppressed by 
as expected from the collinear-improved resummation formalism (see Eqs. ()3.92p and 

4.3.3 Numerical results 

We take current masses and widths of the electroweak gauge bosons, electroweak 
mixing angle, electromagnetic fine structure constant and CKM matrix elements, 
as in Sec. 12.3.31 [149J, and the physical masses of the SUSY particles and the mixing 
angles are obtained with the computer program SuSpect [20^ . We choose one 
mSUGRA point, SPS la, and one GMSB point, SPS 7, with their associated model 
lines, as benchmarks for our numerical study ^5]. SPS la is a typical mSUGRA 
point with an intermediate value of tan (3 = W and ;U > 0, and the model line 
attached to it is specified by mo = — = 0-4 '^-1/2- Foi' "^-1/2 = 250 GeV, this 
SUSY-breaking scenario leads to light sleptons fi, ei, f2, 62, i^r and with masses 
of 136.2, 146.4, 216.3, 212.3, 196.1 and 197.1 GeV and to heavy squarks with 
masses around 500-600 GeV. However, even if the top-squark mass eigenstate ti is 
slightly lighter, it does nonetheless not contribute to the virtual squark loops due 
to the negligible top-quark density in the proton (we consider flavour-conserving 
SUSY-loops in Fig. 14. 6p . SPS 7 is a GMSB scenario with a fi as the next-to-lightest 
SUSY particle and an effective SUSY-breaking scale A = 40 TeV, A^mes = 3 
messenger fields of mass M^es = 80 TeV, tan /3 = 15, and /x > 0, which leads 
again to light sleptons with masses of 114.8, 121.1, 263.9, 262.1, 249.5 and 249.9 
GeV, respectively, and even heavier squarks with masses around 800-900 GeV. Its 
model line is defined by Mmes = 2A. The slepton masses and mixing angles are 
actually quite similar for the SPS la mSUGRA and SPS 7 GMSB points, so that 
the corresponding production cross sections will not differ significantly. Slepton 
detection will, however, be slightly different in both scenarios, as the sleptons decay 
to a relatively massive neutralino LSP at SPS la, but to a very light gravitino 
LSP at SPS 7. The lightest tau slepton thus decays into a tau lepton and missing 
transverse energy. 

Our cross sections are calculated for the Tevatron pp-collider, as well as for the 
LHC pp-collider. For the LO (NLO and NLL) predictions, we use the LO 2001 
[212] (NLO 2004 [207]) MRST-sets of parton distribution functions. For the NLO 
and NLL predictions, Os is evaluated with the corresponding value of ^ = 255 
MeV at two- loop accuracy. We fix the unphysical scales fip and /j^r equal to the 
invariant-mass M of the slepton (slepton-sneutrino) pair. 

The invariant-mass distribution M^da/dM for first- (and equal-mass second-) 
generation sleptons at the Tevatron is shown in the top panel of Fig. 14. 9[ and the 
one for (slightly lighter) third-generation sleptons at the LHC in the bottom panel 
of the same figure. In both cases, we have chosen the SPS 7 GMSB benchmark 
point. The differential cross section da/dM has been multiplied by a factor in 
order to remove the leading mass dependence of propagator and phase space factors. 
As is to be expected for P-wave production of scalar particles, the distributions 
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Figure 4.9: Invariant-mass distribution da/dM of cl pairs at the Tevatron (top) 
and fi pairs at the LHC (bottom) for the benchmark point SPS 7. We show the 
total NLL+NLO matched and the fixed-order NLO SUSY-QCD and LO QCD re- 
sults, including the respective scale uncertainties as horizontally hatched, vertically 
hatched and shaded bands. 
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rise above the threshold at ^/s = 2mi with the third power of the slepton velocity 
(3, see Eq. (|4.27p . and peak at about 100 GeV above threshold (at 370 GeV for 
M^da/dM and 310 GeV for dcr/dM for the Tevatron; 410 GeV and 300 GeV for the 
LHC), before falling off steeply due to the s-channel propagator and the decreasing 
parton luminosity. Furthermore, it can also be seen that the QCD corrections do 
not alter the P-wave velocity dependence close to threshold. At the Tevatron, the 
total and NLO SUSY-QCD predictions exceed the maximal LO cross section by 36 
and 31%, respectively, whereas at the LHC, the maximal cross section increases 
by 28 and 27%. Threshold resummation effects are thus clearly more important at 
the Tevatron, where the hadronic centre-of-mass energy is limited and the scaling 
variable r = M^/s/j is closer to one, and they increase with M to the right of both 
plots. 

The maximal theoretical error is estimated in Fig. 14.91 by an independent 
variation of the factorization and renormalization scales between M/2 and 2M. 
It is indicated as a shaded, vertically, and horizontally hatched band for the LO, 
NLO SUSY-QCD, and the total prediction. At LO, the only dependence comes 
from the factorization scale. It increases with the momentum-fraction x of the 
partons in the proton or anti-proton and is therefore already substantial for small 
M at the Tevatron, but only for larger M at the LHC. At NLO, this dependence is 
reduced due to the factorization of initial-state singularities, but a strong additional 
dependence is introduced by the renormalization scale in the coupling as{fJ-R). 
After resummation, this dependence is reduced as well, so that the total scale 
uncertainty at the Tevatron diminishes from 20%-35% for NLO to only 16%-17% 
for the matched resummed result. The reduction is, of course, more important in 
the large-M region. At the LHC, where ag is evaluated at a larger renormalization 
scale and is thus less sensitive to it, the corresponding numbers are 18%-25% and 
15%~17%. 



For the mSUGRA scenario SPS la, we show in Figs. 14.101 and 14.111 the cross 
section correction factors 

daLO/dM' ^ ' 

where i labels the corrections induced by NLO QCD (Eqs. ()4.45p and ()4.46p '). 
additional NLO SUSY-QCD (Eqs. (|i38]) and g39D), resummation (Eq. (fOT]) ). and 
the matched total contributions (Eq. ()4.64p ) as well as the fixed-order expansion 
(Eqs. ()4.65p and (j4.66p ) of the resummation contribution as a function of the 
invariant-mass M. As one can see immediately, the mass-dependence of these 
corrections for charged-current associated production of sleptons and sneutrinos 
(lower parts of Figs. 14.101 and Fig. 14. lip does not differ substantially from the 
mass-dependence of the neutral-current production of slepton pairs (upper parts). 

At the Tevatron (Fig. I4.10p . where we are close to the threshold, resummation 
effects are already important at low M (4%) and increase to sizeable 16% at large 
M. The NLO QCD result is thus dominated by large logarithms and coincides with 
the expanded result at the permille level. In addition, the relative importance of 
the (finite) SUSY-QCD contributions is reduced, and the total prediction coincides 



79 



l.4r 



1.35 — 



1.3 — 



p p Z e,e* at the Tevatron 




1.25 — 



1.2 - 



SPS1a, irij =146.4 GeV 



Total = Resummation 

NLO SUSY-QCD 



300 



NLOQCD = Expansion 



400 



500 



I \ \ I \ \ \ I I L 

600 700 
M- . [GeV] 



800 



900 



1.25 — 



1.2 



p p ^ W ^ V^* + c. c. at the Tevatron 




Total = Resummation 
NLO SUSY-QCD 
NLOQCD = Expansion 



500 



600 700 
M. . [GeV] 



Figure 4.10: K-factors as defined in Eq. (j4.67p for pair (top) and associated ei^l 
production (bottom) at the Tevatron for the benchmark point SPS la. We show the 
total NLL+NLO matched result, which is almost identical to the purely resummed 
result at NLL, as well as the fixed-order NLO SUSY-QCD and QCD results. The 
latter practically coincides with the resummed result expanded up to NLO. 
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Figure 4.11: i^- factors as defined in Eq. (j4.67p for fi pair (top) and associated fiv* 
production (bottom) at the LHC for the benchmark point SPS la. We show the 
total NLL+NLO matched result, the resummed result at NLL, the fixed-order NLO 
SUSY-QCD and QCD results, and the resummed result expanded up to NLO. 
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with the resummed prediction, since fixed-order and expanded contributions can- 
cel each other in Eq. (j4.64p . We have also verified that exponentiating the finite 
(A^-independent) terms collected in the coefficient function C^^|,)j as proposed in 
|193j (see Eq. I4.63p . leads only to a 0.6%-0.8% increase of the matched resummed 
result. The Tevatron being a pp-collider, the total cross section is dominated by 
(^g-annihilation, and ggr-scattering contributes at most 1% at small M (or small x), 
where the gluon density is still appreciable. Integration over M leads to total cross 
sections for the neutral (charged) current processes in Fig. 14.101 of 4.12 (3.92) fb in 
LO, 5.3 (4.96) fb in NLO (SUSY-)QCD, and 5.55 (5.28) fb for the matched resummed 
calculation. The corresponding (global) X-factors 

= ^ = JdMdg- /dM 

are then 1.29 (1-27) at fixed-order and 1.35 (1.35) with resummation. 

At the LHC, sleptons can be produced with relatively small invariant-mass M 
compared to the total available centre-of-mass energy, so that z = T/{xaXb) = 
M"^ / s <C 1 and the resummation of (1 — 2;)-logarithms is less important. This is 
particularly true for the production of the light mass-eigenstates of mixing third- 
generation sleptons, as shown in Fig. 14. Hi In the low-M (left) parts of these plots, 
the total result is less than 0.5% larger than the NLO (SUSY-)QCD result. Only at 
large M the logarithms become important and lead to a 7% increase of the i^-factor 
with resummation over the fixed-order result. In this region, the resummed result ap- 
proaches the total prediction, since the NLO QCD calculation is dominated by large 
logarithms and approaches the expanded resummed result. However, we are still far 
from the hadronic threshold region, and a consistent matching of both resummed 
and fixed-order contributions is needed. At low M, where finite terms dominate, 
the resummed contribution is close to its fixed-order expansion and disappears with 
M. In the intermediate-M region, one can observe the effect of SUSY-QCD con- 
tributions, in particular the one coming from the qqg-vertex correction (lower left 
diagram in Fig. 14. 6|) . As M > 2mq, one crosses the threshold for squark pair produc- 
tion and observes a resonance in Fig. 14.111 As for the Tevatron, exponentiating the 
finite (A^-independent) terms collected in the coefficient function C^i|,) leads only 
to a 1% increase of the matched resummed result. The LHC being a high-energy 
pp-collider, it has a significant gluon- luminosity, in particular at small M (or x), 
and indeed the gf^-subprocess changes (lowers) the total cross section by 7% at small 
M and 3% at large M. After integration over M, we obtain total cross sections of 
27 (9.59) fb in LO, 34.3 (12.3) fb in NLO SUSY-QCD, and 34.6 (12.5) fb for the 
resummed-improved result, corresponding to global X-factors of 1.28 for fixed-order 
and 1.29 for the matched resummed cross section for both processes. Resummation 
of large logarithms is thus not as important as for the Tevatron at the benchmark 
point SPS la. 

4.4 Total cross section at NLO 

In this section, we study the dependence of total, i.e. invariant-mass integrated, 
slepton pair and slepton-sneutrino cross sections on three different scales: first the 
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Figure 4.12: Total cross sections (top) and i^'-factors as defined in Eq. ()4.68p (bot- 
tom) for first- (and second-) generation slepton pair and slepton-sneutrino associated 
production at the Tevatron along the model line attached to the SPS 7 benchmark 
point (vertical dashed line). We show the total NLL-I-NLO matched and the fixed- 
order NLO (SUSY-)QCD and LO QCD results. 
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Figure 4.13: Total cross sections (top) and i^'- factors as defined in Eq. ()4.68p (bot- 
tom) for third-generation slepton pair and slepton-sneutrino associated production 
at the LHC along the model line attached to the SPS 7 benchmark point (verti- 
cal dashed line). We show the total NLL+NLO matched and the fixed-order NLO 
(SUSY-)QCD and LO QCD results. 
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Figure 4.14: Dependence of the total cross section for first- (and second-) generation 
slepton pairs at the Tevatron on the factorization scale (top), renormalization scale 
(middle), and both scales (bottom) for the SPS 7 benchmark point. 
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Figure 4.15: Dependence of the total cross section for third-generation slepton pairs 
at the LHC on the factorization scale (top), renormalization scale (middle), and both 
scales (bottom) for the SPS 7 benchmark point. 
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dependence on the effective SUSY-breaking scale A as defined in tlie GMSB model 
line SPS 7, and then the dependence on the renormalization and factorization scales 
Hr^f at the point SPS 7. The benchmark point at A = 40 TeV will be indicated in 
the figures where A varies as a vertical dashed line, and we show in addition to the 
scale A the mass scale of the produced charged slepton, ranging from 80 (87.5) to 
280 (385) GeV for ei (h). 

In Fig. 14.121 total cross sections (top) and ii'-factors (bottom) are shown at the 
Tevatron, which is expected to produce a total integrated luminosity of 4-8 fb~^. 
Even for bl pair production the mass-range is limited to masses below 280 GeV, 
where the cross section reaches 0.1 fb and at most one event would be produced. 
NLO and resummation corrections are clearly important, as they increase the LO 
prediction by 18 to 28% (lower part of Fig. I4.12p . At the SPS 7 benchmark point, 
the corrections would thus induce a shift in the selectron mass as deduced from a 
total cross section measurement by about 8 GeV (cf. the two dashed lines in the 
upper part of Fig. I4.l2|) . By comparing the NLO and total predictions, one observes 
an increased importance of threshold resummation for heavier sleptons, as expected. 

We show in Fig. 14.131 the total cross sections in the restricted range A < 135 
TeV. The NLO and resummed corrections are again large (25-30%), but the 
resummation corrections only become appreciable for large SUSY-breaking scales 
(or slepton masses). The largest cross section is obtained for pair production of 
the light stau mass eigenstate, even though it has a large right-handed component. 
Conversely, the heavier stau mass eigenstate has a large left-handed component, 
so that its cross section is less suppressed. At the SPS 7 benchmark point, the 
corrections would again induce a shift in the slepton mass as deduced from a total 
cross section measurement by about 8 GeV, as can be seen for fi on the upper part 
of Fig. 14.131 cf. the two dashed lines. 

Finally, we consider the theoretical uncertainty of invariant-mass integrated to- 
tal cross sections at the Tevatron (Fig. I4.14p and the LHC (Fig. I4.15|) as induced by 
variations of the factorization scale (top), renormalization scale (middle), or both 
(bottom). The /iij-dependence (middle), which is absent in LO, is first introduced 
in NLO, but then tamed by the resummation procedure. On the other hand, the 
logarithmic /iir-dependence (top), already present through the PDFs at LO, is over- 
compensated (reduced) at NLO for the LHC (Tevatron) and then (further) stabilized 
by resummation. This works considerably better at the LHC, where at least one 
quark PDF is sea- like and the PDFs are evaluated at lower x, than at the Tevatron, 
where both PDFs can be valence-like and are evaluated at relatively large x. In 
total, the theoretical uncertainty at the Tevatron (LHC), defined by the ratio of the 
cross section difference at fip = fiji = 171^/2 and fJ-F = fJ-R = '^itT'I over their sum, 
increases from 20 (7) % in LO to 29 (17) % in NLO, but is then reduced again to 23 
(8) % for the resummed-improved prediction. 

4.5 Jointly resummed results 

We present here a joint treatment of the recoil corrections at small qt and the 
threshold-enhanced contributions near partonic threshold, allowing a complete un- 
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derstanding of the soft-gluon effects in differential distributions for slepton pair pro- 
duction at hadron colliders ^80j. We use the joint formalism described in Sec. 13.4 
and compare it to the qt- and threshold resummation formalisms of Sees. 13.2.41 and 
13.3.21 The fixed-order perturbative and the qt- and threshold-resummed results for 
the transverse-momentum and invariant-mass distributions are those presented in 
the previous subsections, while the jointly-resummed matched cross section (see Eq. 

dsns])) is 
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The resummed contribution is given by Eq. ()3.108p 
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while the expansion of the resummed component is (see Eq. ()3.12ip ) 
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dM2 dg2 
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vr 
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M M\ 



(4.71) 



The 0{as) coefficients of the perturbative functions 7i, Q and S, can be found in Sec. 
13. 4i We remind the reader only of the definition of the argument in the logarithm 
(see Eq. (iMD ). 



x(iV,6) = 6 + 



N 



l + rjb/N 



(4.72) 



For the masses and widths of the electroweak gauge bosons, the electroweak 
mixing angle and the electromagnetic fine structure constant, we use the same values 
as in the previous section, i.e. those given by the latest version of the PDG review 
|149j . We focus our study on the production of a right-handed selectron pair at the 
LHC, 



qq^j,Z^^eRe%. (4.73) 

We use the MRST (2004) NLO set of parton distribution functions [207], as is 
evaluated at two-loop accuracy, and we allow fip and fiR to vary between M/2 and 
2M to estimate the perturbative uncertainty. We choose the mSUGRA benchmark 
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point BFHK B (see Tab. 12. 3p which gives, after the renormaUzation group evolution 
of the SUSY-breaking parameters performed by the SPheno computer program 
[156) . a hght en of mass = 186 GeV and rather heavy squarks with masses 
around 800-850 GeV, except for top-squark mass eigenstate ti which is shghtly 
hghter, but which nevertheless does not contribute to the virtual squark loops due 
to the negligible top-quark density in the proton. 

In Fig. \4:.16\ we integrate the previous equations with respect to M^, taking 
as lower limit the energy threshold for eRe|j production and as upper limit the 
hadronic energy, ^/sf^ = 14 TeV. We plot the LO result (dashed line) and the total 
NLL+LO matched result (solid line) with their uncertainty band relative to scale 
variation (yellow band for the LO result and green band for the NLL+LO result). 
The effect of resummation is clearly visible at small and intermediate values of qt, 
the resummation-improved result reaching a value that is even 40% higher than the 
pure fixed-order result at qx = 80 GeV . 

The asymptotic expansion of the resummation formula at LO (dotted line) is in 
good agreement with LO at small values of qt, from which we can conclude that 
the cross section is clearly dominated by the logarithms that we are resumming 
in this kinematical region. In the intermediate-gy region, we can see that the 
expansion is slightly smaller than the fixed-order calculation. Since this effect was 
not present in g^-resummation (see Fig. 14.51 in Sec. I4.2.2p . we deduce that it is 
purely related to the threshold-enhanced contributions important in the large-M 
region. It is also presented in Fig. 14.171 where we directly compare jointly- and 
gT-matched results, the main difference between the two approaches relying indeed 
in this intermediate-^T region where the jointly-resummed cross section is 5%-10% 
lower than the gy-resummed one for 50 GeV < qx < 100 GeV. 

In Fig. 14.161 we estimate the scale dependence by an independent variation of 
the factorization and renormalization scales between M/2 and 2M, and is clearly 
improved using resummation rather than pure fixed-order calculations. In the small 
and intermediate gT-region the effect of scale variation is 10% for the LO result, 
while it is always less than 5% for the NLL+LO curve. 

We also study the dependence of the total NLL+LO matched result on non- 
perturbative effects. We show the quantity A defined in Eq. ()4.39p . 

as a function of the transverse momentum of the slepton pair. The parameter A gives 
thus an estimate of the contributions from the different NP parameterizations (LY-G, 
BLNY, KN) that we included in the resummed formula, which can be found in Eqs. 
(j3.32p . (j3.33p and (j3.34p . They are under good control, since they are always less 
than 5% for qt > 5 GeV and thus considerably smaller than the resummation effects. 

The invariant-mass distribution M'^dcr/dM for e/j-pair production at the LHC 
is obtained by integrating the equations given above with respect to the transverse- 
momentum qt, and is shown in Fig. 14.181 As in the previous section, the differential 
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Figure 4.16: Transverse-momentum distribution for the process pp — > e/j e|j at the 
LHC. NLL+LO matched result, LO result, asymptotic expansion of the resummation 
formula and A-parameter are shown. 
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Figure 4.17: Transverse-momentum distribution for selectron pair production at the 
LHC ion the framework of joint and resummation. 
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cross section da/dM has been multiplied by a factor in order to remove the 
leading mass dependence of propagator and phase space factors. We can again 
see the P-wave behaviour relative to the production of scalar particles, since the 
invariant-mass distribution rises above the threshold at ^/s = 2me^ with the third 
power of the slepton velocity and peaks at about 200 GeV above threshold (both for 
M^da/dM and the not shown da/dM differential distribution), before falling off 
steeply due to the s-channel propagator and the decreasing parton luminosity. Let 
us note that we use the MRST LO 2001 |212j set of parton distribution functions 
for the LO predictions. 

In the large-M region, the resummed cross section is 30% higher than the 
leading-order cross section, but this represents only a 3% increase with respect to 
the NLO SUSY-QCD result. In the small-M region, much further then from the 
hadronic threshold, resummation effects are rather limited, inducing a modification 
of the NLO results smaller than 1%. 



The shaded, horizontally and vertically shaded bands in Fig. 14.181 represent the 
theoretical uncertainties for the LO, NLO SUSY-QCD, and the jointly-matched 
predictions. At LO, the dependence coming only from the factorization scale 
increases with the momentum-fraction x of the partons in the proton (i.e. with 
M), being thus larger in the right part of the figure, but this dependence is largely 
reduced at NLO due to the factorization of initial-state singularities in the PDFs. 
However, we add the dependence due to the renormalization scale in the coupling 
Os(/i_R), leading to a total variation of about 7%-ll%. After resummation, the 
total scale uncertainty is finally reduced to only 7%-8% for the matched result, 
the reduction being of course more important in the large-M region, where the 
resummation effects are more important. 

In Fig. 14.191 we show the cross section correction factors 

daLO/dM ^ ' 

as a function of the invariant-mass M. i labels the corrections induced by NLO 
QCD, NLO SUSY-QCD, joint- and threshold-resummation, these two last calcula- 
tions being matched with the NLO SUSY-QCD result. 

At small invariant mass M, the resummation is less important, since we are quite 
far from the hadronic threshold, as shown in the left part of the plot. At larger M, 
the logarithms become important and lead to a larger increase of the resummed K- 
factors over the fixed-order one. We also show the difference between threshold and 
joint resummations, which is only about one or two percents. It is due to the choice 
of the Sudakov form factor Q and of the W-function, which reproduces correctly 
transverse- momentum resummation in the limit of 6 — > cx), being fixed, but which 
presents some differences with the pure threshold limit 6 — > and A'^ ^ oo, as 
it was the case for joint resummation for Higgs and electroweak boson production 
[77\ I78j. However, this effect is under good control, since it is much smaller than the 
theoretical scale uncertainty of about 7%. 
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Figure 4.18: Invariant-mass distribution da/dM of en pairs at the LHC. We show 
the total NLL+NLO matched and the fixed-order NLO SUSY-QCD and LO QCD 
results, including the respective scale uncertainties as vertically hatched, horizontally 
hatched and shaded bands. 



1.36 



1.34 



1.32 



^ 1.3 



1.28 



1.26 



1.24 



pp -^y,z 



?^ Cjf* at the LHC 



■ Jointly-matched 
Threshold-matched 
NLO SUSY-QCD 
NLO QCD 




500 



1000 



1500 



2000 



2500 



3000 



M. . [GeV] 



Figure 4.19: ii'- factors as defined in Eq. (|4.75|) for en pair production at the LHC. 
We show the total NLL-I-NLO jointly- and threshold-matched result, as well as the 
fixed-order NLO SUSY-QCD and QCD results. 
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Chapter 5 



Squark and gaugino production 
and decays at hadron colliders 



In many models, the non-coloured charginos and neutralinos belong to the class 
of the lightest supersymmetric particles, their low masses counterbalancing the 
small weak cross sections for direct production. Theoretically, full NLO SUSY-QCD 
calculations have been performed ||40j, while experimentally, their subsequent decay 
to gold-plated trilepton signatures [21'i\ I214j have been exploited in several CDF 
and DO analyses at the Tevatron [2l5l Eig Em EH EH] and will be investigated 
at the LHC [22^ 1221] . Due to their strong coupling, squarks should be abundantly 
produced at hadron colliders, and hadroproduction cross sections and decay has 
therefore been studied in detail at NLO SUSY-QCD [3l]. The production of top 
[35j and bottom [95] squarks with large helicity mixing has received particular 
attention, and both QCD one-loop and electroweak tree-level contributions have 
been calculated for non-diagonal, diagonal and mixed top and bottom squark pair 
production ^96j. Very recently, flavour violation has been considered in the context 
of collider searches |97) . Concerning the associated production of squarks and 
gauginos, NLO SUSY-QCD cross sections in cMFV SUSY have been calculated 
some times ago |37 j and generalized to NMFV scenarios [97\ . 

In the following, we use for the sake of simplicity the generic notation 



for the couplings defined in Sec. 12.41 while the propagators appearing as mass- 
subtracted Mandelstam variables read 
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Figure 5.1: Tree-level Feynman diagrams for the production of charged squark- 
antisquark pairs in quark-antiquark collisions. 



5.1 NMFV squark-antisquark pair hadroproduction 

5.1.1 Analytical results 

In NMFV SUSY, the production of charged squark-antisquark pairs 

q{ha,Pa)q'{h,Pb) Ui{Pi)d*j{p2), (5.3) 

where i,j = 1, 6 label up- and down-type squark mass eigenstates, ha^b helicities, 
and Pa,b,i,2 four-momenta, proceeds from an equally charged quark-antiquark initial 
state through the tree-level Feynman diagrams shown in Fig. 15.11 The corresponding 
cross section can be written in a compact way as [97j 
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Figure 5.2: Tree-level Feynman diagrams for the production of neutral squark- 
antisquark pairs in quark-antiquark (top) and gluon-gluon collisions (bottom). 



which combine coupling constants and Dirac traces of the squared and interference 
diagrams. In cMFV SUSY, superpartners of heavy flavours can only be produced 
through the purely left-handed s-channel ly-exchange, since the f-channel diagrams 
are suppressed by the small bottom and negligible top quark densities in the 
proton, and one recovers the result in Ref. [96] . In NMFV, f-channel exchanges can, 
however, contribute to heavy-flavour final state production from light-flavour initial 
states and even become dominant, due to the strong gluino coupling. 

Neutral squark-antisquark pair production in NMFV proceeds either from equally 
neutral quark-antiquark initial states 



q{ha,Pa) q'{h,pb) qi{pi) q*{p2 



(5.6) 



through the five different gauge-boson/gaugino exchanges shown in Fig. 15.21 (top) or 
from gluon-gluon initial states 
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(5.7) 



through the purely strong couplings shown in Fig. 15.21 (bottom). The differential 
cross section for quark-antiquark scattering [97J 
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since only very few interferences (those between strong and electroweak channels of 
the same propagator type) are eliminated due to colour conservation. On the other 
hand, the gluon-initiated cross section 
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involves only the strong coupling constant and is thus quite compact. In the case of 
cMFV, our results agree with those in Ref. ^96j for diagonal and non diagonal squark 
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helicities in the final state. Diagonal production of identical squark-antisquark mass 
eigenstates is, of course, dominated by the strong quark-antiquark and gluon-gluon 
channels. Their relative importance depends on the partonic luminosity and thus 
on the type and energy of the hadron collider under consideration. Non-diagonal 
production of squarks of different helicity or flavour involves only electroweak and 
gluino-mediated quark-antiquark scattering, and the relative importance of these 
processes depends largely on the gluino mass. 

5.1.2 Numerical results 

We present numerical predictions for the production cross sections of squark- 
antisquark pairs in NMFV SUSY at the CERN LHC. Again, we use the QCD fac- 
torization theorem to compute the total unpolarized hadronic cross sections 

^ = / dr/ dy dt fa/A{xa,Ml) h,B{xb,Ml) (5.11) 

where the relevant partonic cross section da/dt is convolved with universal parton 
densities. For consistency with our leading order QCD calculation in the collinear 
approximation, where all squared quark masses, the top excepted, are much lower 
than the partonic centre-of-mass energy, we employ the LO set of the latest CTEQ6 
global parton density fit |222j . which includes nj = 5 "light" quark flavours and 
the gluon, but no top-quark density. For gluon initial states and gluon or gluino 
exchanges, the strong coupling constant asifiR.) is calculated with the corresponding 
LO value of A^q ^ = 165 MeV. We identify the renormalization scale ij-r with the 
factorization scales Ma = Mf, and set the scales to the average mass of the final 
state SUSY particles i and j, m = (mj + mj)/2. 

We take current masses and widths of the electroweak gauge bosons, electroweak 
mixing angle, electromagnetic fine structure constant and CKM matrix elements, 
while the soft SUSY-breaking masses at the electroweak scale are computed thanks 
the computer program SPheno 2.2.3 |156) . Then we introduce NMFV as in Sec. 
12.3.31 through the parameter A defined in Eq. (j2.53p . The physical masses of the 
SUSY particles and the mixing angles, taking into account NMFV, are obtained 
with the computer program FeynHiggs 2.5.1 |157j . 

The numerical cross sections for charged squark-antisquark and neutral up- and 
down-type squark-antisquark pair production are shown in Figs. 15.31 \5A\ 15.51 and 
15.61 for the benchmark scenarios A, B, C and D described in Sec. 12.3.31 respectively. 
The magnitudes of the cross sections vary from the barely visible level of 10~^ fb 
for weak production of heavy final states to large cross sections of 10^ to 10^ fb 
for the strong production of diagonal squark-antisquark pairs. Unfortunately, these 
processes, whose cross sections are largest, are practically insensitive to the flavour 
violation parameter A, as the strong gauge interaction is insensitive to quark flavours. 

Some of the subleading, non-diagonal cross sections show, however, sharp 
transitions, in particular down-type squark-antisquark production at the benchmark 
point B (bottom panel on Fig. 15. 4p . but also other squark-antisquark production 
processes. At A = 0.02, the cross sections for didg and dsd^ switch places. The 
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Figure 5.3: Cross sections for charged squark-antisquark (top) and neutral up-type 
(centre) and down-type (bottom) squark-antisquark pair production at the LHC in 
our benchmark scenario A. gg 
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Figure 5.6: Same as Fig. 15.31 for our benchmark scenario D. 
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concerned down-type squark mass differences are ratlier small (see lower right 
panel of Fig. 12. 5p . and this flip is then mainly due to the interplay between the 
different strange and bottom quark densities in the proton and the level reordering 
phenomenon discussed in Sec. 12.3.41 since at this value of A, the flavour contents of 
di and are exchanged (see Fig. I2.12p . 

At A = 0.035, the flavour content of the eigenstate ^3 switches from sr to dji, 
leading to a sharp increase of the cross sections for d^dg and ^1^3, since these final 
states can then be produced from down-type valence quarks. For A > 0.035, the 
cross section of the latter process, ^1^3 production, increases with the strange squark 
content of di (see Fig. I2.12p . Let us note that all of these non-diagonal squark pair 
production cross sections are mainly dominated by the exchange of strongly coupled 
gluinos despite of their large mass. 



5.1.3 cMFV limit 

We investigate in this subsection the production of diagonal (same mass- 
eigenstates), non-diagonal (different mass eigenstates) and mixed (different flavours) 
third generation squark-antisquark pairs within the limit of the constrained Mini- 
mal Flavour Violation, where the squark mixing is parameterized only by the two 
mixing angles 9^ and 9j^. The cMFV limit of Eq. (j5.8p in the case of diagonal squark- 
antisquark pair production, keeping only strong contributions and showing explicitly 
the mixing angle dependence, is given by [SB] 

- ^"""^ [1 - KW 



dt 9s2 ' " s2 

Ana'^ iu — mj- 

^ [1 ± {h, - K) cos 29^ - hah,] 6,, 6-^ 



2 

18^2 



•-3 

{l+hahb) {1- COS 49q)mh+{l-hahb) (3-Fcos A9q) {tu-m^^ 



tl 



in — TTl^ 
±4 cos 29g{hb-ha) 

9 



S,t6gi, (5.12) 



where the upper sign holds for 61 and the lower sign for 62 production. Stops are 
produced only through s-channel gluon exchange due to negligible top quark PDFs 
in the proton. Even for sbottom production, t-channel gluino contributions are 
suppressed by small bottom PDFs and the heavy gluino propagator. The cross 
section for the gluon-fusion initiated subprocess in Eq. ()5.10p is left unchanged, 
since it is mixing independent. In the case of no squark mixing, our results agree 
with the double-polarized and unpolarized cross sections in Ref. 



For non-diagonal production, s-channel strong diagrams are not present. How- 
ever, we take into account the PDF-suppressed gluino t-channel and the quark- 
induced tree-level electroweak diagrams, proceeding either through an s-channel Z- 
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boson exchange or a f-channel neutralino exchange, 

2 sin2(2^,-) [(1 - ha){l + h,)Ll^z + (1 + - h,)Rl^z 



dt 



2,2xl,{l-xwYsl 



X [tu 

naas sin2(20g) [(1 - + hb)Lqgz - (1 + ^a)(l - hb)Rqqz] 



9xw {1- xw) Sz tg 



[tu 



+ 



9s2 



(1 + /ia/ife)(3 + COS 46*5) - 4(/i„ + hb) cos 26*^] m?s 



+ 



2sin2 26'g(l - /ia/i6)(tw - rn|rn?J 



vra X ^ 



X 



vra \ - 



|(l + /la/i6)(3 + COS A9q 



(1 + 5ij) (1 - xvk)^ t^O t^O 
4(/ia + /ife) COS 29g^Lgg^oLgg^oR^g^oR^^^o m^om^os + (tu - m?^m, 

sin^2^,{(l-/^J(l + /.,)Lj-^-oLj-^,o + (l + /..)(l-/..)i^;-,,^ 



where we have summed over qiq2 + ^20"! final states. The couphng strengths 



-tt- -0 



(5.13) 

(5.14) 
(5.15) 



correspond the the quark-squark-neutrahno interactions in the hmit of non-mixing 
left- and right-handed (s)quarks, provided that we neglect the relatively small 
Yukawa couplings, which is the case here since top (s) quarks do not take part in the 
considered subprocess. 

Finally, mixed squark-pair production, i.e. a pair of one bottom plus one top 
squark, proceeds at tree-level only through an s-channel exchange of a charged W^- 
boson, since the t-channel exchanges of Fig. 15.11 are PDF-suppressed due to the 
negligible top density inside the proton. The ^i^i production cross section is obtained 
from Eq. (|5.4p by keeping only the W contribution. 



At 



™ IT/ |2|T/ |2 
—2- \ Vqq'\ \Vtb\ 



tu — mi mi 

ti bi 



COs2 0^ COs2 6*^ (l-/la)(l + /lfe) 

4:^2^(1 -m2,/s)2 



, (5.16) 



where we show explicitly the dependence on the CKM matrix elements. For the 
mixed production of the heavier squark mass eigenstates, the corresponding index 1 
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has to be replaced by 2 and the squared cosine of the mixing angle by the squared 
sine. 

In the following, we use older values of the SM parameters |197j , and the phys- 
ical masses of the squarks and mixing matrices are calculated using the computer 
program SuSpect Version 2.3 [202j . As a first benchmark, we choose the model line 
SPS la [in] described in Sec. 13.31 and show in the top panel of the Fig. 15.71 that 
cMFV non-diagonal and mixed third generation squark production will be difficult 
to discover at the Tevatron due to its limited centre-of-mass energy. As one can see, 
only diagonal production of the lighter top squark mass eigenstate will be visible in 
the full region of the mSUGRA parameter space shown here with the expected final 
integrated luminosity of 8.9 fb^^. For diagonal sbottom production, the accessible 
parameter space is already reduced to mi/2 ^ 225 GeV. Non-diagonal and mixed 
squark production could only be discovered, if the common fermion mass rax 12 is 
not much larger than 100 GeV. As expected for a collider, the cross sections 
are very much dominated by qq annihilation, even for the diagonal channels, and 
gg initial states contribute at most 15% in the case of diagonal light stop production. 

The LHC with its much larger centre-of-mass energy of and design luminosity of 
300 fb^^ will, in contrast, have no problem in producing all combinations of squarks 
in sufficient numbers. The hierarchy between the strong diagonal production 
channels of ©(a^) and the electroweak non-diagonal and mixed channels of 0{c?^ 
is, however, clearly visible in the bottom panel of Fig. 15.71 the latter being about 
two orders of magnitude smaller. The LHC being a "pp collider and the average 
x-value in the PDFs being considerably smaller, the diagonal channels are enhanced 
by the high gg luminosity, which dominates their cross sections by up to 93 %. 
Among the electroweak 0{a'^) processes, mixed production of top and bottom 
squarks is favoured over non-diagonal top or bottom squark production by the 
possibility of two light masses and a positive charge in the final state, which is more 
easily produced by the charged pp initial state. 

We show in Figs. 15.81 numerical results for the LHC, varying niQ independently to 
test the sensitivity of the cross section on the squark masses and mixing, starting from 
the benchmark point SPS 5 [15] , which is a mSUGRA scenario with low tan {3 = 5, 
large mi/2 = 350 GeV, and large negative Aq = —1000 GeV, leading to heavier 
sbottoms of 566 and 655 GeV, a heavy t2 of 651 GeV, but also a light ti of 259 
GeV. We see a clear hierarchy between the dominant pair production of the lighter 
stop, strong pair production of the heavier stop and sbottoms, charged and neutral 
electroweak production of final states involving at least one light squark, and finally 
charged and neutral electroweak production of the heavier squarks, which may only 
be visible up to mo < 600 GeV. The more pronounced hierarchy can be explained 
by the considerable squark mass differences in SPS 5, which lead to additional phase 
space suppression for the heavier squarks. 

5.1.4 QCD one-loop contributions for non-diagonal squark pair pro- 
duction 

Within QCD and cMFV SUSY, we have also investigated the possibility to pro- 
duce non-diagonal squark pairs by the rescattering of diagonal squark pairs through 
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Figure 5.7: Production cross sections for top (full), bottom (dashed), and mixed top 
and bottom squarks (dotted) at the Tevatron (top figure) and at the LHC (bottom 
figure) as a function of the common fermion mass mi/2 in the mSUGRA model line 
SPS la [45j. 
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pp^q. q'. * LHC, = 14 TeV 
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Figure 5.8: Production cross sections for top (full), bottom (dashed), and mixed top 
and bottom squarks (dotted) at the LHC in the mSUGRA model SPS 5 as a function 
of the scalar mass niQ . |45j . 



four-squark vertices in the final state, in the limit of a decoupled gluino. The cor- 
responding Feynman diagrams, shown in Fig. 15.91 have one-loop topology and are 
therefore suppressed by additional squark propagators and/or annihilations within 
the squark loop. The squared helicity amplitude for the production of non-diagonal 
stop pairs in gluon-gluon collisions is given by [96j 



dt 



(l+hah 



37ai sin2 (4%) , ^ ,2 , 5af cos^ (26ig) sin^ (2%) 



276487rs4 



5ai cos (26*5-) cos (2%) sin^ {29 
23047rs4 



*'* Re (A luq A In 



30727rs4 



5a.^os(2<>,)cos(2^,)s.n^(2..) . 



q,q'7^t;q=iq' 



I A In,- 1' 



(5.17) 



where 



Aln, = m? In^ (-XqJ - m? In^ (-XgJ, with x,^ - 



l+Pq 



(5.18) 



and = y 1 — 4m? is the velocity of the i*^ squark mass eigenstate. In the 
limit of degenerate light squarks, only top and bottom squark loops survive loop 



106 



Figure 5.9: Subdominant loop-level QCD Feynman diagrams for non-diagonal (i ^ j) 
squark production at hadron colliders. 

annihilations, and the squared helicity amplitude simplifies to 



(1 + hah) 



ha,hi, 

dt 



37a^in2^ „ ^ 5«f cos^ (2g-J sin^ (20,-) , ^ ^ 

(5.19) 



5a^os(2.-Jcos(2.,)sin^(2%) 

23047rs4 ^ 



These expressions have been summed over t\t2 + ^2^* final states and generalize 
the corresponding result in Ref. [35], where only top squark loops were taken into 
account. For non-diagonal sbottom production, top and bottom squark indices have 
to be exchanged in the equations above. 

For the non-diagonal third generation squark production ^1^2+^2*1 and 6162+^2^*) 
the only contributing diagrams of Fig. l5.2l is the electroweak Z-boson exchange, since 
the QCD g<^-initiated processes and the gluon s-channel are only open for diagonal 
squark pair production, and the neutralino and gluino diagrams are suppressed by the 
negligible top and small bottom densities in the proton. Indeed, for the mSUGRA 
scenario SPS 5, their contributions are found to be six to eight orders of magni- 
tude smaller than the Z-contribution. Within the SPS 5 scenario and at the LHC, 
Fig. 15.101 shows that for non-diagonal stop production, the QCD loop-contributions 
from Fig. 15.91 are smaller than the electroweak tree-level ones by about one order of 
magnitude, whereas we naively expect the 0{a^) and 0{a^) diagrams to contribute 
with roughly equal strength. As already mentioned, this can be easily explained by 
the presence of additional heavy squark propagators in the loop diagrams. Here, 
we consider not only loops involving top, but also bottom squarks, which do not 
cancel in Eq. (|5.19|) . if the masses of the two bottom squarks are unequal. However, 
the non-diagonal elements in the squark mass matrices are proportional to the rel- 
evant SM quark mass and rrif, <C mt, so that mixing effects are less important for 
sbottoms than for stops. Consequently, sbottom loops contribute about one order 
of magnitude less than stop loops, as can also be seen in Fig. 15.101 SUSY-QCD 
loop diagrams involving gluino exchanges have not been calculated here, as they 
are of 0{a^) and require in addition the presence of heavy top quark and gluino 
propagators in the loop. In the SPS 5 scenario, we have indeed a heavy gluino of 
mass rrig = 725 GeV. The dependence on the cosine of the stop mixing angle in the 
electroweak cross section of Eq. ()5.16p and in the the sbottom loop contribution in 
Eq. (I5T9I1 through sin2(2%) = 4cos2 6',- (1 — cos^ 9i) is clearly visible in Fig. 15.101 In 
contrast, the stop loop contribution in Eq. (|5.19|) has a steeper dependence through 
siv?{4,6i) = 16 cos^ 9i (1 — cos^ %) (1 — 2 cos^ Oif' , which is also visible in the figure. 
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Figure 5.10: Contributions of tree-level electroweak and loop-level QCD processes to 
non-diagonal stop production at the LHC in the mSUGRA model SPS 5 05] with a 
light top squark as a function of the cosine of the top squark mixing angle, together 
with their sums. 
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Figure 5.11: Contributions of the qq (dotted) and gg (dot-dashed) initial states 
to the longitudinal double-spin asymmetry All at RHIC, together with their sum 
(full), for stop pair production as a function of the light top squark mass. The total 
asymmetry using GRSV-VAL (dashed) instead of GRSV-STD (full) parton densities 
[201| is also shown. 
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5.1.5 Double-spin asymmetries 

As a possible application of the polarization dependence of our analytical results, 
we show in Fig. 15.111 the double-spin asymmetry of diagonal light stop production 
at RHIC (within cMFV SUSY). As this polarized pp collider has only a rather small 
centre-of-mass energy of -^/s/^ = 500 GeV, the observable stop mass range is obviously 
very limited. Already for m^-^ > 106 GeV, the unpolarized cross section drops below 
1 fb, while stop masses below 96 GeV are most likely already excluded [149j . This 
leaves only a very small mass window of 10 GeV for possible observations. In Fig. 
15.111 one can clearly see the rise of the asymmetry for qq and gg initial states, as 
the stop mass and the correlated x-value in the PDFs grows. However, as the two 
asymmetries are approximately of equal size, but opposite sign, the total observable 
asymmetry rests below the 5% level in the entire mass range shown. This is true 
for both choices of polarized parton densities, GRSV 2000 standard (STD) as well 
as valence (VAL) [201j . For consistency, the unpolarized cross sections have been 
calculated in this case using the GRV 98 parton density set |200J, as in Sec. 14.11 

5.2 NMFV squark pair production 
5.2.1 Analytical results 



Cji 




Figure 5.12: Tree-level Feynman diagrams for the production of one down-type 
squark (iji) and one up-type squark {q'j) in the collision of an up-type quark (q) 
and a down- type quark {q'). 

While squark-antisquark pairs are readily produced in pp collisions, e.g. at the 
Tevatron, from valence quarks and antiquarks, pp colliders have a larger quark-quark 
luminosity and will thus more easily lead to squark pair production. The potentially 
flavour-violating production of one down-type and one up-type squark 

q{ha,Pa)q'ihb,Pb) di{pi)uj{p2), (5.20) 

in the collision of an up-type quark q and a down-type quark q' proceeds through 
the ^-channel chargino or li-channel neutralino and gluino exchanges shown in Fig. 
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15.121 The corresponding cross section within NMFV SUSY [9; 
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involves the form factors 
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(5.22) 



where the neutralino-gluino interference term is absent due to colour conservation. 
The cross section for the charge-conjugate production of antisquarks from antiquarks 
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Figure 5.13: Tree-level Feynman diagrams for the production of two up-type or 
down- type squarks. 



can be obtained from the equations above by replacing ha^b ~^ ~ha,b- Heavy-flavour 
final states are completely absent in cMFV due to the negligible top quark and 
small bottom quark densities in the proton and can thus only be obtained in NMFV. 

The Feynman diagrams for pair production of two up- or down-type squarks with 
non-minimal flavour violation 
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(5.23) 



are shown in Fig. 15.131 In NMFV SUSY, neutralino and gluino exchanges can lead 
to identical squark flavours for different quark initial states, so that both t- and 
M-channels contribute and may interfere. The cross section [97j 
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depends therefore on the form factors 
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Gluinos will dominate over neutralino exchanges due to their strong coupling, and 
the two will only interfere in the mixed t- and u-channels due to colour conser- 
vation. At the LHC, up- type squark pair production should dominate over mixed 
up-/down-type squark production and down-type squark pair production, since the 
proton contains two valence up-quarks and only one valence down-quark. As before, 
the charge-conjugate production of antisquark pairs is obtained by making the re- 
placement ha^b — *■ —ha,b- If we neglect electroweak contributions as well as squark 
flavour and helicity mixing and sum over left- and right-handed squark states, our 
results agree with those of Ref . 01] . 



112 



5.2.2 Numerical results 



As in Sec. 15.1.21 we employ the LO set of the latest CTEQ6 global parton 
density fit [222], with five active flavours and the gluon, and the strong coupling 
constant is calculated with the value of A^q ^ = 165 MeV. The renormalization and 
factorization scales are set to the average mass of the final state SUSY particles, 
and the SUSY masses and mixings are again computed with the help of SPheno 
and FeynHiggs. 



The numerical cross sections for mixed, up- and down- type squark-squark pair 
production are shown in Figs. ISTT^ 15.151 15.161 and 15.171 for the benchmark scenarios 
A, B, C and D described in Sec. 12.3.31 respectively, where we show their dependence 
on the flavour violating parameter A (see Eq. ()2.53p ). While the diagonal channels 
are practically insensitive to flavour violation, some of the non-diagonal cross sec- 
tions show sharp transitions, as for example at the benchmark point C (Fig. I5.16P 
where they occur between the U2/un and ^2/^4 states, which are pure charm/strange 
squarks below/above A = 0.035. As a side-remark we note that an interesting per- 
spective might be the exploitation of these t-channel contributions to second- and 
third-generation squark production for the determination of heavy-quark densities 
in the proton. This requires, of course, efficient experimental techniques for heavy- 
flavour tagging. 

5.3 NMFV associated production of squarks and gaug- 
inos 



5.3.1 Analytical results 

The associated production of squarks and neutralinos or charginos 

q{ha,Pa)g{h,Pb) Xj{pi)qi{p2) 



(5.26) 



is a semi-weak process that originates from quark-gluon initial states and has both 
an s-channel quark and a t-channel squark contribution. They involve both a quark- 
squark-gaugino vertex that can in general be flavour violating. The corresponding 
Feynman diagrams can be seen in Fig. 15.181 The squark-gaugino cross section is 
given by [97] 



dt 



vr a 



-u 



(1 - ha){l - h) \Lq^qxf + (1 + ha){l + h) 



+ 



+ 



+ 



ti 



2('ut-m|m|^.; 

Stg^ 

^Xj (^Xj ~ '^Qi ) 



'1 - ha) \ Lqiqxj \ + (1 + K) \ Rg^qXJ \ 



{l-ha){l-h)\Lq^q^^\^ + {l + ha){l + h)\Rq^q^f 



Stq^ 



(1 - ha) \Lq^qXj \ + (1 + ha) \Rq,qxj \ 



(5.27) 



where = 6x14/(1 — xw) for neutralinos and = 12x\y for charginos. Note that 
the t-channel diagram involves the coupling of the gluon to scalars and does thus 
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Figure 5.14: Cross sections for mixed (top), up-type (centre) and down-type (bot- 
tom) squark-squark pair production at the LHC in our benchmark scenario A. 
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Figure 5.17: Same as Fig. 15.14] for our benchmark scenario D. 
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Figure 5.18: Tree-level Feynman diagrams for the associated production of squarks 
and gauginos. 
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Figure 5.19: Cross sections for associated production of squarks with charginos (left) 
and neutralinos (right) at the LHC in our benchmark scenario A. 



10 



sS 

1 

D 



0.7 





PP - 



























, < 1 , < < 1 < < ]" \ 



10 F 



PP -^lj<lj + C-C. 



1 r 



0.1 



0.01 





:A 








A A 






, ,,/, 1 





0.02 0.04 0.06 0.08 0.1 0.01 0.02 ^ 0.03 0.04 0.05 

Figure 5.20: Same as Fig. 15.191 for our benchmark scenario B. 
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Figure 5.21: Same as Fig. IS.liJl for our benchmark scenario C. 
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Figure 5.22: Same as Fig. 15.19) for our benchmark scenario D. 
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Figure 5.23: Tree-level Feynman diagrams for the production of gaugino pairs. 

not depend on its helicity hf,. The cross section of the charge-conjugate process can 
be obtained by taking ha — > —ha- For non- mixing squarks and gauginos, we agree 
again with the results of Ref . [31] . 

5.3.2 Numerical results 

Again, we employ the LO set of the latest CTEQ6 global parton density fit 
|222j . with five active flavours and the gluon, and the strong coupling constant 
is calculated with the value of A^q = 165 MeV. The renormalization and 
factorization scales are set to the average mass of the final state SUSY particles, and 
the SUSY masses and mixings are computed with the help of SPheno and FeynHiggs. 

The numerical cross sections for associated production of squarks with charginos 
and neutralinos production as a function of the flavour violating parameter A are 
shown in Figs. \5.19\ I5.20^ 15.211 and 15.221 for the benchmark scenarios A, B, C 
and D described in Sec. 12.3.31 respectively. The cross sections for the semi-strong 
production of very light gaugino and not too heavy squark vary from 10~^ fb to 
10^ fb and are quite sensitive to the flavour violation parameter A thanks to the 
quark-squark-gaugino vertex in diagrams of Fig. 15.181 Smooth transitions are 
observed for the associated production of third-generation squarks with charginos 
and neutralinos, and in particular for the scenarios A and B. 

For benchmark point A (Fig. 15.1^ . the cross section for ^4 production decreases 
with its strange squark content (see Fig. I2.10p . while the bottom squark content 
increases at the same time. For benchmark point B (Fig. I5.2U|) . the same (oppo- 
site) happens for dg (di), while the cross sections for uq increase/decrease with its 
charm/top squark content. Even in constrained Minimal Flavour Violation, the as- 
sociated production of stops and charginos is a particularly useful channel for SUSY 
particle spectroscopy, as can be seen from the fact that cross sections vary over 
several orders of magnitude among our four benchmark points (see also Ref. [223] ). 

5.4 NMFV gaugino pair production 

5.4.1 Analytical results 

We consider the purely electroweak production of gaugino pairs 

q{ha,Pa)q'{h,Pb) Xi{Pl)Xj{P2) (5.28) 

from quark-antiquark initial states, where flavour violation can occur via the quark- 
squark-gaugino vertices in the t- and tt-channels (see Fig. I5.23p . However, if it were 
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not for different parton density weights, summation over complete squark multiplet 
exchanges would make these channels insensitive to the exchanged squark flavour. 
Furthermore there are no final state squarks that could be experimentally tagged. 
The cross section can be expressed generically as [97] 



dt 



ira 
3^ 
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Wrl^ h.hj+MQRLQRL]{ut-m\m%) 
l^?Liil'^x.^x. + |Qii?|'*x>*x.+Re[gMii?](«i-^.^,)l, (5.29) 



i.e. in terms of generalized charges. For x^-production, these charges are given 
by 
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(5.30) 



Note that there is no interference between t- and u-channel diagrams due to 
(electromagnetic) charge conservation. The cross section for chargino pair 
production in e"^e~-collisions can be deduced by setting Cg — > = — 1, 
Lqq'z Leez = {2Tf - 2 ei Xw) and Rqq'z ReeZ = -2eixw Neglecting 
all Yukawa couplings, we can then reproduce the calculations of Ref. |224j . 
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The charges of the chargino-neutrahno associated production are given by 
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(5.31) 



The charge-conjugate process is again obtained by making the replacement 
ha,b —ha,b in Eq. (j5.29p . In the case of non- mixing squarks with neglected 
Yukawa couplings, we agree with the results of Ref. [IQ], provided we correct a sign 
in their Eq. (2) as described in Ref. |225j . 



Finally, the charges for the neutralino pair production are given by 
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Figure 5.24: Cross sections for gaugino pair production at the LHC in our benchmark 
scenario A (left) and B (right). 
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Figure 5.25: Same as Fig. 15.2^ for our benchmark scenario C (left) and D (right). 
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(5.32) 



which agrees with the results of Ref. [226| in the case of non-mixing squarks. 



5.4.2 Numerical results 

As in the previous sections, the SUSY masses and mixings are computed with 
the help of SPheno and FeynHiggs, and the LO set of the CTEQ6 parton density are 
used. The numerical cross sections for gaugino pair production are shown in Figs. 
15.241 and 15.251 for our benchmark scenario A, B, C, and D. Their cross sections are 
rather large, but practically insensitive to the flavour violation parameter A, since 
they are summed over exchanged squark flavours. 
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Table 5.1: Dominant s-, t-, and u-channel contributions to the flavour violating 
hadroproduction of third-generation squarks and/or gauginos and the competing 
dominant flavour-diagonal contributions. 
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5.5 Impact of flavour violation on squark and gaugino 
production 

In the previous sections, we have performed an extensive analysis of squark and 
gaugino hadroproduction and decays in non-minimal flavour violating supersymme- 
try. Within the super-CKM basis, we have taken into account the possible misalign- 
ment of quark and squark rotations and computed all squared helicity amplitudes for 
the production and the decay widths of squarks and gauginos in compact analytic 
form, verifying that our results agree with the literature in the case of non-mixing 
squarks whenever possible. In order to obtain numerical predictions for hadron col- 
liders, we have implemented all our results in a flexible computer program, allowing 
us to discuss in detail the dependence of the cross section on flavour violation. An 
illustrative summary of flavour violating hadroproduction cross section contributions 
for third-generation squarks and/or gauginos is presented in Tab. 15.11 together with 
the competing flavour-diagonal contributions, which are the only contributions in 
cMFV SUSY. 

5.6 NMFV decays of squarks, gluino and gauginos 

5.6.1 Squark decays 

We turn now from SUSY particle production to decay processes and show in 
Fig. 15.261 the possible decays of squarks into gauginos and quarks (top) as well as 
into electroweak gauge bosons and lighter squarks (bottom). Both processes can in 
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Figure 5.26: Tree- level Feynman diagrams for squark decays into gauginos and 
quarks (top) and into electroweak gauge bosons and lighter squarks (bottom). 

general induce flavour violation. The decay widths of the former are given by 



a 

2mlxwil- xw) ^ ^ ^' 



,2 ^2 ^2 



+ 



4 m,-o m„,. Re 



HQkXj QiQkXj 



\\l2i 2 2 2 



(5.33) 



^ \ \ 2 2 2 

- — o m^, - m_± - m / 





2 






+ 









4m-± m„' Re 



XJ- Ik' 



(5.34) 



li^mk 



2a 



3 m? 



(5.35) 



while those of the latter are given by 
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The usual Kallen function is 

A(x, y, z) = + + — 2{x y + y z + zx). 
In MFV, our results agree with those of Ref. [227]. 



A3/2(m|,m2j„m|). (5.37) 



(5.38) 
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Figure 5.28: Tree- level Feynman diagrams for gaugino decays into squarks and 
quarks (left) and into lighter gauginos and electroweak gauge bosons (centre and 
right). 



5.6.2 Gluino decays 



Heavy gluinos can decay strongly into squarks and quarks as shown in Fig. 15.271 
The corresponding decay width 
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can in general also induce flavour violation. In MFV, our result agrees again with 
the one of Ref . [227]- 



5.6.3 Gaugino decays 

Heavier gauginos can decay into squarks and quarks as shown in Fig. 15.281 (left) 
or into lighter gauginos and electroweak gauge bosons (Fig. 15.281 centre and right). 
The analytical decay widths are 
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for charginos and 
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(5.45) 



for neutralinos, respectively. Chargino decays into a slepton and a neutrino (lepton 
and sneutrino) can be deduced from the previous equations by taking the proper hm- 
its, i.e. by removing colour factors and down-type masses in the coupling definitions. 
Our results agree then with those of Ref. |228| in the limit of non- mixing sneutrinos. 



Note that the same simplifications also permit a verification of our results for squark 
decays into a gaugino and a quark in Eqs. (|5.33|) and (I5.34|) when compared to their 
leptonic counterparts in Ref. |228| . 
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Chapter 6 

Conclusion and outlook 



The Standard Model of particle physics provides a successful description of 
presently known phenomena, except for neutrino physics. However, despite of its 
success, a set of conceptual problems do not have a solution within the framework 
of the SM, such as the origin of mass, gauge coupling unification, or the hierarchy 
problem. Several attempts have been made to solve these problems, leading to 
various theories beyond the SM, even if there is still no experimental evidence of 
their existence. Each of these theories predicts new particles, with masses lying in 
the TeV-range, i.e. the discovery reach of present and future hadron colliders, the 
Tevatron and the LHC, which will then be able to put constraints on these models, 
or conclude about their (non-)existence. 

In this thesis, we have considered the production of sicptons, squarks and 
gauginos of the Minimal Supersymmetric Standard Model. Cross sections for SUSY 
particles production at hadron colliders have been extensively studied in the past 
at leading order and also at next-to-leading order of perturbative QCD, since they 
arc expected to receive important contributions from radiative corrections. These 
corrections include large logarithms, which have to be rcsummcd in order to get 
reliable perturbative results. We have thus performed a first and extensive study 
of the resummation effects for SUSY particle pair production at hadron colliders, 
focusing on Drell-Yan like slepton-pair and slepton-sneutrino associated production 
in minimal supergravity and gauge-mediated SUSY-breaking scenarios. We have 
presented accurate transverse-momentum and invariant-mass distributions, as well 
as total cross sections, resumming soft-gluon emission contributions to all orders in 
the strong coupling. 

Monte Carlo event generators are also commonly used, especially by experi- 
mentalists, to calculate observables depending on the soft-gluon emission from the 
initial state partons in hadronic collisions. This allows to compare experimental 
data to theoretical predictions and to simulate experimental signatures, when 
there are no experimental data yet. These programs usually implement the hard 
scattering process at the leading order, matching it with parton showering, the latter 
parameterizing the initial- and final-state radiation. Let us note that recently, event 
generators using next-to-leading order calculations have also been developed. It is 
expected that both Monte Carlo parton showering and resummation calculations 
should accurately describe the effects of soft-gluon emission from the incoming 
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partons. A comparison between the predictions given by the two approaches 
would certainly be useful in testing their reliability. Finally, other SUSY particle 
production processes, such as gaugino- or squark pair hadroproduction should also 
be considered. 

The quest for super symmetric particles at hadron colliders will nevertheless 
rely on our ability of predicting both the SUSY signal and the SM processes being 
the backgrounds for these searches. In this work, wc have focused on presenting 
accurate predictions for a slcpton-pair signal, i.e. two highly energetic lepton 
and a large amount of missing energy. We have not considered the background 
consisting mainly in lepton pairs coming from WW and ti decays. A detailed 
phenomenological study, including all background contributions, remains to be 
performed, in order to propose proper experimental cuts to enhance the signal over 
the background ratio. Furthermore, the sensitivity of these cuts to a complete 
experimental environment should also be investigated. 

In non-minimal supersymmetric models, novel effects of flavour violation may 
occur. In this case, the flavour structure in the squark sector cannot be directly 
deduced from the trilinear Yukawa couplings of the fermion and Higgs supermulti- 
plets. We have performed a precise numerical analysis of the experimentally allowed 
parameter space, considering minimal supcrgravity scenarios with non-minimal 
flavour violation, looking for regions allowed by low-energy, electroweak precision, 
and cosmological data. Leading order cross sections for the production of squarks 
and gauginos at hadron colliders have been implemented in a flexible computer 
program, allowing us to study in detail the dependence of these cross sections on 
flavour violation. 

A full experimental study including heavy-flavour tagging efficiencies, detector 
resolutions, and background processes would, of course, be very interesting in order 
to establish the experimental significance of NMFV. While the implementation of 
our analytical results in a general-purpose Monte Carlo generator should now be 
straight-forward, such a detailed experimental study is beyond the scope of this 
work. Other possible extensions of our work include the investigation of other 
SUSY-breaking mechanisms or the computation of the radiative corrections for all 
of these processes within the NMFV framework. 
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